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ABSTRACT 



Model Independent Explorations of Majorana Neutrino Mass Origins 



James Phearl Jenkins Jr. 



The recent observation of nonzero neutrino mass is the first concrete indication of 
physics beyond the Standard Model. Their properties, unique among the other fermions, 
leads naturally to the idea of a Majorana neutrino mass term. Despite the strong the- 
oretical prejudice toward this concept, it must be tested experimentally. This is indeed 
possible in the context of next generation experiments. Unfortunately, the scale of neu- 
trino mass generation may be too large to explore directly, but useful information may 
still be extracted from independent experimental channels. Here I survey various model 
independent probes of Majorana neutrino mass origins. A brief introduction to the con- 
cepts relevant to the analysis is followed by a discussion of the physical ranges of neutrino 
mass and mixing parameters within the context of standard and non-standard interac- 
tions. Armed with this, I move on to systematically analyze the properties of radiatively 
generated neutrino masses induced by nonrenormalizable lepton number violating effec- 
tive operators of mass dimensions five through eleven. By fitting these to the observed 
light mass scale, I extract predictions for neutrino mixing as well as neutrinoless double 
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beta decay, rare meson/tau decays and collider phenomenology. I find that many such 
models are already constrained by current data and many more will be probed in the near 
future. 1 then move on demonstrate the utility of a low scale sec saw mechanism via a 
viable 3 + 2 + 1 sterile neutrino model that satisfies all oscillation data as well as solves 
problems associated with supernova kicks and heavy element nucleosynthesis. Prom this 
I extract predictions for tritium and neutrinoless double beta decay searches. This is sup- 
plemented throughout by descriptions of practical limitations in addition to suggestions 
for future work. 
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Preface 

While there exists many astrophysical hints of physics beyond the Standard Model, 
the observation of nonzero neutrino mass is the only undisputed terrestrial evidence. It 
is set apart from the rest because it is not (directly) related to poorly understood grav- 
itational interactions. If it were just the fact that neutrinos have mass, this observation 
would be important but not earth shattering... After all, the quarks and charged lep- 
tons have nonzero mass values that span nearly six orders of magnitude (this is puzzling 
by itself). The neutrinos' mystique is enhanced because they are so different from the 
charged fermions in both structure and magnitude. Specifically, we find that the neutri- 
nos possess unprecedented tiny masses, nearly twelve orders of magnitude smaller then 
the largest quark mass, or equivalently the electroweak scale (another puzzle). Addition- 
ally, we observe large (perhaps maximal) mixing which is quite distinct from the observed 
quark mixings. It is likely not a coincidence that neutrinos are the only known electrically 
neutral fundamental fermion. They are not protected by this U{1)q gauge symmetry and 
consequently may have new types of interactions of the Majorana type. On the grander 
scale, neutrino data unambiguously confirms that the pattern of fermion masses begs for 
an explanation. Fitting the known Yukawa couplings as well as neutrino masses (of either 
Dirac or Majorana type) into the Standard Model can be easily accomphshed, but there 
is likely something more. Flavor constrains are some of the strongest in all of particle 
physics and, enhanced by information in the neutrino sector, are currently poised to select 
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the next set of fundamental laws that govern the universe. The next step on this journey 
is the exploration of neutrino mass origins. 

It is possible that this goal is beyond the reach of next generation experiments due to 
the possibility of potentially ultra high energy scales involved in neutrino mass generation. 
Still, it is my biased hope that, by combining information from distinct terrestrial, astro- 
physical and cosmological probes, the physics community can at least narrow down the set 
of model possibilities, if not choose the correct one. Of course completeness dictates that 
such analysis should be as model independent as possible. This spirit is the inspiration 
for the majority of my research interests and, in particular, the present analysis. Here, 
I embark on a model independent survey of (potentially) testable models of Majorana 
neutrino mass. 

This thesis is taken from three previously published papers done in collaboration with 
Andre de Gouvea and Nirmala Vasudevan [U, [2], [3]. They were selected and modified to 
reflect the unifying themes outlined above. While this work does not cover all possible 
models of neutrino mass, it is easily broad enough to aid in the analysis of current and 
future neutrino data. This analysis will help constrain large classes of new physics models 
as well as lead the way for other theoretical/phenomenological studies. 
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CHAPTER 1 
Introduction 

The Standard Model (SM) of particle physics has enjoyed tremendous success in re- 
cent decades. Indeed, consistency between experiments and theoretical calculations has 
been demonstrated in a variety of terrestrial channels over a wide range of energy scales 
exceeding twelve orders of magnitude!- Despite extensive searches, the majority of exper- 
iments have not yet uncovered a shred of (convincing) evidence of new physics Beyond 
the Standard Model (BSM). On the contrary, most results agree with SM predictions to 
an unprecedented and often puzzling degree. See appendix|A]for a brief review of the SM. 

Despite its success, there is still strong reason to believe that there exists new physics 
BSM. These are of both an indirect and direct nature. Indirectly, one is naturally led 
to the idea of new physics above the electroweak scale to stabilize the scalar Higgs mass 
to radiative corrections [4]. Additionally, the running and approximate intersection of 
the SM gauge couplings at high energies hint at a possible grand unification of particle 
physics within a single (broken) gauge group [5]. The details of this unification requires 
the existance of supersymmetry for full consistency. Direct evidence for phenomenon BSM 
is also available from astrophysical/cosmo logical observations. Specifically, an analysis of 
the energy budget of the universe reveals that only a small fraction of the universe is 
composed of the SM particles with the majority residing in dark energy (70%) and dark 



^This conservative range is taken from the sub-eV scale associated with atomic fine structure to the 
ultra-GeV scale accessable to present day colliders. 
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matter (26%) [6|. Finally, as will become clear shortly, the discovery of neutrino mass 
via flavor oscillations is a clear sign of physics BSM and is the motivation of this thesis. 
Originally, the evidence of this phenomenon was also of astrophysical origin but has been 
convincingly confirmed by terrestrial experiments. All of these "problems" suggest the 
existence of new physics for full resolution. 

The observation of nonzero neutrino mass is the first unambiguous evidence for physics 
BSM. In other words, neutrino mass can not be accommodated within the SM as currently 
formulated. This is easy to understand by considering the "types" of masses available to 
the neutrino. First, neutrinos may be Dirac particles and thus posses a mass term similar 
to that of the charged fermions, of the general form 

CDirac = MD{vN + Nv) (1.1) 

where is a new right-handed SM singlet dubbed "right handed neutrino." This mass 
term arises after electroweak symmetry breaking from the Yukawa coupling of the vN 
combination to the neutral component of the SM Higgs doublet. Notice that these terms 
involve two different fields, the familiar left-handed neutrino and the new right-handed 
state. The addition of is already an indication of new physics, but as is often over- 
looked in such discussions, it is not the only BSM effects guaranteed to arise from Dirac 
neutrinos. Since the singlet A^ carries no unbroken charge, it should combine to form a 
mass term N'^N, which ultimately spoils the Dirac nature of the physical neutrino state 
upon diagonalization. A new global symmetry must be imposed on the system by hand 
to forbid it whereas it is accidental within the SM. Indirectly, the Dirac nature of the 
neutrino suggests more substantial modifications to the SM in order to explain the tiny 
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Yukawa coupling constant needed to yield the observed sub-eV neutrino mass, roughly 



The second type of neutrino mass term is allowed since the neutrino carries no unbro- 
ken gauge quantum numbers. This is of the Majorana form written as 



Notice that this mass is written entirely in terms of a single field u and so no new field 
content is required for its construction. Physically, such an interaction states that left- 
handed neutrinos have a small but nonzero probability (proportional to M|,j) of being 
observed as right-handed antineutrinos. In this case, neutrinos and antineutrinos are 
equivalent. The only reason the distinction between them is convenient is that the nu- 
merical coefficient Mm is so small. This mass term, however, is not invariant under the 
full unbroken electroweak symmetry of SM SU{2)l x f/(l)y where the neutrino is charged 
under both weak isospin and hypercharge. The only means of such a construction is via 
the Higgs boson, but unlike the Dirac case, the Higgs quantum numbers do not allow the 
required renormalizable interaction. This necessarily requires the addition of new physics 
of various types. One can form the required term by coupling to a pair of Higgs doublets 
as {1 / A) H H u^u , but this is of dimension five and as such, nonrenormalizable. It requires 
an ultraviolet completion at or below the scale A to preserve unitarity. This is in fact the 
effective operator induced by the seesaw mechanisms and yields a mass of order f^/A 
after electroweak symmetry breaking. The seesaw neutrino mass mechanism will be dis- 
cussed at length in chapter |H The case of neutrino mass generation by more general high 
dimensional operators is discussed in chapter [3l Alternatively, to construct Eq. (11. 2p . one 



10-12. 
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could simply introduce a new scalar 0, in addition to the SM Higgs that transforms as a 
triplet under SU (2)^ to form the SM singlet (t)i''^v. It will yield a Majorana mass after the 
neutral component of (p acquires a vacuum expectation value either directly, or indirectly 
via a trilinear Higgs coupling (j)HH. In the first case the vacuum expectation of (j) must 
be unnaturally small to accomidate the tiny observed neutrino masses. However, in the 
latter case, provided the mass of (j) is large, the system reduces down to the previously 
mentioned dimension five seesaw operator after electroweak symmetry breaking with the 
scalar mass acting as the suppression scale A. 

It is natural to wonder about the consequences of introducing both Dirac and Majorana 
masses. It is well known that any such combination will yield physical Majorana neutrinos 
after diagonalization to the mass basis. This is the principle behind the type-I see- 
saw mechanism. Prom this it seems that the Majorana nature of the neutrino is the 
most natural choice and indeed it is a majority prejudice within the theoretical neutrino 
physics community. I emphasize that the construction of Dirac-type neutrinos requires 
the existence of a new symmetry imposed to protect lepton number (L) or more precisely 
Baryon minus Lepton Number (B-L). While B-L is accidental within the SM there is no 
reason to believe that it will be so with the addition of new physics. Ultimately, the 
nature of the neutrino is one we must answer experimentally. The overall purpose of my 
thesis is to aid in this endeavor by exploiting the connection between Majorana neutrinos, 
Lepton Number Violation (LNV), and mixing phenomena and thereby extract predictions 
for a variety of physics experiments. 
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CHAPTER 2 

The Physical Range of Neutrino Parameters 

In order to properly study Majorana neutrino interactions, it is important to under- 
stand how to fully parameterize neutrino mixing both within the SM, minimally extended 
to include neutrino mass, and in the context of more general BSM interactions. This in- 
volves questions that arise regarding both the number of free parameters and their physical 
ranges. Parameter counting arguments are well known in both the Dirac and Majorana 
cases. In words, these go as follows for the general case of n Majorana neutrinos. The 
transformation from the flavor bases, where the gauge interactions are diagonal, to the 
mass basis is accomplished by an n x n complex matrix. These 2n^ degrees of freedom 
minus ri^ unitarity conditions and n field rephasings leaves room for n{n — 1) parameters. 
For Dirac neutrinos, one is allowed to rephase both the left-handed lepton doublet and 
right handed neutrino fields separately, thus reducing the parameter count by n — 1 where 
the additional "1" arises from the universal abihty to rephase the entire system. The 
freedom is shared between real mixing angles, Majorana phases, and Dirac phases. 

While this much is known, the complex relationship between the physical ranges of 
these parameters has not been explored fully. It turns out that one may limit the range of 
a particular parameter set by simply extending the range of another in a non-trivial way. 
Furthermore, this relationship is modified in the presence of new physics. In what follows, 
I develop a framework in which to study these ranges via symmetries of the mass matrix. 
The reasoning behind this is simple: If a symmetry exists between the mixing values, it 
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implies the existence of two or more degenerate regions within the parameter space. One 
may choose to exclusively populate any one of them while at the same time taking heed 
of the consequences imposed by this choice on other parameters. This is accomplished 
using an extension of the SU{2) algebra defined in Appendix [Bl Admittedly, the use of 
this system is not optimal for a specified neutrino system. It is natural to use the well 
known algebra of SU {n) to study the symmetries of an n neutrino system. The drawback 
of this approach is that one must use a different algebra to analyze each different system. 
The utility of the approach adopted here is that it works for all cases, as will become 
apparent. 

In this chapter, I explore the physical ranges of neutrino mixing parameters within 
the SM for the minimal case of two neutrinos, the realistic case of three neutrinos, and 
the general case of n neutrinos. Special emphesis is placed on the parameter implications 
of additional sterile neutrinos and their mixing. This is followed by a discussion of the 
modifications to these relationships under the influence of new physics. 

2.1. Within the Minimally extended standard model 

With the introduction offered in Appendix |Al it is easy to see that the neutrino sector 
of the SM extended to accommodate Majorana neutrino mass may be expressed by the 
Lagrangian 

C^SM D ^ {l^'aYLW;^ + fa^iyaW^) + \ ^al^^aZ^, + l^M^p^p. (2.1) 
Y 2 ^ COS oy^ / 

Where Map is the Majorana neutrino mass matrix. This expression is written in the flavor 
bases where the charged lepton masses and neutrino interaction terms are diagonal. I now 
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determine the transformations of the neutrino mixing matrix that leaves this invariant. 
The reader should note for future use that an overall sign change of any neutrino field, 
either in the flavor or mass basis, is unphysical, in that it will not affect any observable 
process and is a common property of all quantum field theories. This followes directly 
from the rephasing freedom already discussed in the previous section. Using this I analyze 
the symmetries of the mass matrix. 



2.1.1. The Two Neutrino Case 

In the case of only two Majorana neutrinos, the mass eigenstates ui and 1^2 are related to 
the flavor eigenstates z/q, and Ujj by the unitary transformation W . This can be parame- 
terized as 

/. \ 



( \ 



(2.2) 



1 

e'^ 



\ 



J 



( 



\ 



COS Q — sin Q 
sin d cos Q 



\ 



J 



( \ 



where the phase matrix and rotation matrix R-*^^ are defined in Appendix [B] Here, all 
physics is contained in four parameters expressible as two real and positive masses mi, 1712, 
one real mixing angle 6, and one CP violating Majorana phase (p. The symmetries of 
Eq. (12.31) may be used to limit the physical ranges of these quantities. To this end, I must 
find transformations up to an overall rephasing that renders 



N^P^(</.)R^^(^)Nf = S^^P2(</)')R'^'^'' 



(2.3) 
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where Nm and Nf are unphysical shifts of the mass and flavor eigenstates respectively, and 
S"*^^ is a discrete permutation matrix that flips 1-2 vector elements. Assuming the Standard 
Model Lagrangian augmented by Majorana neutrino masses, these are just simple field 
sign redefinitions. All possible physical transformations are contained in the continuous 
shifts of mixing parameters P^(0 + 6(f))TV''^{6 + 66), as well as the discrete interchange of 
mass eigenstates. This is accomplished by the operation S"*"^ = P^(7r)R-^^(— 7r/2) which 
changes the sign of the mass squared difference Am^ = rri^ — ml. To perform a systematic 
symmetry search, I put both sides of Eq. (12.31) into the same form for easy comparison. 
Beginning on the right hand side, I find that the case without the mass eigenstate flip is 
trivially P^(0)P^((50)R'''^(56')R'^^(^^). Adding the ui 1/2 operation I commute to find 

SP^(0 + (50)R^='(^ + (5^) = P^(7r)R^2(-|)p2(0)p2(50)R^2(5^)Ri2(0) (2.4) 

= p2(7r)P^(0)P^((50)R^^(-^)R^2((50)R^^(^) 

= P^(0)P^(vr - 0)pi(0 + Sm'' {S0 - I) ^''(0) 

= P2(0)p2(7r - 20 - 50)R^2 (^66 - ^) R^^(0). 

The last step utilizes a total rephasing by e~^^'^~^^'^^ to yield the requisite form. On the 
left of Eq. (12.31) I must consider eight cases defined by the matrix pairs 



(N^, Nf) = {(I, I), (I, -I), {P\n), I), {P\n), I), {P\n), P\n)), {P\n), P^n)), 

(I,P^(7r)),(I,pi(7r))}. (2.5) 
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Matrix 


P^(<50)R'^(56») 


V^{-K -1(i>- &^)BJ-'^{&Q - 7r/2) 


P''(0)R^''(0) None 
P^(0)R^^(7r) Vc, ^ -> -Vfi 
P^(7r)R'^(0) V2 -V2 
p2(7r)R^^(7r) -vx 

P^(0)R^^(-2e) -V2,Vli^ -Vji 

P^(0)R^^(-26» + 7r) V2 ^ -V2,va. ^ -vo, 

P^(7r)R^^(-20) Vf, ^ -Vf, 
P^(7r)R^^(-2e' + tt) Vo,^ -Vo, 


(</>,e)^(0,0) 
(<^,e) ^ (<^,0 + 7r) 
(<^ + 7r,0) 
(<^,0)^(<^ + 7r,e + ^) 

(<^,e) ^ (0,-e + 7r) 


(0,e)^(-,^ + 7r,e + 7r/2) 
(0,e)-^(-,^ + 7r,e-7r/2) 
(<^,e) ^ (-0,e + 7r/2) 
(<^,0) ^ (-0,e-7r/2) 
(<^,0) ^ (-0 + 7r,-e + 7r/2) 
(0,9) ^ (-0 + ^,-9-^/2) 
(0,9) ^ (-0,-9 + 7r/2) 
(0,9) (-0,-9-7r/2) 



Table 2.1. Exhaustive summary of two neutrino mixing symmetries. See 



text for details. 



The only difficulty arises from a nontrivial Nf commuting with the rotation matrix. Since 
P-'-(7r) = — P^(7r) in this two dimensional case, it is enough to consider only 

N^P^(</.)R^2(^)P^(7r) = p2(0)N^R^2(0)P^(7r)R^2(-^)R^2(e) (2.6) 

= p2(0)N^P^(7r)Ri2(-2^)Ri2(0) 
= p2(0)N^p2(7r)Ri2^-20 + 7r)Ri2(0) (2.7) 

and change overall signs as needed, with the operation R-'-^(7r), to match the intended 
structure. Once this is done to both sides of Eq. fl2.3l) . the outer factors of P^(0) and 
R-'^^(^) cancel, revealing simplified equations that may be easily solved for the physical 
parameter shifts 50 and bQ. 

Table 12.1.11 shows the solutions to these equations for every possible case. The first 
column and row displays the simplified matrix equations obtained from manipulations 
of the left and right side of Eq. fl2.3p . respectively. For convenience, the second col- 
umn displays the unphysical field redefinitions associated with each corresponding row. 
The remaining columns contain an exhaustive list of physical transformations that leave 
the two neutrino Majorana mixing matrix invariant. All of these take the general form 
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(0, 9) (±0 + (50, ±6' + 50). Notice that many of the hsted transformations yield equiva- 
lent information. In particular, I find that the phase — > —0 for all cases where the mass 
eigenstates are exchanged since z/i ^ z/2 ^ P^(0) P"'"(0)) which must be countered by 
a total rephasing by —0. Additionally, I see that is invariant under translations by vr in 
cases defined by the number of unphysical field redefinitions. In a similar way, it is clear 
that the mixing angle 6 is only invariant under discrete sign changes as well as shifts by vr 
and 7r/2. Each invariance listed in the table may be interpreted as a single constraint on 
the physical parameter ranges. I follow [8] , choosing independent entries, and limit these 
as: 

(1) z/2 -z/2, and + vr. 

Here I find that the phase is invariant under shifts by vr, given a nonphysical 
mass eigenstate redefinition. This suggests degeneracies between the parameter 
regions of size vr. I choose G [— 7r/2, 7r/2] without loss of generality. This result 
conforms via direct symmetry arguments to the common conception of Majorana 
phase ranges. 

(2) 6 ^ 6 + -n, Ui — z/i, and z/2 — > — z/2- 

The mixing angle shift suggests parameter space degeneracies of size tt. Hence, 
is limited to span vr radians, half its original range, to the interval 9 G [— 7r/2, 7r/2], 
chosen for convenience. 

(3) 9 -9, — > -Z//3, and z/2 -z/2. 

Thus, neutrino mixing is invariant under reflections about 6* = 0, implying a 
degeneracy between positive and negative values. This suggests a halving of the 
mixing angle physical range to G [0, 7r/2]. 
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(4) 9 ^ 9 + -K V2 ^ —^2: ^ ^2 and —(f). 

This symmetry relates three distinct physical transformations, and it is clear from 
Table 12.1.11 that no smaller set will yield the same results. Hence, this relation 
may be interpreted in one of three ways. First, if 9 and are allowed to move 
within their full range as constrained above, one may choose a particular order- 
ing of mass eigenstates, or equivalently, the sign of the neutrino mass squared 
difference. This is typically done in full, two neutrino oscillation analysis where 
Am^ is taken positive. If both positive and negative Am^ are allowed, one has 
the freedom to limit either 9 G [0,7r/4] or G [0,7r/2]. In the case of standard 
vacuum oscillations, the sign of Am^ is unphysical, leading to the notion that the 
mixing parameter ranges are automatically limited as described here. While true 
for pure neutrino mixing, the SM charged current interactions shown in Eq. (12. ip . 
relevant for matter effects, break this degeneracy and one must explicitly choose 
how to limit and interpret the resulting parameter space. 

I point out that the Majorana neutrino mass matrix, which governs such lepton number 
violating processes as neutrinoless double beta decay, is invariant under these transfor- 
mations by construction. It is also easy to check that other processes such as neutrino 
flavor oscillations and neutrino-antineutrino oscillations share this property since their 
amplitudes are strongly related to the mass matrix. I illustrate some of the above obser- 
vations by example in this simple two neutrino context. Most conclusions are apparent 
by inspection here, but carry over to the more general case of n neutrino flavors in a much 
less obvious way. 
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First, consider the two flavor neutrino oscillation probability of a being measured 
as a Up after traversing a baseline L with energy E. In vacuum, this is given by 

P:™iL, E) = sin^ 26 sm' . (2.8) 

Here, the Majorana phase, as well as the sign of Am^, is unphysical. It is clear that the 
physical range of 9 is only in the [0,7r/4] interval from the factor sin^2^^. To introduce 
the Am^ sign as a physical degree of freedom, I must introduce matter effects. For 
neutrinos propagating in a constant electron density background, the mixing parameters 
of Eq. (12. 8p must be replaced by effective matter quantities. In terms of the dimensionless 
parameter A = 2\^EGpNe/ Am'', where A^'e is the local electron number density, the 
modified oscillation probability is 

Pr^ZiL, E) = sin2 ( ^^Vl + A2-2Acos20l . (2.9) 

^"^^''^ ' ^ l + A2-2Acos2^ \ AE J ^ ' 

Since cos 2^ now appears independently, I have lost the freedom to constrain 9 beyond 
[0,7r/2]. However, if I allow for a mass eigenstate flip, which induces an A ^ —A trans- 
formation, the angular degeneracy is restored. I see that one may choose to confine either 
Am? positive or G [0, 7r/4]. 

To see how the Majorana phases enter into this discussion, I must move to lepton 
number violating processes where they contribute to important physical effects. These are 
much harder to observe/constrain since all rates must be proportional to neutrino mass 
values, as opposed to the interferometric dependence in standard neutrino oscillations. 
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Consider the rate for neutrino-antineutrino oscillation which goes like [9], llOj 



sin^ 29 
4^2 



ml + ml — 2mim2 cos 



( 



2E 



- 20 



) 



(2.10) 



This is similar to the neutrino vacuum flavor oscillation formulae, except that this rate is 
directly proportional to the neutrino mass scale and shows a Majorana phase dependence. 
Here the physical range of 6 is still of size 7r/4, but now the sign of Am^ and are 
physical. I may limit the range for one of these by noting that a negative (p value can be 
compensated by mass eigenstate flip. This is the physical manifestation of the ambiguity 
noted in the final entry of the above symmetry list. The next logical step is to consider 
neutrino-antineutrino oscillations in matter, where, as in the flavor cases, the mass matrix 
diagonalization is modified by an effective matter potential. The process carries through 
as before, yielding no new information except for an illustration of the interplay between 
6, Am^ and (p at the same time. I leave this to the reader and move on to the more 
realistic three neutrino case. 

2.1.2. The Three Neutrino Case 

The case of three Majorana neutrinos is more complicated than the two neutrino analysis 
performed in Subsection 12.1.11 Nevertheless, the methodology and many features of the 
previous example carry over directly. Here, neutrino mixing is defined in terms of nine 
parameters, conventionally chosen as three mass eigenvalues {mi, m2, m^} taken real and 
positive, three real mixing angles {Ou, On, O23}, and three CP violating phases {6, 02, ^s}- 
The Majorana phases 02 and 03 are not physical, meaning that they may be phased away, 
when neutrinos are Dirac particles. The mixing between the mass eigenstates ui, 1/2 and 
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z/3 and the flavor eigenstates z/g, i/^, and Vr may be parameterized by 
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where, to conserve space, I employ the shorthand = cos% and Sij = sin 9ij. 

As in the previous case, I search for physical transformations of the mixing variables 
that leave Eq. fl2.12p invariant and interpret these as degeneracies within their parameter 
spaces, which lead directly to range limitations. Here, it is not useful to list every possible 
symmetry due to the overwhelming number of possibilities, many of which yield equivalent 
physical information. Rather, I enumerate in Table 12.1.21 a subset of simple cases that 
best reveal the underlying physics. The first column labels the transformation number 
referred to in the text, followed by the unphysical field redefinitions and the physical shifts 
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Unphysical 



Physical 
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^23 — ^23 + TI", ^13 — ^ ~^13? ^12 —012 
013 + TT, 012 012 + TT, 6^23 ^ 6*23 + VT, 5 - 



6 + 77 



Table 2.2. Summary of selected three neutrino mixing symmetries. Al- 
though this list is not exhaustive it is a complete representation of the 
symmetry structure. See text for details. 



in columns two and three, respectively. For convenience, the symmetry ordering is set 
such that the first physical transformation in each row may be used to limit a particular 
parameter provided the constraints listed above it are imposed. In this way, one may run 
down the list systematically to reveal the most constrained parameter spaces available 
to the system. One should note a few key features of this table. First, for simplicity, I 
only include those transformations that lead to parameter space degeneracies and skip 
nontrivial parameter redefinitions of the form Pj — > Pi(Pi, P2...Pn) for arbitrary mixing 
variables Pj. Such cases that mix parameters are addressed separately. Furthermore, 
the entries listed here are in no way unique, but they do constitute a complete set, 
since combining them in various ways will produce all other elements. These variations, 
however, yield no new physical insight and are therefore neglected. 

From this, the physical mixing parameters may be limited in the following way. The 
first two entries tell us that I may halve the ranges of the Majorana phases to tt with 
impunity. The next entry constrains the 6*13 mixing angle and its associated phase 5. Here, 
one may chose either variable to be limited to a range of tt provided the other occupies 
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the full 271 region. This suggests that one may have either 6 G [0, vr] and ^13 G [— vr, tt] 
or 5 G [0,27r] and 613 G [0,7r], but not both. Typically, the latter option is chosen by 
experimental analyses to maintain consistency with the two neutrino case where the Dirac 
phase, and therefore this ambiguity, is not present. It will become clear that, since these 
two symmetries are independent of unphysical transformations they are thus robust to the 
effects on non-standard interactions. Entries 4 and 5 allow for the similar range limitation 
of the mixing angles 612 and 623 to [0, vr]. Notice that while naively requiring both positive 
and negative 6^13 values, this is independent of the previously described ambiguity. This 
is because 6*13 —>■ —613 and 6 ^ 6 + n are interchangeable transformations, such that the 
mixing angle sign change may be compensated by a phase shift. Entries 6, 7 and 8 permit 
one to further limit each mixing angle to a range of 7r/2, which I take to be [— 7r/2,7r/2] 
without loss of generality. I see here that these actions utilize the nonphysical sign change 
of the mass eigenstates. Finally, making these standard choices, I see that the physical 
ranges of each mixing angle is limited to 6ij G [0, 7r/2] and the Majorana phases are limited 
to (pi G [—tt/2, 7r/2], while the Dirac phase is unconfined within its full 6 G [0, 27r] range. 

Unfortunately, without an exhaustive search, it is impossible to know if and when I 
include all relevant transformations. To this end, I must submit to plausibility arguments 
supplemented by numerical examples. The invariance criterion is simple if I consider each 
of the six parameters of the mixing matrix Eq. fl2.12p to transform as Pi ^ Pi + 6Pi 
and formally commute the result into the form Nmf/^'(p.j, 5Pj) = U'^ {Pi) M {Pi ^5 Pi). One 
must then find all possible 5Pi = 6Pi{Pj) that renders the matrix M proportional to some 
product of P''(7r)s, which are just nonphysical field redefinitions. Working through this 
manipulation, I make the observation that only factors of P'^(7r) commute through the 
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rotation matrices trivially, only inducing parameter sign changes. All other attempts do 
not commute and lead to additional terms directly proportional to some of the mixing 
parameters. It follows that M is a product of such factors, so it is natural to conclude 
that solutions 6Pi{Pj) can not exist in general for all parameter values. This statement 
holds true in the case of three parameter variations, which may be checked analytically, 
but is not obvious when all six parameters are allowed to shift. To establish confidence 
that no symmetries are missed in the general case, I perform a comprehensive numerical 
scan for solutions. Due to the large number of possible variations, this operation is 
coarse by necessity. That being said, our methodology did find the class of solutions 
already explored with no indication of other possibilities. I therefore conclude that the 
symmetries spanned by the transformations listed in Table [2?L2] are complete. 

In addition to the symmetries outlined in Table 12.1.21 physically distinct transforma- 
tions also exist that do not limit parameter spaces, but should rather be interpreted as 
variable redefinitions. These are cases where parameters mix with each other in a nontriv- 
ial way. The simplest example that still supplies physical insight is the three neutrino ana- 
log of symmetry four in the two neutrino case involving the exchange of mass eigenstates. 
Here I see that the Lagrangian is invariant under 612 6i2—tt/2, 02 — ^ ~4>2, 03 03—02, 
z/i ^ z/2 provided the unphysical sign redefinition 1/2 — > —1^2- No phase space limitations 
may be extracted from this due to the mixing of the Majorana phases 02 and 03. It 
requires a nontrivial relabeling of the 03 parameter. The utility of this transformation is 
clear when the Majorana phases are unphysical, as is the case for Dirac neutrinos, where 
the troublesome shifts may be phased away, or in all lepton number conserving processes 
where 02 and 03 simply do not show up. Under these circumstances, one may use this 
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Figure 2.1. Mapping of the inverted neutrino mass hierarchy between the 
popular (312) and (123) schemes. The relationship between mixing elements 
is clear from this plot, but the mapping between the specific mixing param- 
eters is highly nontrivial as given in the text. 



symmetry to limit parameters in much the same way as in the two neutrino scenario. 
Namely, one may choose either positive Am^g? ^12 ^ [0,vr/4]. The central position of 
the 1 — 2 parameters in this discussion is unambiguously chosen by our standard mixing 
matrix parameterizations. The general conclusion is that the i — j plane may enjoy this 
constraint if it contains the first Euler rotation of the mass eigenstates in the transforma- 
tion to the flavor basis. The noncommutative nature of orthogonal rotations renders the 
effects of any other i/j/ ^ vy much more complex. In particular, all other mass eigenstate 
exchange symmetries require mixing/relabeling of the mixing angles and Dirac phase that 
cannot be phased away in any reasonable situation. Therefore, most such transformations 
are of little practical importance and are not considered further. 
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The only exception to this involves the permutations of two eigenstate pairs ui ^ 1^2 
followed by z/2. This transformation is of interest because it defines the mapping 

between two popular parameterizations of three neutrino mixing within the context of 
the inverted mass hierarchy displayed in Figure 12.11 On one hand, as shown in the 
leftmost spectrum, one may assume the mass states ordered, from smallest to largest, 
as 171^,1711,1712 (henceforth the 312-case), so as to maintain the mixing angle definitions 
obtained for the normal hierarchy. On the other hand, one may wish to maintain the 
mass ordering scheme mi,m2,m^ (the 123-case) as shown in the rightmost spectrum 
It is natural to explore the mapping between the mixing parameters in these two 
cases. Figure 12.11 which indicates the flavor composition of each state by color coding, 
represents this transformation between the mixing matrix elements f/*j via explicit la- 



beling. From this alone I see, for example, that U, 
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jy^(i23)* implies that 



"^12 '-^23 ' + '^12 '^23 '^13 



c{2 ', where the parameter superscripts indi- 



cate the mass ordering. This is one of the simplest cases taken from the figure, but it still 
provides a very nontrivial relation between the variables. I proceed formally, using the 
language of Appendix [B], to complete the mapping. Specifically, I study the effect of the 
following operation 
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(2.12) 



on the mixing matrix. Here, one of the phase rotations P^(vr) may commute through 
the real rotations to cancel out, leaving only a simple sign change. Acting on the mixing 
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matrix, this has the following effect: 

H^^Ut = Ri=^(|)R2=^(-^)p2(02)P'(03)R''(^^i2)P'(-5)P'(5)R''(^i3) (2.13) 

xP\-5)R'^i9,s-l) 
= P1(02)P'(03)R''(^^12 + ^)P^(-5)P'(5)R^'(^^i3)P'(5) X P\-6) 

R^^(^23-|). 

Notice that the entire transformation may be absorbed into shifts of the physical mixing 
parameters leaving no nonphysical sign change requirements. Formally, the transforma- 
tions shown in Eq. fl2.14p are very simple. The largest change arises from the relabeling 
of the mixing planes which leads to a situation where the complex rotation governed by 
6*13 and 6 now acts in the 1 — 2 plane, while 612 governs the 1 — 3 rotation. As intuitively 
expected, it exchanges the roles of 612 and 613. The main difficulty arises operationally 
while massaging the Euler rotations into the conventional order, since they do not com- 
mute. To put the result of Eq. fl2.14p into the form of Eq. (12.121) one must commute 
the 1 — 2 and 1 — 3 rotations. A little thought reveals that this alone is not enough, as 
no consistent parameter mapping may exist for all relevant angles in this case. Thus, to 
proceed, one must also search for highly nontrivial transformations of 623- Commuting 
the matrices formally for the full transformation, I find the following mapping of mixing 
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parameters from the 312 to the 123 spectral cases. 
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(2.14) 

These are unique up to trivial sign changes, provided the overall nonphysical rephasing 
by e^^^~^'^^\ The transformations are relatively simple for the 9i2 and ^13 parameters, but 
the redefinition is less attractive for $23 due to its induced dependence on 5. In the case 
of Dirac CP conservation, the properties 5 = and 5 = tt are maintained in the mapping 
(as they should) and the troublesome transformations are simplified to 02 — ^ ~03 + (~7'')) 
03 ^ 02 - 03 + (-tt) and cj^ (C23S12S13 + 012523)^ / (c?3Si2 - 1). 



2.1.3. Comments on the n Neutrino case 



The results of the previous subsections may be directly generalized to the case of n 
neutrinos. Here, all physics may be parameterized by n real and positive masses rrii and 
n{n — l) mixing angles/phases. This arises from the 2n^ degrees of freedom in an arbitrary 
complex n X n matrix minus unitarity conditions and n field rephasing. Choosing the 
standard parameterizations, the remaining freedom may be categorized into n(n — l)/2 
real mixing angles (n — l)(n — 2)/2 Dirac phases % and n — 1 Majorana phases 0^. Due 
to the limited rephasing freedom, I may define select elements of the mixing matrix real. 
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Following |12j I choose to remove all Dirac phases from the three row main diagonal band 
so that rotations between the neighboring (a — 1) — a planes are taken to be real, while the 
rest are not, and thus, include a CP violating Dirac phase. Taking these complex rotations 
in the i-j plane to be R'j(%, Sij) = P'(-5y/2)pj(5y/2)R'j(%P'(5,j/2)pj(-5,j/2), I find 
that an arbitrary n x n Majorana mixing matrix may be written as: 



n r i—2 
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= p"(0„)[/Lin[^^"(^^-M 
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where the products are taken on the right such that nr=i ~ ^i^2^3---^n- The U^_i 
in the second line is actually of dimension n x n, but is suggestively written to indicate 
that the upper — 1 x n — 1 diagonal block is the reduced n — 1 neutrino mixing matrix. 
This reproduces the standard mixing matrices for the previous n = 2 and n = 3 cases 
and one may, in principle, compute this structure for an arbitrary number of states. For 
example, assuming n = 4 Eq. (I2.16P takes the form 



ul = P^(02)R''(^12)P'(03)R''(^13, (5l3)R'=^(^23)P'(04)R''(^14, 5i4)R''(^24, ^24)R''(^34) 

(2.16) 

In order for this analysis to make contact with observed reality, I now describe the 
parameterization of the 3 + n sterile neutrino scenario. Assuming the standard three light 
neutrinos that mix among themselves via the 3x3 matrix Ui and n sterile states that 
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mix via the approximately unitary n x n matrix Uh- The total mixing matrix [Utot\n+3 
to first order in the light-heavy mixing is 
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where 
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(2.18) 



Here, \UH,o]n = \UH,i\n = 1 and Qj is zero for all elements a, h except the row a = i where 
it is 

db,a+3e'''"^+\ (2.19) 

Note that all Dirac phases are present, per my outlined convention, except ^34 = (since 
I hold rotations between adjacent planes real). 

In the special 3 + 1 sterile neutrino case I can work out the mixing matrix to second 
order. 
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The schematic form of the sterile-active mixing 

(2.21) 



is clear in all cases and will be very useful in Chapter |H 

I may now search for symmetries of U\. To limit the parameter space in a systematic 
fashion, I present the symmetries in the same order as in the three neutrino case. The 
transformations outlined here are equally complete. First, notice that all Majorana phases 
may be trivially commuted to the left of the expression. As such, for all integers i G (2, n), 
a phase shift 0^ — > 0^ + vr may be absorbed by a sign redefinition of the mass eigenstate 
z/j —Vi. Hence, as was recently shown in [13j using different methods, each Majorana 
phase may be limited to a range of tt. Next, I find that Qij — — %, followed by bij — » 5jj+7r 
is also an invariance for integers i + 1 7^ j. Thus, for the complex rotations one may 
compensate for negative angle values, provided a corresponding shift vr of the Dirac phase. 
In a similar way, the angles governing real rotations participate in this invariance via the 
symmetries da-\,a -da-i,a, Oia -Oia, 9aj Saj for all integers i G (1, a - 2) and 
j G (a + l,n) provided the unphysical redefinitions z/^ — > — z/^ and ~^ ^^a^ of the 
mass and flavor states, respectively. Here, one may systematically limit the angles ^a-i,a 
since all other angles are either constrained by the previous transformations, or show 
up to the right of the real rotation in question. In this latter case one may proceed 
inductively from the rightmost rotation R"^-'-'"(^^„_i_„). Finally, I see that our system 
obeys the symmetry Oab — > Oah + tt. Oh, — >• —On,, 9aj — > —9aj for all integers i < a and 
a < j < b, provided the unphysical redefinitions Ua —Va and Vh —Uh of the mass and 



37 



flavor states, respectively. Hence, I obtain the expected result that the physical range of 
all mixing angles, Majorana phases, and Dirac phases are 9ij G [0,7r/2], G [— 7r/2, 7r/2], 
and 6ij G [0, 27r], respectively. As in the three neutrino case, no simple invariance may be 
formed with mass eigenstate permutations. 

That being said, some relations still hold among the mixing angles and phases in 
these most general cases. Such instances are highly nontrivial and a full exploration 
is beyond the scope of this analysis and not very enlightening. In any case, all of the 
examples outlined in subsections 12.1.11 and 12.1.21 are still valid in the full n neutrino 
case. To see this, I show that if Ul_i is invariant under some symmetry up to al- 
lowed state sign changes, then Ul is also invariant. In this context, invariant is taken 
to mean f/^'„i = Nmf/^_iNf where, as in subsection 12.1.11 and Nf are simply 
sign change matricies acting on the mass and flavor eigenstates, respectively. Both of 
these may be written as products of P'(vr) and handled with the algebra outlined in 
Appendix [Bl To begin, I note that Nm can always be absorbed by mass eigenstate 
sign changes, as can a vr phase shift of P"(0„). Thus, it is enough to show that an 
arbitrary P**(7r) from Nf commutes through the remaining rotations to the flavor state 
vector. Since a < n — 1, only the commutation with R*^" is nontrivial. For a < n — 1 I 
have P-(7r)R-"(0,„,5,„) = R-"(-^,„, 5,„)P-(7r) = R-"(^^,„,5,„ + 7r)P-(7r), which yields 
only a 6an shift by vr. For the case a = n — 1 1 am left with P"^-'-(7r)R"^-'-'"(^^„_i „) = 
P"-i(7r)P"-i(7r)P"(7r)R"-i'"(a„_i,„ + tt) = P"(7r)R"-i'"(0„_i,„ + vr), where the remain- 
ing phase shift is absorbed by a Un sign change after commuting through each R'", which 
induces a 5i„ ^ 5i„ + vr for all i < n. The 9n-i,n Gn-i,n + vr is absorbed by a separate 
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^ce„-i ~^ ~'^an-ij'^a^ ~^ ~^a^ transformation. Hence, I see that the invariance of f/„_i is 
preserved, provided each 5j„ occupy its full physical range of 27r. 

2.2. Nonstandard Interactions 

Many conclusions derived in the previous sections will be modified when nonstandard 
neutrino interactions are introduced into the SM Lagrangian of Eq. (12. ip . These can take 
the form of lepton number conserving interactions given by 

^nT = ilpV^^.v^Z'^ + ilp^^^,^pW'^ + h.c. (2.22) 

and lepton number violating interactions such as 

= rilpV-^vpSn + r]lpU-'JpS, + h.c. (2.23) 

Here, each line contains a neutral current type and charged current type interaction 
denoted by the superscripts on the coupling constants ^ and rj. To be concrete, these 
are written in terms of renormalizable expressions where neutrinos couple to vectors Z' 
and W , as well as charged and neutral scalars denoted by Sc and Sn, respectively. All 
that is needed for my purpose is the neutrino structure, so one may substitute a general 
fermion current in place of the bosons to yield effective operators if desired. All four of 
the couplings listed are general complex matrices in flavor space defined by the flavor 
subscripts a and (3. These are not diagonal in either the flavor or mass basis. I point out 
that terms in Eq. (12.231) violate lepton number by two units and as such, will themselves 
lead to radiatively generated Majorana neutrino masses at some order in perturbation 
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theory pLj. If this is the dominant source of neutrino masses generation, i]"'^'^^ would then 
be approximately aligned with the mass basis. 

Upon transforming to the mass basis with the unitary matrix U, I find that the 
coupling constants shift to 

C U^CU (2.24) 
rf U^rfU 

Using the invariance condition U"^ N^f/^Nf as before, I search for the symmetries of 
these interactions. First, by recalling that both N matrices are Hermitian, it is clear that 
both the lepton number violating and conserving interactions yield the same conditions, 
so it is enough to discuss the transformation of ^. For the neutral current term I want 
Nmf/^Nf.^'^Nff/Nm = U'^E^^U up to field sign redefinitions which can clearly absorb the 
Nm factors. Thus, a given symmetry of the mixing matrix is also a symmetry of neutral 
current nonstandard interactions if the resulting Nf commutes with Since Nf is 
some product of P'^'(7r), for some set of integers Oj, this translates to the condition that 
iaij = ^ja-i = for 3 1^ cii- The charged current case is even less restrictive in that all 
that is required is Nmf/"''Nf^'^ = U'^S,^ up to field sign redefinitions. Here, this reduces 
to Nf^'^ = ^'^Nf', where Nf' is some other product of P'^'(7r). In component form, this 
translates to the condition that for each hi there exists a q such that j ~ ^ 
i 7^ Q and = for all k ^ hi. This is a much looser set of conditions than for the 
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neutral current case. I point out that any symmetries of the mixing matrix broken by the 
nonstandard interactions may be restored by allowing for coupling constant sign changes 
of the form C ^ Nf^Nf or ^ Nf^^ 

In terms of the specific symmetries found in Section 12.11 I see that most parameter 
limitations still hold in the face of general new physics. These are transformations whose 
Nf may be commuted and absorbed entirely into a redefinition of the mass eigenstates. 
In particular, within the n neutrino scenario, the invariance involving 9a-i,a —da-i,a is 
broken while the following symmetries still hold: 

(1) 0i — > 0i + TT, z/j — ^ — z/j for all integers 1 <i. 

(2) 6ij — > — 5ij 5ij + vr for all integers i ^ j — 1. 

(3) 9ij — > 9ij + TT, 9i,ij^i —6*4,4+1, 9ia —9ia, 9bj ^(^bj, — > —^i, — > —Vj for all 
integers h < i and i + 1 < a < j . 

In the last entry the 6*^ j+i — transformation is factored out of the other operations, 
as it may not be independently addressed by another unbroken symmetry. Therefore, in 
the case where the positive-negative symmetry is broken, one must choose which 

degeneracy to exploit when limiting the parameter space. In either case, this leads to 
an expanded mixing angle range % G [— 7r/2, 7r/2]. This occurs when the nonstandard 
interactions, written in the flavor basis, have a structure where there exists at least one 
nondiagonal ^° entry in each of the i^^ and (z + 1)**^ row or column for the neutral current 
case. The situation is similar for charged current new terms, except that here one must 
have a nondiagonal ^'^ entry in both the i^^ and {i + 1)**^ rows, or at least one in each 
column. 
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Figure 2.2. Neutrino electron quasi-elastic scattering diagrams via the ex- 
change of massive W and Z' vector bosons. 

It is useful to explore this result with examples taken from the two neutrino scenario. 
Of course many processes may be used to illustrate these points, but for simplicity, I focus 
on neutrino electron scattering z/jC VjC. The Standard Model contribution to this via 
t-channel Z and W exchange is easy to understand and is responsible for coherent forward 
neutrino scattering in dense media that leads to matter effects in neutrino oscillations. 
For this, one may subtract the diagonal Z couplings from the effective matter potential 
as a common factor and consider only the charged current interactions. In a similar 
way, nonstandard neutrino couplings will modify matter affected neutrino oscillations. 
One may calculate the modified probabilities, but for my purpose it is enough to simply 
consider the relative amplitudes. I consider the SM contributions as well as new neutral 
and charged current interactions mediated by Z' and W , respectively. The relevant 
diagrams are shown in Figure 12.21 along with explicitly labeled vertex couplings. In this 
example, I assume that new physics has nonzero coupling to first generation charged 
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leptons, but this is not a general requirement. The amphtude for this process is 

A, = FzS.j + FwKUe, + Fz^KjCpUpi + Fw'K^Cae&pj, (2.25) 

where Fi is some function containing process-dependent kinematical factors and irrele- 
vant coupling constants. For my purpose, it is enough to know that these functions are 
independent of neutrino mixing parameters and couplings. The full expression may be 
calculated within the context of a specific model if needed. I point out that all Majorana 
phase dependency washes out of this expression since it conserves lepton number. For 
illustrative purposes, I evaluate the i = j = 1 elastic scattering case since it has contribu- 
tions from all four interaction types identifiable by their Fi factors. In component form, 
the amplitude is 

An = Fz + Fw cos' e + ^Fz^{Ce + C,, + {Ce-C,,)cos2e-{Q^ + C,e)^i^^0) 

+Fw>{CeCos9-C,,smeY. (2.26) 

I explore this term by term, beginning with the SM contributions. Clearly, the Z exchange 
terms do not effect the symmetry structure of neutrino mixing, but W mediation breaks 
the 6 degeneracy between [0, 7r/4] and [7r/2, 7r/2] due to the cos^ 6 factor. This is the reason 
behind the physical interpretation of the Am^ sign discussed by example in Subsection 
12.1.11 Moving to the new physics contrabutions I see that W mediation maintains 9 — > 
6 + 71 while breaking the 6 —9 symmetry, unless compensated by a corresponding sign 
flip of the off diagonal coupling constant ^^g. The same is true for the Z' interaction 
which breaks 6^—6 unless SJ^ + ~^ ~i.iea + ^L)- Therefore, I see explicitly that 
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flavor off diagonal nonstandard interactions expand the physical mixing angle range to 
9 e [— 7r/2, 7r/2], as expected. 

2.3. Conclusions 

I conclude with a brief summary of the key issues explored in this chapter. From 
exploiting the symmetries of the neutrino mixing matrix, I explore degeneracies within 
the neutrino parameter space. There is a significant amount of freedom in choosing the 
physical ranges of these quantities and care must be taken to ensure consistent analysis of 
neutrino data. This is particularly true when comparing the results of experiments that 
rely on different "physical" parameter combinations as in, for example, lepton number 
violating and lepton number conserving phenomena where Majorana neutrino phases 
are physical and non-physical, respectively. With this in mind, it is best to choose a 
parameter scheme and stick with it. Following the majority of the literature, I chose 
to restrict all Majorana phases 0, e [— 7r/2, 7r/2], Dirac phases Sij e [0, 27r], and mixing 
angles % e [0, tt/2] for the remainder of this thesis. Such a convention can accommodate 
all of the physics of neutrino mixing, assuming SM interactions. Interestingly, this also 
holds in the case of arbitrary BSM effects, provided that some select new physics coupling 
constants are allowed to take on both positive and negative values. 
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CHAPTER 3 

A Survey of Lepton Number Violation Via Effective Operators 

As previously emphasized, the discovery of neutrino masses via their flavor oscillations 
over long baselines constitutes the flrst solid evidence of physics beyond the standard 
model (SM) of particle physics pHj. While this is an important flrst step toward a 
deeper understanding of nature, it poses many more questions than it answers. A number 
of theoretically well-motivated models have been proposed and explored to address the 
origin of the neutrino mass but, strictly speaking, these represent only a handful out of 
an inflnite set of possibilities. The question of how well future experiments can probe and 
distinguish different scenarios arises naturally and is quite relevant given the current state 
of high energy physics. The coming years promise detailed explorations of the terascale 
with the Large Hadron Collider (LHC) and the more distant International Linear Collider 
(ILC) or variants thereof. Expectations are that combined information from these two 
facilities, coupled with high precision, low energy results and cosmological observations 
will shed light on some of the current mysteries of physics, including that of the neutrino 
mass. 

Here, I concentrate on the possibility that the neutrino masses are generated at some 
high energy scale A where U{1) b~l, the only non-anomalous global symmetry of the stan- 
dard model, is broken. Such a scenario is well-motivated by the observed properties of the 
light neutrinos including tiny masses, large mixings and the fact that neutrinos are the 
only electrically neutral fundamental fermions. More speciflcally, once U{1)b-l is broken. 
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neutrinos are not protected from getting non-zero Majorana masses after electroweak sym- 
metry breaking. On the other hand, since the renormahzable minimal standard models 
preserves U{1)b-l, B — L breaking effects will only manifest themselves at low energies 
through higher dimensional operators. This being the case, one generically expects neu- 
trino masses to be suppressed with respect to charged fermion masses by (f/A)", n > 1, 
where v is the Higgs boson vacuum expectation value. 

By further assuming that all new degrees of freedom are much heavier than the weak 
scale, I am guaranteed that, regardless of the details of the new physics sector, all phe- 
nomena below the weak scale are described by irrelevant, higher dimensional operators. 
In this spirit, the observable consequences of all high energy models that lead to small 
Majorana neutrino masses can be catalogued by understanding the consequences of irrel- 
evant operators that break B — L hj two units. With this in mind, I will survey all such 
non-renormalizable effective operators for phenomenological signatures at future and cur- 
rent experiments. I restrict myself to operators that will lead to lepton number violation 
(LNV), as these will be directly connected to the existence of small Majorana neutrino 
masses. This means that 1 do not consider operators that conserve L but violate B, and 
hence also B — L, by two units (such operators lead to, for example, neutron-antineutron 
oscillations), nor do I include operators that respect B — L. Most of the time, the latter 
will not mediate any observable consequences for large enough A, except for operators of 
dimension-six and above that can mediate proton decay. 



Throughout, I will assume that the weak scale degrees of freedoms are the known standard model fields, 
plus a minimal Higgs sector. Hence, I assume that there are no gauge singlet "right-handed neutrino" 
fermions or higher SU{2)l Higgs boson representations, such as Higgs boson triplets. 
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To begin, I systematically name and classify all relevant LNV operators. Fortunately, 
this has already been donqj in |15] up to and including operators of mass dimension 
eleven!^ For each operator I then calculate/estimate the analytic form of the radiatively 
generated neutrino mass matrix. Upon setting this expression equal to the experimen- 
tally measured neutrino masses, I extract the energy scale A associated to the new LNV 
physics. Armed with these scales, I proceed to calculate each operator's phenomeno- 
logical signatures at a variety of experimental settings. Additionally, having explicitly 
calculated the operator-induced neutrino mass matrices, I may also verify, under some 
generic assumptions, whether one can account for the observable lepton mixing pattern. 
After such a general survey, one is adequately equipped to take a step back and select 
phenomenologically/theoretically interesting operators for further detailed study by "ex- 
panding" effective vertices to reveal particular ultraviolet completions. In this way, one 
can use the results presented here as a means of systematically generating renormalizable 
models with well-defined experimental predictions. 

This chapter is organized as follows. Sec. 13.11 is devoted to an introduction to the 
effective operators and methods. In Sec. l3.1.H I derive and comment on the scales A of new 



physics that are used throughout the remainder of the text. In Sec. 13.21 I survey various 
experimental probes of LNV for each operator, and address if and when our analysis breaks 



^The authors of [15] discuss all possible effective operators of dimensions up to and including eleven, but 
only explicitly list those deemed unique in the sense that they cannot be written as the product of any 
previous operator with a Standard Model interaction. I append their list and naming scheme to include 
these into our analysis. 

'^I will argue later that irrelevant operators with mass dimension thirteen and higher, if related to neutrino 
masses, will require new physics below the electroweak scale so that I would have already observed new 
physics if neutrino masses were generated in this way. Furthermore, from a model building perspective, 
it is difhcult to develop models that predominantly yield effective operators of very high mass dimension. 
The probability that such scenarios are both theoretically well-motivated and evade all observations 
appears to be slim. 
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down due to added model structure or additional assumptions. Specifically, I study both 
current constraints and future prospects for neutrinoless double-beta decay experiments in 
Sec. 13.2. H followed, in Sec. I3.2.2[ by a similar analysis of other rare decay modes, including 
those of various mesons and W/Z gauge bosons. In Sec. 13.2.31 I present collider signatures 
of LNV as they apply to future linear collider facilities running in the e~e~ collision mode, 
and describe extensions of our analysis to include associated 77 collisions. I also comment 
on searches for LNV in future hadron machines. Sec. l3.3l describes current constraints from 
neutrino oscillation phenomenology due to the general structure of the derived neutrino 
mass matrices. In Sec. 13.41 I highlight a number of "interesting" operators, defined by low 
cutoff scales and prominence of experimental signatures, which are still allowed by current 
constraints on LNV. I undertake a slightly more detailed discussion of their characteristics 
and signatures and present some sample ultraviolet completions. I conclude in Sec. 13.51 
with a summary of our assumptions and results, augmented by commentary on future 
prospects for LNV searches. Our results are tabulated by operator name in Table 13.11 for 
easy reference. 

I hope that this analysis will prove useful to various audiences on a number of distinct 
levels. In the most superficial sense, the casual reader should note the general features 
of LNV as well as the diversity of model variations. Such information is best expressed 
in terms of the operator distribution histograms scattered throughout the text. These 
are color-coded by operator dimension or cutoff scale, and typically contain additional 
information, including current experimental prospects. On the more technical side, those 
interested in specific neutrino mass generating models will find detailed, operator specific, 
information that may be utilized as crude model predictions. Additionally, as already 
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alluded to, one may even "hand-pick" operators for model development based on specific 
phenomenological criteria. Finally, I urge experimentalists to search for new physics in 
all accessible channels. It is our ultimate goal to provide motivation for experimental 
considerations of non-standard LNV effects, beyond neutrinoless double-beta decay. 



Here I analyze SU (3)c x SU (2) l x U{1)y invariant AL = 2 non-renormalizable effective 
operators of mass dimension up to and including eleven. They are composed of only 
the SM field content as all other, presumably heavy, degrees of freedom are integrated 
out. As already emphasized, I do not allow for the existence of SM singlet states (right- 
handed neutrinos) or any other "enablers" of renormahzable neutrino masses, such as 
Higgs SU (2) L triplet states. I therefore assume that all lepton number violation originates 
from new ultraviolet physics and that neutrino masses are generated at some order in 
perturbation theory. 

A d-dimensional operator O"^ is suppressed by d — 4 powers of a mass scale A that 
characterizes the new physics, in addition to a dimensionless coupling constant A: 



where I sum over all possible flavor combinations that make up the same "operator-type," 
as defined below. For each operator, A/A is approximately the maximum energy scale 
below which the new perturbative ultraviolet physics is guaranteed to reside, and A is 
used as a hard momentum cutoff in the effective field theory. Among all d-dimensional 



3.1. The lepton number violating scale 
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operators, I define A so that the largest dimensionless couphng A is equal to unity. Unless 
otherwise noted, I will assume that all other A are of order one. 

In the first two columns of Table 13.11 I exhaustively enumerate all possible lepton 
number violating operators of mass dimension less than or equal to eleven. All together, 
this amounts to 129 different types of operators, most of which, 101 to be exact, are of 
dimension eleven and consist of six fermion and two Higgs fields. Remaining are 21, 6 and 
1 operator of dimension nine, seven, and five, respectively. The dimension-nine operators 
can be of two different kinds, as defined by their respective field content. They either 
contain four fermion and three Higgs fields or simply six fermion fields with no Higgs 
field content. For consistency, I use the notation of reference |15j , where such a listing 
was first introduced. Our operator naming scheme is also derived from the same list, 
which I trivially extend to include 21 elements only mentioned in that analysis. These 
are the dimension-nine and dimension-eleven LNV operators that can be constructed 
from the "product" of the previously listed dimension-five and dimension-seven operators 
with the SM Yukawa interactions. I individually identify those operators with the same 
field content but different SU{2)l gauge structure with an additional roman character 
subscript added onto the original designation from |l5j. This is done in order to render 
our discussion of the various operators clearer, since specific gauge structures can play an 
important role in the derived energy scale and predictions of a given operator. Note that 
I neglect effective operators that contain SM gauge fields, since, as argued in [15j . these 
are not typically generated by renormalizable models of new physics. 
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Our notation is as follows. 
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(3.2) 



are the left-handed lepton and quark SU{2)l doublets, respectively, e'^, u'^ and d!^ are the 
charge-conjugate of the SU{2)l singlet right-handed charged lepton and quark fermion 
operators, respect ivelty. Conjugate fields are denoted with the usual "bar" notation (Z, 
Q, e'^). For simplicity, I am omitting flavor indices, but it is understood that each matter 
fermion field comes in three flavors. All matter fields defined above are to be understood 
as flavor eigenstates: all SM gauge interactions, including those of the l^-boson, are 
diagonal. Without loss of generality, I will also define the L and e'^ fields so that the 
charged-lepton Yukawa interactions are flavor-diagonal. 
I take the SU{2)l doublet Higgs scalar to be 



H 



(3.3) 



and assume that, after electroweak symmetry breaking, its neutral component acquires 
a vacuum expectation value (vev) of magnitude v ~ 0.174 Tevjll thus spontaneously 



breaking the electroweak gauge symmetry SU{2)l x U{1)y f^(l)cm- In Table ISTTj the 
components of the SU{2)l doublets are explicitly listed and labeled with k, . . . = 1,2. 
In order to form gauge singlets, operators are contracted either by the antisymmetric 
tensor e^, defined such that ei2 = 1, or by trivial contractions with a conjugate doublet 



""^Our numerical value for v is distinct from many treatments of the SM where v is taken to be 0.246 TeV. 
These are equivalent up to a factor of V2 and are both valid provided a consistent treatment of the 
interaction Lagrangian. 
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field. Different gauge contractions are partially responsible for the wide variety of operator 
structures encountered in this study. 

In order to avoid unnecessarily messy expressions, several features are missing from the 
operators as listed in Table 13. 1[ To begin, SU{3)c color indices are suppressed in these 
expressions. Color contractions are only implied here because SU{3)c is an unbroken 
symmetry of the SM and hence there is no sense in distinguishing the various quark field 
components. I assume that the parent ultraviolet completion to each operator treats 
the color gauge symmetry properly by introducing appropriately chosen heavy colored 
particles to render the theory gauge invariant. Slightly more serious is the omission of 
flavor indices to label the fermion generations. For most of this analysis, I assume that 
all new physics effects are generation universal and thus, flavor independent. This is not 
guaranteed to be the case, as is painfully obvious within the SM. One will also note that, 
depending on the SU{2) structure of the effective operator, different flavor-dependent 
coefficients will be strictly related. For example, including flavor dependent couplings A^^, 
Oi should read \]^pU^L-'f^H^H^eikeji, where A^^ = A^„ (symmetric) for all a,(3 = e, /i,r. 
On the other hand, should read (for fixed Q and df' flavors) y^^pU^V^Q^ df" eije^u 
where A^^ = — A^^ (antisymmetric) for all a,(3 = e, /i,r. Large differences among the 
various flavor structures of each operator may very well exist. Flavor is an important 
facet of LNV phenomenology, and is addressed where relevant within the text. 

The final feature missing from our notation is explicit Lorentz structure. Each operator 
must, of course, form a Lorentz scalar, but there are numerous field configurations that 
can bring this about. The Higgs field is a scalar, and as such, transforms trivially under 
the Lorentz group and is thus of no relevance to this discussion. The fermions, however. 
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transform non-trivially and their contractions must be accounted for in each operator. 
Simple combinatorics dictate that there are at most 45 such possibihties for the six- 
fermion operators that comprise the bulk of our sample, 3 in the four-fermion case and 
only 1 for the lone dimension five operator. Additionally, each contraction can be made 
in a variety of ways, corresponding to the bilinear Dirac operators 1, and a^'^ — 
'^[l^iY] of the scalar, vector and tensor types, respectively. Since I am dealing with 
chiral fields, the addition of the 75 matrix to form the pseudoscalar and axial- vector 
bilinears is redundant. While this helps reduce the number of possibihties, the task of 
listing, categorizing, and analyzing all possible Lorentz structures for each operator is still 
quite overwhelming and is not undertaken in this general survey. Fortunately, different 
Lorentz structures for the same operator-type lead to the same predictions up to order 
one effects. This is especially true for the "interesting" operators characterized by TeV A 
scales. I shall quantify this statement and mention specific structures when relevant. That 
being said, the Lorentz structure of an effective operator can suggest a lot of information 
about its parent renormalizable model. For example, it can suggest the spin of the heavy 
intermediate states and the forms of various vertices. 

Armed with these operators, I can calculate the amphtude of any AL = 2 LNV process. 
It is important to emphasize that when addressing the phenomenological consequences of 
any particular operator O, I assume that it characterizes the dominant tree-level effect 
of the new heavy physics, and that all other effects - also characterized by other LNV 
effective operators of lower mass dimension - occur at higher orders in perturbation theory. 
Our approach is purely diagrammatic, in that I begin with an operator-defined vertex and 
then proceed to close loops and add SM interactions as needed to yield the correct external 
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state particles. In this sense, special care must be taken to respect the chiral structure 
as defined by each operator. In order to reach the intended external states, to couple 
to particular gauge bosons, or to close fermion loops, one must often induce a helicity 
flip with a SM mass insertion. I express these inserted fermion masses in terms of the 
respective Yukawa couphngs, yf {f — £,u, d) and the Higgs vev, v. The Higgs field can be 
incorporated into this procedure in a number of ways. I treat the two charged and single 
neutral Nambu-Goldstone Higgs bosons, and respectively, within the Feynman- 
't Hooft gauge as propagating degrees of freedom with electroweak scale masses. The 
physical neutral Higgs, /io, can be either chosen to propagate as a virtual intermediate 
state, or couple to the vacuum with amplitude v. 

In order to avoid the task of explicitly evaluating a huge number of multiloop Feynn- 
man diagrams, I succumb to approximate LNV amplitudes based on reasonable assump- 
tions and well-motivated rules. Our methodology is motivated by exact computations 
with one-loop, dimension 7 operators where the work is analytically tractable, as well as 
on general theoretical grounds. For select operators, I have also checked our assumptions 
against predictions from ultraviolet complete models with success. In order to perform a 
particular calculation, I draw the appropriate diagram(s), taking care that no momentum 
loop integral vanishes by symmetry reasons. This step is potentially quite involved, as 
multiple diagrams can give sizable amplitude contributions depending on the character- 
istic energy transfer in the system, not to mention the cumbersome Dirac algebra within 
the respective loops. Given the high, often super- TeV, mass scale associated with our 
calculations, it is often convenient to work in the gauge field basis where each boson state 
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is associated with a single SM group generator, as is natural before electroweak symme- 
try breaking. In a similar sense, all fermions, including those of the third generation, 
are taken to be masslcss to zcroth order. All masses are included perturbatively where 
needed via mass insertions. At first guess, it would seem that our results are only valid 
in the rather subjective limit v. By direct comparison with other more complete ap- 
proximations, however, I find that our predictions are very reasonable at all scales above 
0.5 TeV. Keeping all of this in mind, I apply the following "rules" to obtain approximate 
amplitude expressions. 

(1) Trivial numerical factors: A number of numerical factors can be read off triv- 
ially from the Feynman diagrams. Specifically, one can extract the presence of 
the suppression scale A"^''"^) directly from the dimension d operator, as well as 
the dimensionless coupling constants A. Generally, A is a generation dependent 
quantity, but for lack of any experimental evidence to the contrary, I take A = 1 
universally unless stated otherwise. In the case of scenarios already constrained 
by current data, I will relax this assumption to "save" the operator and comment 
on the phenomenological consequences of the change. 

Furthermore, various factors of the electroweak scale v may be extracted 
from the operator's Higgs field content, in addition to fermion/gauge boson mass 
terms. In this way, I may also include the various Yukawa and gauge coupling 
factors Uf (/ = i, u, d) and Qi, respectively, where i runs over the three SM gauge 
groups. For simplicity, I neglect the gauge subscript i in further analytic expres- 
sions. Finally, a color factor of 3 associated with each quark loop should also 
be included in our computations, but can (and will) be neglected for simplicity 
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from algebraic expressions where order one factors are irrelevant and only serve 
to render expressions more cumbersome. I note that all coupling constants are 
subject to renormalization group running. In particular, those occurring within 
a loop should be evaluated at the scale A. I neglect this order one effect since 
it is most important at large A scales where operators tend to have less of a 
phenomenological impact due to the (1/A)" suppression. 
(2) Loop factors: In all of our calculations, I assume that each operator defines 
an effective field theory, characterized by the scale A. This imphes that all 
momentum integrals are effectively cut off at A, above which new states will 
emerge to regularize the theory. Divergences in such loops tend to cancel the large 
scale suppressions inherent to the bare operators, and thus enhance predicted 
LNV rates. Specific divergences can be determined by simple power counting of 
momentum factors. Of course, multiple loop integrals are often convoluted to 
the point where substantial simplification is needed to determine the dominant 
divergent term. Such a complication is in part due to the numerator of the 
Dirac propagators, which include single momentum factors and must therefore 
be present in pairs to contribute effectively to an ultraviolet divergence. The 
process of adding loops to induce A power law divergences should only be pursued 
to the point where the suppression of the induced effective term is no less than 
A~^. Any further divergent contribution must be treated as a renormalization to 
lower order terms, and hence, can only add small finite corrections to the total 
amplitude. In any case, those diagrams with the smallest scale suppressions are 
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not always the most dominant, as will become clear later when I discuss specific 
results. 

In addition to power-law divergences, each loop is also associated with a 
numerical suppression factor. This arises from the proper normalization of the 
loop four- momentum integral as a factor of (27r)~^, the characteristic phase space 
"volume" of a quantum state. It allows one to view the integral as a coherent 
sum over all possible intermediate configurations in a consistent way. Partially 
evaluating these integrals for a number of examples, one quickly finds that two 
powers of tt cancel with the four dimensional Euclidian space solid angle J dQ^. 1 
introduce a suppression factor of (IGvr^)^^ ~ 0.0063 for each diagram loop, which 
tends to offset enhancements from associated divergent factors. A quadratically 
divergent loop diagram is often proportional to the lowest order contribution 
times (l/167r^)(A/f )^ to the power n (number of loops in the diagram). This 
contribution is larger than the leading order one if A > Anv ~ 2 TeV for any 
number of loops. The situation is often more involved, as many loops turn out 
to be logarithmically divergent or even convergent. The important conclusion 
is that adding loops is not an efficient way to enhance LNV rates at the low 
scales accessible to future experiments. This fact is demonstrated by example in 
Sec. [321 

Finally, as already alluded to, many diagram loops will exhibit logarithmic 
divergences, as is the standard case in renormalizable theories involving fermion 
and vector fields. This occurrence typically refiects the differences between the 
two characteristic scales inherent to the system, namely A and v, and are of 
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the general form ^log'^(A/t') for some power n. ^ is a small, loop suppressed, 
dimensionless coupling coefficient. Numerically, these logarithms are much softer 
than their quadratically divergent counterparts seen elsewhere in the diagrams 
and can safely be neglected. 

(3) Intermediate states: I treat all virtual intermediate states, outside of loops, as if 
they carry the characteristic momentum of the interaction Q and neglect Dirac 
structure, unless stated otherwise. In particular, goldstone bosons are assigned 
the propagator (Q^ — M|)"^ and fermions are assigned {Q — Mf)~^. In the case of 
an intermediate neutrino, this reduces to a simple factor oiQ~^ for all realistic Q 
values. Hence, for very low energy processes {Q <^ 100 MeV), neutrino exchange 
diagrams tend to dominate LNV rates. 

(4) Lorentz structure: For the purposes of our analysis, I assume that all Lorentz 
contractions between fermions are scalar-like. As previously mentioned, the ab- 
solute magnitude of most LNV amplitudes is robust under this assumption up 
to order one factors. The only qualitative exception to this occurs in some cases 
involving fermion bilinear terms with a tensor Lorentz structure {i)' (J n^ip) . This 
factor, when coupled between two fermions contracted in a loop, will yield a van- 
ishing rate due to its antisymmetry inside of a trace, since Tr{(T^y) — 0. This 
can be bypassed by introducing a new momentum vector into the trace, implying 
the addition of another loop. In most cases, this is most efficiently accomplished 
with a new gauge boson line, which is accompanied by a logarithmic divergence. 
The combination of both factors leads to a marginal amplitude suppression (with 
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respect to the same operator where all fermion bilinears are Lorentz scalars) for 
all energies of interest. 

With these approximations in hand, it is a simple matter to estimate the amplitude 
associated to any given diagram. Still, one must wonder about the uncertainty induced 
onto the calculations by such varied assumptions. Can results obtained by such methods 
supply valid physical predictions? The answer, of course, depends on the question that 
is being asked. Here, I will only be interested in estimating order of magnitude effects, 
including what value of A is required in order to explain the observed neutrino masses 
and, once A is so constrained, what is the order of magnitude of other related observable 
effects. 

One may wonder whether a more detailed estimate of the effects of each individual 
operator would lead to more reliable results. The answer is negative. It is easy to show 
that different renormalizable theories that lead to the same effective operator at tree- 
level will mediate different processes at the loop-level with order one different relative 
strengths. Furthermore, the derived cutoff scales inherit the uncertainty from the absolute 
value of the heaviest neutrino mass, which is only loosely bounded between 0.05 eV and 
1 eV by the extracted atmospheric mass squared difference |16] and tritium beta decay 
kinematic measurements |17L I18j . This is an order of magnitude uncertainty that cannot 
be avoided even if one were to perform a detailed computation within a well-defined 
ultraviolet complete theory. 

In summary, given all approximations and uncertainties, our results are only valid up 
to ± an order of magnitude. In this spirit, one need not explicitly consider order one 
factors that will necessarily yield negligible corrections by these standards. Such a large 
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error tolerance supplies the need for care when interpreting results. In particular, one 
should not place too much emphasis on any one bound or prediction, unless it is very 
robust, I.e., able to withstand variations of at least a factor of ten. Of course, for those 
operators constrained by several different independent sources one can, and should, take 
more marginal results seriously. 

3.1.1. Neutrino Masses and the Scale of New Physics 

Having defined the set of LNV operators, I now extract the scale of new physics from the 
direct comparison of radiatively generated neutrino mass expressions to their observed 
values. Since there are three light neutrino masses, I will use the heaviest of these to set 
the overall mass scale. Neutrino oscillation data, currently providing the only evidence 
for neutrino masses, constrain the relative magnitudes of the mass eigenstates but not the 
overall scale |14j . Such data only supply a lower bound on the heaviest neutrino mass, de- 
rived from the largest observed mass squared difference Am^3 = \ml — mW ^ 0.0025 eV^, 
the atmospheric mass squared difference |16] . At least one neutrino mass must be greater 
than Am^3 ^ 0.05 eV. Neutrino oscillations also teach us that the next-to-heaviest 
neutrino weighs at least a/Attt^ ^ 0.009 eV (the solar mass-squared difference), in such 
a way that the ratio of the heaviest to the next-to-heaviest neutrino masses is guaranteed 
to be larger than, approximately, 0.2. No lower bounds can be placed on the lightest 
neutrino mass. An upper bound on the heaviest neutrino mass is provided by several 
non-oscillation neutrino probes. Cosmology provides interesting constraints on the sum 
of light neutrino masses, but these are quite dependent on unconfirmed details of the 
thermal history of the universe and its composition |19L I20L I21L 122] . Most direct are 
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kinematic measurements of the tritium beta decay electron endpoint spectrum [23j . Both 
types of probes provide upper bounds near 1 eV, hkely to improve in coming years. I 
choose to perform our calculations assuming the mass scale my ~ 0.05 eV, corresponding 
to the experimental lower bound. In this way, each extracted operator scale A, inversely 
related to the neutrino mass, represents a loose upper bound. Since most rates for LNV 
observables are proportional to some inverse power of A, this choice implies the added in- 
terpretation that, all else remaining equal, our results for such rates should conservatively 
reflect lower limit predictions. 

LNV neutrino masses are nothing more than self-energy diagrams evaluated at van- 
ishing momentum transfer. These must couple together the left-handed neutrino state 

with the right-handed anti- neutrino state vp, as shown schematically in diagram (a) 
of Fig. 13. 1[ Here the flavor indices a and (3 can accommodate any of the three lepton 
flavors = e, /i,r). The derived Majorana masses rriap = are generally complex. 
The large grey circle in this diagram represents all possible contributions to the neutrino 
mass. Specifically, it contains the underlying AL = 2 operator along with all modifica- 
tions needed to yield the correct external state structure. This includes such objects as 
loops, additional gauge boson propagators and SM coupling constants. Generally, several 
diagrams can contribute to this mass generation, but special care must be taken that these 
are not proportional to any positive power of Q, the momentum carried by the neutrino 
legs, as this would not lead to a nonzero rest mass correction. 

Diagrams (6) — (e) of Fig. [STTl are examples that serve to illustrate some typical features 
encountered in our effective operator induced self-energy calculations. The underlying 
LNV operators shown in each diagram contain six fermion fields and are therefore of 
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id) (e) 

Figure 3.1. Sample diagrams that radiatively generate Majorana neutrino 
masses. Diagram (a) is representative of all operators that can generate the 
needed external state neutrinos. This usually proceeds via loop contrac- 
tions and other couplings, hidden within the light gray region. Diagrams 
(6) — (e) help illustrate the methodology of this analysis. Despite obvi- 
ous differences, all of these generate effective dimension-five interactions of 
calculable strength. See the text for more details. 



dimension nine or eleven, as are the majority of the analyzed operators. Each of these 
diagrams generates an effective dimension-five interaction 

_ (5) (L"/f)(L/^g) 

^5 -tap ' ^ ' 

(5) (5) 

where — is a generation dependent coupling constant that is calculable, given 
the structure of the original operator. This is easily verifiable via direct power counting, 
despite differences in dimension, loop number, field content, and helicity structure. It 
turns out that most operators, especially those characterized by super-TeV scales, possess 
this property. 

I describe each sample diagram in turn to point out important features. A subset of 
the subtleties described below is encountered when estimating the neutrino masses niap 
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for the entire effective operator set. Diagram (6) is a simple two-loop radiatively generated 
mass term proceeding from dimension-nine operators, such as Oi^ — ULPQ^d^Q^dFeikeji, 
containing the fermion structure fLfRfLfiii where the fields / and /' arc contracted into 
loops with mass insertions that supply the needed field chirality flip. Masses arising from 
such operators are proportional to two powers of fermion Yukawa couplings. Strictly 
speaking, allowed fermions from all three generations traverse the closed loops and con- 
tribute to the mass. However, assuming universal new physics coupling constants, third 
generation fermions will strongly dominate the induced neutrino mass. In cases such as 
these, one can freely suppress couplings to the lighter two generations without modifying 
the expected value of A. Since diagrams arising from the majority of our operator set 
contain at least one loop of this kind, this property proves quite useful when attempting 
to avoid low energy nuclear physics constraints, as will be discussed in more detail later. 

Diagram (rf) involves an operator of dimension eleven, such as 
O22 — ULP L^e'^Lke'^H'-H'^eiiejm-i but has a similar structure to Diagram (&) since both 
neutral Higgs flelds Hq couple to the vacuum, yielding a factor. In this case, the parent 
underlying operators contain the fermion structure fL{R)fL(R)f'L(R)fL{R)i simple variants 
thereof. Prom this, it is clear that such operators will create and annihilate the same field 
and one can close the fermion loops without mass insertions. A little thought reveals that 
such loops, if left on their own, will vanish by symmetry, since J d'^klk^ / gik"^)] — for 
all functions g{k^)- Hence, non-zero neutrino masses only appear at a higher order in 
perturbation theory [i.e., I need to add another loop). To maintain the chiral structure of 
the diagram a gauge boson line insertion is always the most effective. The specific gauge 
field required in this step depends critically on the quantum numbers of the fermions / 
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and /' contained in the operator itself and must be determined on a case-by-case basis. 
The absence of Yukawa dependence renders the estimated value of the cutoff scale A 
insensitive to the values of the dimensionless operator couplings (A), given the way A 
is defined. Notice that this three- loop diagram, like Diagram (6), predicts an anarchic 
neutrino Majorana mass matrix, currently allowed by the neutrino oscillation data |24j . 
That is, up to order one corrections, all entries are of the same magnitude, m,^ ^ 0.05 eV. 
This is in contrast to the remaining sample diagrams ((c) and (e)), which both suggest 
flavor-structured mass matrices. 

Dimension-nine operators, such as Oig = UQ^ d'^df'e^u'^eij, yielding diagram (c) have 
the peculiar property that, upon expanding out the various SU{2)l contractions in terms 
of component fields, no Vai-'p content is present to form the external legs of a mass di- 
agram. Here the LNV is introduced via the fermion structure VaGRp-, which annihilates 
a left-handed neutrino and creates a left-handed positron. Hence, to tie in the needed 
antineutrino line, one must both flip the charge lepton helicity and carry away the excess 
charge with some bosonic state. Of course, such a charged boson is guaranteed by charge 
conservation to be needed elsewhere in the system to close some ///' loop. In this par- 
ticular example, the process is illustrated by the exchange of a charged Higgs goldstone 
boson to clearly visualize the chiral structure, but of course, one may equivalently 
think in terms of a transverse W~ boson. The crucial point is that this mass is necessarily 
proportional to a charged lepton Yukawa coupling yi^, of the same flavor as the external 
neutrino, since I am working in the weak eigenbasis where the gauge couplings and the 
charged-lepton Yukawa couplings are flavor diagonal. By symmetry, the contribution to 
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the ma/B entry of the neutrino mass matrix is proportional to Ve^+yep, which reduces to the 
largest coupling y^^ in the realistic case of hierarchial charged lepton Yukawa couplings. 

Finally, Diagram (e) yields a five-loop suppressed neutrino self-energy originating from 
a dimension-eleven LNV operator such as O^q = e'^e'^Q'^d'^Q^ df^H^H^eikeji. This represents 
the most complicated structure considered in this analysis. As in Diagram (c), no explicit 
Va^js structure is available in the underlying operator, but in this case all of the LNV 
arises from e^„e^^-type interactions. Curiously, this interaction already flips helicity as 
it annihilates a left-handed positron and creates a right-handed electron. Unfortunately, 
being an SU{2)l singlet, cr only couples to the neutrino via a charged Higgs induced 
Yukawa interaction; therefore, this amplitude must be proportional to the product y^^yi^ 
to yield a legitimate neutrino mass contribution. One might imagine that the Higgs fields 
contained in the LNV operator could be used to produce the needed neutrino legs, but this 
is not possible since the resulting loop would have a structure of the form J d'^^^j^y 
which vanishes in the "rest mass" limit. It is clear that both Higgs fields must again couple 
to the vacuum and the needed flip must come from the other fermion loops. The resulting 
loop integrals can be separated into two convoluted pieces corresponding to both loop/leg 
pairs. A little thought reveals that each loop set contains three fermion lines whose 
associated integral is again proportional to the momentum of the external neutrino, and 
thus is not a valid mass correction. To fix this last problem without further complicating 
the chiral structure, one can add a gauge boson exchange between the fermion loops, as 
was also done in Diagram {d). 

Despite the dominance of the generated dimension-five interactions described by Eq. fl3.4p 
for the majority of the studied LNV operators, I find that this need not be the case for all of 
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them. For some operators, the dimension-five neutrino-mass effective operator Eq. (13.41) 
occurs at higher order in perturbation theory than the dimension-seven neutrino-mass 
effective operator (schematically, {LH)^H'^). For these, the neutrino mass matrix is gen- 
erated after electroweak symmetry breaking from 

_ iLH)iLH) , ^^r) iLH)iLH)HH 

- + , (3.5) 

(7) 

where are new calculable coefficients. This type of structure is present in the following 
operators 

O7, 021^ ,,, O22, O23, O25, C26(,, 027^, O29,C'30,, O31, O44,, C'sT- (3.6) 

In general, they are associated with dimension-eleven operatorgj whose mass diagrams 
are found trivially by connecting the external fermion loops and coupling the neutral 
Higgs fields to the vacuum. This adds two factors of the electroweak scale to the mass 
expressions. Dimensional analysis dictates that the fermion loops must conspire to yield 
an addition factor of f ^, usually from mass insertions utilized to fiip helicities. For the 
dimension-seven operators in Eq. (13. 5p . the resulting neutrino mass expression is pro- 
portional to f^/A^. If I assume, as is usually the case, that most of the dimensionless 
factors of Eq. (13.51) are common to both C.^^ and ^.^p, I find 771^3 oc l/167r^ + t>^/A^. In 
such cases, the dimension-seven contribution is only relevant for operator cutoff scales 
A < Attv ^ 2 TeV. Such low scales are seldom reached considering that these opera- 
tors are efficient at mass generation at low orders and consequently do not possess the 
necessary suppression factors. Still, for completeness, I include these terms when relevant. 



^This also occurs with operator O7, which is of dimension nine. This operator is the exception, in that 
it exphcitly contains three Higgs bosons which naturally aids in building the needed factors in a way 
similar to that discussed in the text. 



Table 3.1: Dimension-five through dimension-eleven LNV operators analyzed in this survey. The first two columns display the operator 
name and field structure, respectively. Column three presents the induced neutrino mass expressions, followed by the inferred scale of 
new physics, Ky. Column five lists favorable modes of experimental exploration. Column six describes an operator's current status 
according to the key U (Unconstrained), C (Constrained) and D (Disfavored). See text for details. 
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The third column of Table I3.H labeled rriaf^, presents our estimate for the operator- 
induced Majorana neutrino mass expressions. These were derived based on the estimation 
procedure discussed earlier. Trivial order one factors, as well as the generation dependent 
coupling constants A have been omitted, as already advertised. Flavor specific charged 
lepton Yukawa couplings are explicitly denoted yg^ and to distinguish them from yg, 
yu and yd, meant to represent a, /3-independent Yukawa couplings. A summation over all 
"internal flavors" is assumed for each entry. For order one coupling constants, this sum 
is strongly dominated by third generation Yukawa couplings. Upon setting these mass 
expressions equal to the observed scale of light neutrino masses (0.05 eV), I extract the 
required cutoff scale A for each operator. This quantity, defined to be A^,, is listed in 
column four in units of one TeV. Numerical results were obtained assuming the current 
best fit values for all SM parameters. Associated errors are negligibly small as far as our 
aspirations are concerned. 

Fig. 13.21 displays the distribution of extracted cutoff scales, K^. The histogram bars 
are color coded to reflect the different operator mass dimensions. The distribution spans 
thirteen orders of magnitude, from the electroweak scale to 10^^ TeV. It is interesting to 
note the general trend of operator dimension with scale: as expected, higher dimension 
operators are characterized by lower ultraviolet scales. For operators associated with the 
lowest ultraviolet cutoffs, the lepton number breaking physics occurs at the same energy 
scale as electroweak symmetry breaking. In this case, one needs to revisit some of the 
assumptions that go into obtaining the bounds and predictions discussed here. Regardless, 
it is fair to say that some of these effective operators should be severely constrained by 
other experimental probes, as will be discussed in the next section. 
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Figure 3.2. A summary histogram of the scale of new physics extracted 
from the 129 LNV operators introduced in Table 13. 1[ I assume a radia- 
tively generated neutrino mass of 0.05 eV and universal order one coupling 
constants. The contributions of operators of different mass dimensions are 
associated to different colors (shades of gray), as indicated in the caption. 



The natural scale for most of the explored operators is well above 10 TeV, and thus 
outside the reach of future experimental efforts except, perhaps, those looking for neutri- 
noless double-beta decay. The remainder, however, should yield observable consequences 
in next-generation experiments. This small subset arguably contains the most interesting 
cases on purely economic grounds, as they naturally predict tiny neutrino masses as well 
as TeV scale new physics, which is already thought to exist for independent reasons. It 
is aesthetically pleasing to imagine that all, or at least most, of nature's current puzzles 
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can arise from the same source, as opposed to postulating various solutions at differ- 
ent energy scales. It is important to note that one can "push" more of the operators 
into the observable TeV window by modifying the coupling of the new physics to differ- 
ent fermion generations. In particular, since many of the induced neutrino masses depend 
upon fermion Yukawa couplings, one can efficiently reduce scales by simply and uniformly 
decoupling the third generation. In most cases this can yield a Ky reduction of several 
orders of magnitude; a factor that can be further enhanced by also decoupling the sec- 
ond generation. Under these conditions, the resulting distribution, analogous to Fig. 13. 2[ 
would show the majority of the operators piled up near and slightly above the electroweak 
scale. A detailed exploration of this possibility would be impractical and is not pursued 
further. I will, however, like to emphasize that this strategy of decoupling the new physics 
from the heavy fermions is very non-standard. In most cases, one is tempted to decouple 
light fermions from new physics both because these lead to the strongest constraints and 
because one tends to believe that the large Yukawa couplings of third generation fermions 
are entangled with the physics of electroweak symmetry breaking. 

Not all extracted cutoff scales are subject to a strong dependency on SM Yukawa 
couplings. In particular, the Aj, values for the majority of dimension-eleven operators in 
the large histogram bar near 10^ TeV would not shift down at all under this hypothet- 
ical decoupling of the third generation from the new physics. These are the operators, 
as shown in Diagram {d) of Fig. 13.11 whose induced neutrino mass matrix is indepen- 
dent of the Yukawa sector. In such cases, m^^ are only functions of the various gauge 
couplings. As such, these constitute the most robust results of our analysis. These op- 
erators all predict an anarchic Majorana neutrino mass matrix of overall scale given by 
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TTLj^ = (7^/(167r^)^f ^/A, implying an energy scale A^, ~ 10^ TeV. The only other "Yukawa 
invariant" cutoff scale estimate arises for the dimension- five operator Oi. Oi captures 
the physics of all versions of the seesaw mechanism [Tj, and is at the heart of most of 
the model building currently done within the neutrino sector. Its ultraviolet completion 
can precede in only three distinct ways [25] . These possibilities are via the exchange 
of heavy gauge singlet fermions (type I seesaw), SU{2)l triplet scalars (type II seesaw) 
[26], SU{2)l triplet fermions (type III seesaw) |27] . or some combination thereof. Its 
popularity is well-founded for a number of reasons, including its underlying simplicity 
in structure as well as the purely empirical fact that it is the "lowest order means" of 
neutrino mass, and as such is easily generated by a "generic" LNV model. Additionally, 
the high scale associated with the seesaw mechanism can be easily incorporated within 
existing theoretical models and serves to help explain the observed baryon antisymmetry 
of the universe via leptogenesis [28] . 

For the purposes of direct observation, Ois high cutoff scale, nearly 10^^ TeV, places 
it well outside of the "detectable region" (A < 10 TeV) and renders it uninteresting for 
the purposes of our analysis. Of course, there always remains the possibility that Oi is 
generate by very weakly coupled new physics (or very finely-tuned new physics [29] ). in 
which case I expect to run into the new ultraviolet degrees of freedom at energies well 
below 10^^ TeV. In the case of Oi, it has been argued that new physics at almost any 
energy scale (from well below the sub-eV realm to well above the weak scale) will lead to 
light neutrino masses [301 131L [2|, 132] without contradicting current experimental results. 
Such possibilities - related to the fact that the new physics is very weakly coupled - are 
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not being explored here, as I always assume that the new degrees of freedom are heavier 
than typical experimentally accessible energy scales. 

Armed with our derived new physics scales Aj^, 1 proceed to plug them back into the 
different irrelevant LNV operators and search for possible means of future observation 
as well as already existing constraints. Generally, those operators that yield the largest 
experimental signals have the lowest cutoff scales. I conclude that, if associated to neutrino 
masses, the effective cutoff scale A,^ of the following effective operators is constrained to 
be less than 1 TeV: 

C34, O3Q, O37, C38, C53, ^54„,(,,e,d) ^56, Oqq, Oqq, Oro- (3.7) 

These may lead to observable effects at future high energy accelerator facilities. Addition- 
ally, such low scales may also indirectly lead to observable effects in "low energy" (but 
high sensitivity) experiments. There are more operators associated with slightly higher 
scales between (1 — 10) TeV that may manifest themselves experimentally via virtual 
effects. These are 

C^16, Ois, OiQ, O35, 043a,i,,e' ^50, C>52, ^55a,b,c^ ^57, ^^58, ^^65, ^68„,i,, ^73a) ^75- (3.8) 

These operators yield finite predictions for more than one observable, such that experi- 
mental efforts in seemingly unrelated fields can help constrain the class of possible LNV 
models or even help identify the true LNV model. 
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3.2. General operator constraints and predictions 



There are, currently, bounds on LNV processes from a number of independent experi- 
mental sources |33l 134] . Many of these are presently too mild to constrain the operators 
listed in Table 13.11 once their ultraviolet cutoffs A are set to the required value indicated 
by the presence of non-zero neutrino masses, Kv The situation, however, is expected 
to improve in the next several years with increased rare decay sensitivities and higher 
collider energies. Here I survey the experimental signatures of these operators in terms of 
the minimal scenarios described above. Specifically, I address the potential of neutrinoless 
double-beta decay (Sec. I3.2.ip . rare meson decays (Sec. I3.2.2p . and collider experiments 
(Sec. 13.2.31) to constrain the effective operators in question, assuming that, indeed, they 
are responsible for the observed non-zero neutrino masses. As before, I will use the ap- 
proximations discussed in Sec. 13. and warn readers that all the results presented are to 
be understood as order of magnitude estimates. The results, however, are useful as far as 
recognizing the most promising LNV probes and identifying different scenarios that may 
be probed by combinations of different LNV searches. 

Most of this section will be devoted to probes of LNV via simple variants of the 
following process, which can be written schematically as 



Greek subscripts run over all different fermion flavors. Given the assumed democratic 
models, coupled with our present lack of experimental information, one would expect 
that all flavor combinations are equivalent to zeroth order. Any indication to the contrary 
would signify important deviations from simple expectations, and thus begin to reveal the 
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flavor structure of tlie new pliysics. Tlie above selected "golden modes" often yield the 
largest LNV rates, but this is not always the case. For example, some operators do not 
allow tree-level charged dilepton events, but rather prefer to include neutrino initial or 
final states. LNV processes with initial and final state neutrinos are extremely difficult 
to identify. The only hope of such discovery channels is, perhaps, via neutrino scattering 
experiments on either electron or nucleon targets, using well understood neutrino beams. 
I point out that any neutrino/anti-neutrino cross contamination induces ambiguity onto 
the total lepton number of the incident beam and would serve as a crippling source 
of background for LNV searches. This reasoning rules out conventional superbeam [35] 
facilities as well as proposed neutrino factories [36] , which contain both neutrino and anti- 
neutrino components, but does suggest modest possibilities for future beta-beams |37j . 
Given projected beta-beam luminosities and energies along with the derived cutoff scales 
Ky, it seems unlikely that LNV can be observed in such experiments. Another possible 
discovery mode involves only two external state quarks and an associated gauge boson as 
in the sample process laf^p —>■ du + W~ . It turns out that the rates for such processes are 
generally suppressed for the majority of operators involving six fermion fields, as I am 
trading a phase space suppression for a stronger loop suppression. For those operators 
with only four fermion fields, the situation is not as straightforward and, in some cases, the 
three particle final state is preferred. Typically, the neutrino mass induced cutoff scales 
of those operators are high (A^ ^ 100 TeV), so it would be quite difficult to observe 
such effects. Of course, any VT-boson final state will either promptly decay leptonically, 
yielding missing energy and unknown total lepton number, or hadronically, reducing the 
reaction back to that of the golden mode. 
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Another possibility is to replace two or more of the external quark states in Eq. (13.11) 
with leptons in such a way as to preserve charge, baryon number, and AL = 2 constraints. 
While many operators favor this structure, a little thought reveals that at least one 
external neutrino state is always present, which leaves only a missing energy signature, 
and little means of lepton number identification in a detector. Such events would not be 
clean, but of course, three final-state charged leptons and missing energy are enough to 
extract the existence of at least AL = 1 LNV, provided that the number of invisible states 
is known to be no greater than one. This last requirement is difficult to achieve in the 
presence of the large backgrounds and the limited statistics expected at future collider 
facilities, but should still be possible given a concrete model probed near resonance (see 
for example [38] ). Therefore, while important and potentially observable, this mode is not 
generally the best place to look for LNV and is neglected in the remainder of our analysis. 
From this perspective, the only other relevant channel of LNV discovery is related to W 
and Z rare decays into final states with non-zero total lepton number. This possibility is 
briefiy addressed in Sec. 13.2.21 

3.2.1. Neutrinoless Double-Beta Decay 

Here I probe the expectations for neutrinoless double-beta decay (/3/?0z/) for each operator 
listed in Table [3TT1 /?/50z/ is the LNV (AL = 2) process where, within a nucleus, two down 
quarks convert into two up quarks with the emission of two electrons but no neutrinos, 
or in the language of nuclear physics {A, Z) {A, Z + 2) + e~e~. See |39L 140] and 
references therein for a comprehensive review. While precise computations of nuclear ma- 
trix elements are essential for making detailed predictions |41] . the minimal parton-level 
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description given above is adequate for the purposes of this study. There is a continuing 
legacy of cutting edge experiments designed to search for PPOu with no success to datej^ 
Currently the ^^Ge half-life for this process is bounded to be greater than 1.9 x 10^^ yr and 
1.57 X 10^^ yr at 90% confidence by the Heidelberg-Moscow [43] and IGEX |44] exper- 
iments, respectively. Future experiments are poised to improve these limits (for several 
different nuclei) by a couple of orders of magnitude within the next five to ten years |45] . 

If one assumes that (3l3Qv proceeds via the exchange of light Majorana neutrinos, its 
amplitude is proportional to the ee element of the Majorana neutrino mass matrix. 



where are the neutrino masses and Uei are elements of the leptonic mixing matrix. With 
this, one can extract the upper bound rriee < 0.35 eV (90% confidence level bound, [340 
from current experiments while next-generation experiments are aiming at niee ^ 0.05 eVc 
|45] . In general, LNV new physics will lead to additional contributions to jSpOu, most of 
which are not proportional to rriee- However, the amplitude for (3(3Qv can still be expressed 
in terms of an effective mee, '"^ee , which is an operator-specific quantity that will be used 
to analyze new models of LNV. 

Here, I define six different "classes" of diagrams one can construct out of LNV ir- 
relevant operators that contribute to /?/50z/ at the parton level. These are illustrated in 

Fig. 13.31 and classified by the dimension of the generated LNV interaction, depicted by 

-"^There is currently a positive report of l3l3Qv at the 4.2ct level by a subset of the Heidelberg-Moscow 
collaboration [42] . They report a measured half-life of 1.74^° ;[g x 10^^ years which maps to rrigf ^ (0.2 — 
0.6) eV. I choose to neglect this controversial result, which is still awaiting independent conformation. 
^The parameter change from half-life to depends heavily on nuclear matrix element calculations. 
Current calculations induce an uncertainty of less then a factor of four on mee for most parent isotopes 
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large gray dots. In order to unambiguously separate the different classes, note that the 
grey circles are defined in such a way that all fermion and Higgs legs that come out of it 
are part of the "parent" operator O (and not attached on via reducible SM vertices), while 
all other interactions are SM vertices. The dots should be viewed as hiding the underly- 
ing LNV interactions. In general, they contain a mixture of coupling constants and loop 
factors that must be evaluated explicitly for each diagram. It is important to emphasize 
that the contribution of a generic operator O to l3j3Qv will consist of contributions from 
all different classes, while usually dominated by one of them. I show the lepton number 
conserving electroweak vertices (point-like) as effective four-fermion interactions, justified 
by the low energy scale of nuclear beta decays. The dotted lines indicate the exchange 
of VT-bosons, labeled by W and H (charged Higgs goldstone boson). Helicity arrows are 
explicitly included where uniquely determined, implying that the arrowless legs can have 
any helicity. 

Dy describes the standard scenario of (3(30i> mediated by light Majorana neutrinos. 
It is simply two electroweak vertices held together by a Majorana mass term on which 
two neutrinos are annihilated. The amplitude for this diagram is proportional to rriee, 
as defined in Eq. (13. 2p . The dependence on such a neutrino mass is intuitively clear 
considering the need for a helicity flip on the internal neutrino line. The remaining 
diagrams are qualitatively different from this standard case. Most importantly, none 
of them require "helicity flips" and are therefore not directly proportional to neutrino 
masses. They are, however, proportional to inverse powers of the new mass scale A^, 
and hence also suppressed. These effects are not entirely independent, since the value of 
Aj, was extracted from the requirement that neutrino masses are small, but correlations 
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Figure 3.3. The parton level Feynman diagrams contributing to neutrino- 
less double-beta decay, labeled by the dimension of the underlying lepton 
number violating interaction, indicated by gray dots. Each diagram is gen- 
erated, at some order in perturbation theory, by all analyzed interactions, 
but estimates of their magnitudes depend heavily on the details of the op- 
erators, including their associated scale Aj,, fermion content and helicity 
structure. 



are relaxed enough to allow nontrivial consequences. It is this partial decoupling from 
neutrino masses that allows larger than naively expected contributions to (3(30i> from some 
of the LNV irrelevant operators. Before proceeding, I make the trivial observation that 
the amplitudes following from D^, D4 and are additionally proportional to two powers 
of CKM matrix elements, namely iKidP, whereas Dq and -D7 are only proportional to one 
power of VudU The tree-level diagram Dq has no CKM "suppression." While this is a 
purely academic fact in the case of PPOu {\Vud\ ~ l'); it leads to important consequences 
for analogous rare decays that depend on the much smaller off-diagonal CKM matrix 
elements. I will return to these in the next subsection. 



"^This is true provided that I assume no flavor structure for the underlying operator, or, equivalently, 
that all dimensionless coupling constants are order one. If one is motivated by experiment to postulate 
a minimally flavor violating scenario, to perhaps ease constraints from flavor changing neutral currents, 
the statement must be modified accordingly. 
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For a given operator, the relative size of each diagram's contribution to the total 
decay rate depends on many factors including the operator's dimension, scale, fermion 
content and helicity structure. The dominant contributions must be calculated on a case 
by case basis. Generally, the high scale operators {K^ > 10 TeV) are dominated by the 
two dimension-four diagrams D^, and since many factors of Aj, will be canceled by 
divergent loops inside the gray dots thereby minimizing the suppression. All else 

being equal, is the strongest of the pair since it is enhanced by ~ Q^^ from the two 
propagating neutrino lines as opposed to only ~ Q^^ for the one neutrino case shown 
in diagram D4. For those operators with no tree-level vv field content, can still be 
very important, but its dominance is nevertheless rare. As discussed in Sec. 13.11 these 
are precisely the operators that have the greatest loop suppressions and consequently 
lower energy scales suggesting the need for diagrams beyond D4. The effects of low cutoff 
scale operators (A^, < 1 TeV) are not severely suppressed by l/A^, (by definition), so the 
dominant diagrams will typically be of the highest dimension allowed by the tree-level 
structure of the operator. For such low scales and for operators of the following schematic 
form dduuee (dimension 9) or dduueeHoHo (dimension 11), Dg always dominates the (3(3Qv 
rate yielding amplitudes proportional to 1/A^ and f^/A^, respectively. For intermediate 
scales, and when the operator's field content does not directly support /5/?0z/ due to 
lack of quark fields, the situation is not as straight forward and one must perform the 
relevant computations to determine the dominant diagrams. Still, it should be noted that 
diagrams containing internally propagating neutrinos are enhanced by inverse powers of 
Q and maintain a slight advantage over their neutrinoless counterparts. One can thus 
generally expect diagram Dq to dominate the decay rates for low scale operators when 
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Dg is suppressed. The opposite is true for interactions taking place at higher energies in, 
for example, next-generation colliders, as discussed in Sec. 13.2.31 

Since each diagram in Fig. 13. 31 can have different external helicity structures, the differ- 
ent contributions to the total rate will be added incoherently, thus eliminating the effects 
of interference. There are some case specific coherent contributions that I neglected in 
our treatment since most rates are dominated by a single diagram. Another potential 
difference among the different contributions is related to nuclear matrix element calcula- 
tions: can the calculations done assuming /3/30z/ via the standard light Majorana neutrino 
exchange scenario of diagram Di^ be applied to the more general cases encountered here? 
I have nothing to add to this discussion except to naively note that there is no obvious 
reason why such rates should be severely suppressed or enhanced relative to the stan- 
dard scenario. I therefore assume that all nuclear matrix elements are identical and can 
be factored out of the incoherent sum. I assume that this approximation is not more 
uncertain than the other sources of uncertainty inherent to our study (likely a very safe 
assumption). 

As drawn, each diagram Di contributes to the amplitude that characterizes (3(3Qv. For 
example, the amplitude associated with is proportional to 

\V 

Ad, = ruee ^ , (3.3) 

where Gp is the Fermi constant. The remaining diagrams will contribute with A^i oc 
({v,Q)A'^~^, up to a dimensionless coefficient containing various numerical/loop factors, 
as well as general scale depencies parameterized by some power of the ratio v/A. The 
function C{v,Q) has mass dimension i — 9 so that all have the same mass dimension. 
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Note that all aspects of are calculable given a LNV operator and diagram. I can 
analyze each operator in terms of an effective m^f , defined in terms of the underlying 
dimension nine amplitude by 



where i runs over the set {z/, 4, 5, 6, 7, 9} that labels the diagrams shown in Fig. 13.31 and 
Q ~ 50 MeV is the typical momentum transfer in P/SOu. rrf^ can be directly compared 
with experiment and used to make prediction for future observations. A few comments 
are in order regarding this quantity. First, it is a useful derived object that has no direct 
connection to a real neutrino mass and is valid to arbitrarily large values. Note that in 
the case of Majorana neutrino exchange, rrtgf = niee only if mee "C Q. When neutrino 
masses are greater than Q, rrf^ oc 1/m. Our definition of m^f also conforms to the use of 
large effective masses in |33] . The second comment is that, unlike the case of mee, which 
is valid for any process involving the exchange of electron-like Majorana neutrinos, mgf 
is case specific. It must be calculated separately for each process, as each one, in general, 
is composed of different diagrams. In particular, the calculations of the effective mass for 
/?/30z/ expressed here are not directly applicable to other LNV processes and should not 
be interpreted as such. 

The distribution extracted from all operators is shown in Fig. 13.41 assuming the 
scales derived in Sec. 13. H and color-coded for convenience within the histogram. Specif- 
ically, I indicate in green the operators that are characterized by sub-TeV scales and thus 
accessible to next-generation experiments via direct production. The blue and red op- 
erators are characterized by scales between (1 — 5) TeV and (5 — 25) TeV respectively. 




(3.4) 
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Figure 3.4. ml^ distribution derived for the neutrinoless double-beta decay 
process as described in the text. The calculations were made assuming 
the scales derived in Sec. 13.11 as well as universally order one coupling 
constants. The histogram bars are labeled explicitly with operator names 
and color-coded by their cutoff scales. Also shown in light gray is the region 
probed by next-generation experiments. The vertical axis is truncated at 
15 operators to best display the relevant features of the plot. 



where virtual effects should be most important for collider searches. The majority of op- 
erators, shown in cyan, are suppressed by scales greater than 25 TeV and are hence quite 
difficult to observe in other search modes. I also explicitly label each operator within the 
histogram bars for easy identification and comparison. One should notice the expected 
general trend that increasing leads to a decrease in ml^ and vice-versa. The vertical 
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axis is truncated at 15 operators, as the bar near 0.05 eV, dominated by the hght Ma- 
jorana neutrino exchange described above, would extend to nearly 100 operators. With 
broken vertical lines, I indicate the current 90% upper bound |34j, mf = 0.35 eV, and 
the potential reach of future experiments. 

This distribution, which spans over six orders of magnitude (from 10~^ eV till 10^ eV), 
reveals many important features of the effective operator set. Beginning at the largest 
mgf values, I find that the twelve operators appearing near 300 eV all have the expected 
common feature of low energy scales, including Oig with only just above the 1 TeV 
mark. Additionally, the contribution of the majority of these operators to (3j3Qv is dom- 
inated by the tree-level Dq diagram. The exceptions are (^54^ ^ and Oyo, all of which 
are characterized by sub 0.5 TeV scales and dominated by diagram Dq. Consequently, 
these are subject to a loop and Yukawa/gaugqj suppression relative to their dominated 
cousins, but the difference is not visible given the resolution of the figure. It is interesting 
to note that these three operators have the correct quark and lepton content for large 
jSpOu, but their SU{2)i gauge structures forbid large tree-level contributions. Similarly, 
operators Oie, (^55^ ^eSa? and C75 are also dominated by diagram Dq accompanied by 
slightly higher cutoff scales. This drives down nf^ significantly considering the leading 
one-loop scale suppressions of and A^^ for the dimension-eleven and dimension-nine 
operators respectively. I point out that operators 054^ ^ ^ ^ and O^^^ ^ ^ yield almost identi- 
cal expressions for their respective /3/30i/ amphtudes (as well as their radiatively generated 
neutrino mass expressions) with up and down quark Yukawa couplings exchanged. While 
this action enhances most of the (3f30i> couplings relative to those of 054, it also 



^As it turns out these are all suppressed by a single bottom quark Yukawa coupling as well as two powers 
of the SU{2)l gauge coupling 5, but this fact cannot be deduced from Fig. 13.41 alone. 
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raises the A;^ scale by nearly a factor of four and thus drives ml^ down by orders of 
magnitude. 

The remaining operators all predict rrtgf < 1 eV, close to current experimental bounds. 
The histogram bar near 0.1 eV is composed of operators of very different A^, scales. O34, 
and O5Q are both characterized by low cutoff energy scales around 0.5 TeV, but, due to 
their fermion and helicity structure, their contributions to jSpOu are dominated by two- 
loop versions of diagram Dq. The neutrino-mass-required cutoff for 073^ is around 7 TeV 
and its contribution to PPOu is also dominated by diagram Dq. In this case, however, 
the two-loop version turns out to be larger than the allowed one-loop amplitude due to 
strong scale suppressions (the added loop reduces the cutoff dependency from A^^ to A"'^). 
This behavior is characteristic of operators with a larger value of A^. The A^, = 40 TeV 
operator O20, defines the lower edge of this histogram bar. It is dominated by the one- 
loop diagram Dq enhanced by a top quark Yukawa coupling and, being a dimension nine 
operator, is only suppressed by A~'^ from the start. The next bar down contains operators 
dominated by D,,. Most of these are suppressed by a very high energy scale, but a small 
subset is characterized by scales A,y < 25 TeV. In particular operators On, Ois and 
are all cutoff at 2 TeV but, due to their fermion content they cannot participate in 
any of the non-standard interactions of Fig. 13.31 at a low enough order in perturbation 
theory. Similarly, the intermediately scaled operators O^^^^^, O50, Oq2, C'es? and 
Oeg^ have either the wrong fermion content or gauge structure to enhance any of the /3/30z/ 
diagrams (other than Di,) to an observable level. These operators are important because 
their minimal forms are experimentally unconstrained yet still potentially observable to 
both next-generation f3f30i' and collider experiments. The remaining histogram bars with 
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are not accessible to PPOu experiments in the foreseeable future. Each of 
these diagrams are dominated by D^, either due to high suppression scales as in the case 
of O7 and Os, or, as in O35 and O58, the operator's fermion content simply disfavors other 
contributions to the f3(30i> amplitude. It is the general form of the neutrino mass matrix 
derived in Table ISTTl where I see that rriee oc i/e, that drives these operators away from their 
peers near ml^ = 0.05 eV. It is unfortunate that the two "low" dimensionality operators 
Oj and Os are cutoff by energy scales A^, in excess of 100 TeV and are hence invisible 
to any direct probe. If either of these operators have anything to do with nature, it is 
unlikely that LNV will be observed in the foreseeable future in any experiment. On the 
other hand, any observation of LNV will rule out these types of scenarios. Additionally, as 
will become clear shortly in Sec. 13. 3[ current neutrino oscillation data already marginally 
disfavor such operators and have ample room to tighten constraints in the near future. 

It is interesting to point out that the lower boundary of the currently excluded region 
falls within the m^f distribution, suggesting exciting prospects for the future. That being 
said, one should not read too much into current and future null results as, for most 
operators, relatively small cancelations and order one factors, not accounted for here, can 
push the relevant rates below the observable level depending on the underlying ultraviolet 
theory. On the other hand, one is allowed to interpret that operators that lead to m^g > 
10 eV are severely constrained (if not ruled out) as proper explanations for neutrino masses 
if one assumes the new physics to be flavor "indifferent" - order one factors cannot be 
evoked to save the scenario. Once this assumption is dropped, however, it is quite easy 
to "fix" these scenarios, since the large f]f30i> rate is a direct consequence of the universal 
order one couplings and the relatively low cutoff energy scale A^. For example, one can 
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suppress the coupling of new physics to first generation fermions (compared to second 
and third generation fermions) , thereby suppressing the worrisome diagrams of Fig. 13.31 
This will have little effect on the relation between A and the neutrino masses, discussed in 
Sec. 13.11 since these are either generation independent or highly reliant on third generation 
Yukawa couplings. Of course, by combining jSpOu searches with other probes I can obtain 
a much better idea of the origins of LNV as well as the relevant model(s), if any, chosen 
by nature. 

3.2.2. Other Rare LNV Decay Processes 

Most of the qualitative discussions of Sec. 13.2.11 devoted to PPOu, can be directly applied 
to other rare decay processes with the same underlying kernel interaction described by 
Eq. ( 13.11) . For such processes one need only analyze simple variants of the diagrams listed 
in Fig. 13.31 using crossing amplitude symmetries to account for the needed initial and 
final state fermions. Other factors must be added to the various electroweak vertices to 
account for quark flavor mixing. The requisite CKM matrix elements can highly suppress 
many diagrams for processes involving cross-generational quark couplings. In fact, only 
tree- level Dg diagrams are safe from such suppressions. Next, and most importantly, one 
must include the appropriate characteristic momentum transfer Q of the new system. 
Specific rates are highly dependent on this quantity as effective operator cross-sections 
typically grow with some power of Q. The particular exponent of the power law depends 
on the diagram, but naive dimensional analysis dictates that F oc Q^"^ for diagram Dg, 
rendering it highly dependent on a reaction's energy transfer. The fact that each diagram 
varies with Q in a different way implies that predicting the dominant contributions to a 
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given process is non-trivial and must be addressed quantitatively. Finally, in the cases of 
hadronic decays, one must also account for initial/final state matrix elements. I assume 
that all factors can be simply estimated on dimensional grounds. 

Unlike the j3j3Qv case, some meson decay modes proceeding via new LNV tensor in- 
teractions are expected to be suppressed. Such processes are one instance in our analysis 
where an operator's Lorentz structure can qualitatively affect expected LNV decay rates. 
One can understand this by considering a meson decay mediated by a new tensor parti- 
cle. The parton level interaction has the form {ua^yd)T^^'^ where the initial state quarks 
are explicitly shown and all other fields are contained in the tensor T^'^. Following the 
standard procedure I factor out the hadronic structure in the form of a free decay con- 
stant and write the amplitude as generally allowed by Lorentz invariance in terms of the 
external state's four-momentum. Due to the antisymmetry of a^v, this amplitude van- 
ishes to first order. Non-zero contributions to this decay mode must necessarily involve 
individual parton momenta and are therefore suppressed relative to the usual vector-like 
decay calculations. From this, it is clear that models of LNV containing tensor couplings 
will often evade the predictions and bounds of this section. Tensor operators will mediate 
LNV meson decays into more complicated final states (one may include, say, initial/final 
state radiation). Associated rates are, however, subject to additional gauge coupling and 
phase space suppression that tend to further reduce the already tiny LNV rates beyond 
any hope of detection. 

Rare LNV meson decays have been experimentally pursued for many years [34] . Here, 
I focus on the AL = 2 processes M' M+i^i^, where M' and M are the initial and final 
states mesons respectively and the is represent like-sign lepton pairs of arbitrary fiavor. 
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H,W 

Figure 3.5. The parton level Feynman diagrams contributing to rare LNV 
meson decay labeled by the dimension of the underlying lepton number 
violating interaction, indicated by gray dots. Each diagram is generated, at 
some order in perturbation theory, by all analyzed interaction, but estimates 
of their magnitudes depend heavily on the details of the operators, including 
their associated scale A^, fermion content, and helicity structure. 

Electric charge conservation dictates that M' and M have equal and opposite charge. 
Here I take each meson to consist of a color singlet up-type/antidown-type bound state 
and factor out all long distance hadronic effects. In this way 1 can view the meson decay 
process as du ia'^f^ + du for all up-type and down-type quark flavor combinations. The 
effective LNV diagrams contributing to this process are shown in Fig. 13.51 with the same 
naming scheme as their analogs in Fig. 13.31 Here, V and V denote potentially distinct 
elements of the quark mixing matrix. I additionally point out the potential dependency on 
all entries of the Majorana neutrino mass matrix elements m^^ in diagram D^, as opposed 
to the PPOu case where Di, depends only on rriee- These processes probe combinations of 
the neutrino masses that are naively unconstrained by (3j3Qu |46j . In general, the varied 
flavor structures encountered in meson decays allow for experimental probes into new 



"'^For simplicity I assume that both the process M' — > M and its conjugate have similar amplitudes 

and therefore treat them symmetrically. Large CP-violating effects can invalidate this assumption. 
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physics couplings across the fermion generations. I pointed out earher that some of the 
LNV operators lead to unacceptably large rates for PPOu unless first generation quarks 
participate in the new interactions with severely suppressed couplings (compared with 
second and third generation quarks). If such a scenario is realized in nature, rare D or B 
decays may be much more frequent than naive expectations. For this reason, improving 
rare decay sensitivities to all channels is essential to completely constrain models of new 
LNV physics beyond the minimal framework analyzed here. 

Reference |33] summarizes LNV upper bounds on all of these processes in terms of 
the effective Majorana neutrino mass matrix element that one would extract from 
observation assuming that all decay rates are dominated by the light neutrino exchange 
shown in D^. Hence, I can compare operator expectations with current experimental 
limits in exactly the same way as was done in Sec. 13.2. 1[ For a given LNV meson decay, 
is defined from the contribution of the different classes of diagrams to the rare 
meson decay in question, exactly as was defined in the previous subsection (see 
Eqs. (13. 313.41) ). Direct estimates for different process reveal m^^ distributions similar to 
that for m^f depicted Fig. 13.41 up to "rescalings" that reflect the different kinematics 
and the presence of small CKM mixing matrix elements. Results are summarized in 
Fig. 13.61 for a representative sample of charged meson decays. Each histogram is labeled 
by its associated decay mode and is color-coded to indicate the neutrino-mass constrained 
cutoff scale A^, of the different LNV effective operators. For simplicity, I refrain from 
listing operator names on the individual histogram bars (as opposed to what was done 
in Fig. 13.41) . The "operator ordering" is very similar to that of Fig. 13.41 especially in 
the low Aj, scale, high effective mass regime where decay rate predictions are particularly 
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Figure 3.6. m^^ distribution for several rare LNV meson and r decays. 
Calculations assumed the charge lepton flavors ia^^ = /^e, while the r decay 
histogram (lower right-hand panel) was obtained assuming the final state 
mesons MM' = KK. The histogram bars are color-coded by suppression 
scale. Current bounds on these processes are typically above 1 TeV and are 
not visible at these small scales. 
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important. Note that the horizontal axes are relatively fixed for easy comparison and 
that the vertical direction is truncated and does not reflect the true "height" of the lowest 
mass bar (order one hundred operators). 

Specifically, I present effective Majorana neutrino mass distributions for the processes, 



reading down the panels from left to right, D — > + D K + i'^i^, Dg ^ -k + 



K + iti% B 



ir + iti% B 



7r + 



V±f± 



as well as the rare r decay. 
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MM' + Here the final state leptons can be of any flavor allowed by energy 
conservation. Since, as previously discussed and explicitly verified numerically, the specific 
details of the distributions are mainly dictated by kinematics and CKM matrix elements, 
these results are robust under changes in the final state lepton flavors. The r decay 
distribution shown in the lower right panel is representative of all possible decay products 
including first and second generation charged leptons and light meson states. One should 
notice the expected general operator trend within each histogram as the characteristic 
cutoff scale is decreased, as well as the expected peaks near 0.05 eV dominated by light 
Majorana neutrino exchange. Additionally, each distribution is much "broader" than 
the one in Fig. 13.41 This observation exemplifies the fact that effective mass calculations 
depend critically on the underlying process. Indeed, maximum m^f values can reach nearly 
10^° eV for the —>■ K~ +e~^e~^ decay but only 10^ eV for jSpOu. Current upper bounds 
for m^^ from these processes, mostly well above one TeV, are well beyond the largest 
operator predictions here, ranging from < 0.09 TeV for the case of 7r~e'^/i^ to 

m^^ < 1800 TeV for the case of K^fi^fx'^ |33j . It is curious that the best meson 

decay bounds come from rare LNV kaon process but, as can be seen in the lower left 



panel of Fig. 13.61 these yield by far the lowest predictions. Future experiments have the 
potential for observing LNV for a select few operators only provided vast improvements 
in meson production luminosities. Current and upgraded B-factories [4Tj are expected 
to provide the most significant improvements, considering the large derived 5-meson 
effective masses shown in Fig. 13.61 Still, the best cases from the figure yield only the tiny 



^The actual calculations displayed in Fig. 13.61 assumed the charge lepton flavors £a(f3 = /^e, while the r 
decay histogram (lower right-hand panel) was produced assuming the final state mesons MM' — KK. 
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branching fraction 1.8 x 10^ for the case of the rare decay iv^e'^jj,'^, nearly eleven 

orders of magnitude below the current experimental limit of 1.3 x 10^^ |34] . 

Another possible search mode involves the decay of the Z-boson into LNV final states. 
The dominant contributions to this process are generally unrelated to the reactions sum- 
marized in Eq. (13. ip and shown schematically in Figs. 13.31 and 13.51 While there is a slight 
connection between them as one can always attach a Z-boson to various fermion lines in 
each diagram, there are potentially large lower order contributions arising within the op- 
erators themselves. The latter, when present, can easily overtake the associated "golden 
mode" counterparts. In this context, such processes can be thought of as the decay of 
the longitudinally polarized Z-boson. Strict bounds exist on such decays from the LEP-I 
j48j and SLD [49j experiments. Each element of the operator set predicts decays into 
final state fermions with total lepton number L = 2. The dilepton pair can be of any 
fiavor and is generally accompanied by two or four additional fermion states, depending 
on the dimension of the operator. I restrict our discussion to the dimension-eleven oper- 
ators comprising the majority of the sample, as these are typically suppressed by lower 
cutoff Ay scales and, equally important, explicitly contain Higgs doublets in their field 
content. In this case, tree- level decays result in a six-fermion final state which suffers from 
a large phase space suppression and cumbersome multiplicities that are likely to render 
even the most sophisticated search ineffective. The only possibility of this type that yields 
a charged dilepton signal is Z — > ^t^t^QQQ (quarks of all allowed flavors implied), but 
many other possibilities exist involving invisible flnal state neutrinos. A little thought 
also reveals that closing fermion loops in an attempt to obtain simpler flnal states and 
thus render the analysis more tractable will necessarily result in flnal state neutrinos. 
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Therefore, the majority of the Z-boson LNV decay channels involve invisible final states 
with practically undetermined total lepton number. The prospect of direct discovery by 
these means seems dismal, but indirect constraints on LNV are still possible from bounds 
on the Z-boson invisible decay width. There is currently a statistically insignificant, but 
nonetheless captivating, 2a deviation between the observed invisible decay width and SM 
expectations assuming three light neutrino species [34j . The experimentally extracted 
branching ratio was found to be slightly smaller than its predicted value so that a new 
LNV contribution of the form Z would push the invisible branching ratio in 

the "wrong" direction. From these bounds the decay width of any new contribution to 
the Z-boson decay is constrained to be less than 2.0 MeV at the 95% confidence level 
[341 150] . A quick estimate reveals that this constrains the dominant LNV amplitudes 
Az < A/47r(2.0 MeV/Mz) ~ 0.53. For the dimension-eleven operators of interest, the 
largest possible amplitude is of order y^/(167r^)^(f /A)'^ where y is an arbitrary fermion 
Yukawa coupling and four powers of the cutoff scale A are removed by divergences in the 
closed diagram loops. The constraint above translates into y'^{v/A)^ < 4.1 x 10^, which 
is easily evaded by even the best case scenario of y = ?/t ~ 1 and A ^ v. Experimental 
bounds on Tin^ must be improved by a factor of a million before they start significantly 
constraining LNV (under the assumptions made here). This result holds for virtually all 
possible flavor structures. I conclude that rare Z-boson decays are not practical discovery 
modes for the LNV effects considered here, but look to future rare Z-boson decay studies 
for more information. 

In a similar way, one can also dismiss the case of rare IV-boson decays as promising 
probes of LNV. As in the Z-boson case, the H^-boson can decay into a variety of L = 2 
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final states proceeding either through couphngs to left-handed fermion lines or explicit 
operator content. Here, however, there is no six-fermion, same-sign dilepton final state 
with no neutrinos due to conservation of charge and weak isospin, so the lowest order 
observable mode is already loop suppressed to W~ — >• ^^^z? + IQ- Current TV-boson decay 
bounds are far too weak to constrain such suppressed LNV (34j and are not likely to 
improve to the level implied by the operators under consideration, which predict the tiny 
decay rate Tl^v < mvi/(47r)/(167r^)^(t>/A)"'^° ^ 10"^ MeV in the best case scenario of 
electroweak scale A^,. I also point out that, contrary to the Z-boson decay limits, there 
are no robust, indirect bounds that can be used to constrain LNV in the case of the W- 
boson. Note that, despite dismal prospects for gauge boson decay driven LNV discovery 
within the minimal framework of "natural" effective operators, one can still construct 
theoretically well-motivated models that will yield observable signals. Particularly, in 
a weak-scale seesaw mechanism {Oi), the new degrees of freedom, comprised mostly of 
Majorana gauge singlet fermions (right-handed neutrinos), can mediate visible, AL = 2, 
VF-boson mediated processes with little or no scale/loop suppression. This class of model 
is analyzed in |51j and is exempt from the discussion outlined here. 

3.2.3. Collider LNV Signatures 

If neutrino masses are a consequence of ultraviolet physics related to cutoff scales around 
the TeV scale, I expect future high energy collider searches to directly access the new 
LNV physics. For example, the direct, resonant, production of new states could lead to 
rather spectacular signals of these models. It would also indicate the breakdown of the 
effective field theory approach undertaken here. To pursue such possibilities, one must 
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assume a specific ultraviolet sector and study its signatures and implications on a case by 
case basis. In the looming shadow of the LHC, [52] and the more distant ILC [53] . such 
an analysis is highly warranted but will not be pursued here. Instead, I assume that the 
masses of new ultraviolet degrees of freedom remain out of the reach of next-generation 
accelerator experiments. Such a situation can be easily accommodated within the context 
of the preceding results, considering the order of magnitude nature of the estimates. 

I will concentrate on the process e^e^ QQQQ (which will usually manifest themselves 
as jets) with no missing energy in an ILC-like environment |53] with a center-of-mass 
energy of 1 TeV and an integrated luminosity of 100 fb^^. I also make the oversimplifying 
assumption that the detector system has equal acceptance to all quark flavors, and the 
ability to efficiently distinguish quarks, gluons and rs. By summing over all possible quark 
final states it is simple to estimate the total LNV cross section for each effective operator, 
assuming it is responsible for neutrino masses. Such searches can be complemented by 
looking at e~e~ — > W~W~, which have been discussed in detail in the literature [54]. As 
discussed in Sec. 13.2.21 the different LNV operators couple to one or more gauge bosons 
via an appropriately closed fermion loop or direct coupling to the Higgs doublet field. 

Charge and baryon number conservation dictate that the two quarks in e~e~ QQQQ 
are down-type quarks, while the two antiquarks are up-type antiquarks. At the parton 
level, the scattering process is similar to /5/50z/, which motivates exploiting simple vari- 
ations of the diagrams in Fig. 13.31 in order to calculate the relevant amplitudes, as was 
done in Sec. I3.2.2[ Here, the extensions are obvious: use crossing symmetry to rotate all 
lepton lines into the initial state and all quark lines to the final state taking special care 
to insert appropriate CKM matrix elements where needed. Due to the large characteristic 
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momentum transfer Q of the e~e^ scattering, one must also "expand" the electroweak 
vertices and account for gauge boson propagation. With this in mind, the amphtude 
calculations can be carried over directly from the previous sections. Specific results are, 
however, quite distinct due to the higher center-of-mass energies involved. In the lan- 
guage of the underlying diagrams mediating this reaction, for diagrams characterized by 
TeV cutoff scales, diagram Dg, if allowed at tree- level, will dominate the rates. As in the 
previous cases, for intermediate to high cutoff scales, general diagram dominance must 
be addressed on a case by case basis. It is important to appreciate that, since these are 
non-renormalizable effective interactions, cross-sections grow with center-of-mass energy. 
For this reason, I expect many of the low cutoff scale operators to yield observably large 
signals at the ILC. 

Fig. 13.71 shows the e~e~ qqqq cross-section distribution, in femtobarns, at the ILC, 
calculated for all 129 of the analyzed LNV operators. Once again, the extracted value 
of the cutoff energy scale assuming constraints from neutrino masses are color-coded 
to indicate operators associated with a low (Aj, < 10 TeV) or high (A,^ > 10 TeV) 
ultraviolet cutoff. Each bar is also labeled with the respective constituent operators, for 
convenience. Note that the vertical axis is truncated at fifteen operators (the left-most 
bin is over 60 operators high) to help clearly display relevant features of the plot. I 
also highlight the potential reach (defined as cross-section greater than the inverse of 
the integrated luminosity) of the ILC with a broken vertical line, assuming 100 fb"^ of 
integrated luminosity. This particular ILC luminosity value should be considered as a 
loose lower bound, introduced to give a feeling for the observable scales involved. It has 
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Figure 3.7. Distribution of total cross-section for the process e~e~ — > qqqq 
and no missing energy at an e~e~ collider with 1 TeV of center-of-mass 
energy. Estimates were obtained assuming the scales A^, derived in Sec. 13.11 
as well as order one coupling constants. The histogram bars are labeled 
with operator names and color-coded by A^, cutoff scale. Also shown (bro- 
ken vertical line) is the reach of such an experiment assuming 100 fb"^ of 
integrated luminosity. The vertical axis is truncated to best display the 
relevant features of the plot. 



recently been argued, for example, that a realistic machine should be able to outperform 
this estimate by over an order of magnitude |53j . 

A glance at Fig. 13.71 reveals that it generally adheres to the expected correlation of 
decreasing A^ scales with increasing LNV rates, similar to what is observed for other 
LNV observable {e.g., Fig. 13.41) . The similarities between the different processes extend 



beyond mere trends to the specific ordering of the operators within each histogram. This 
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reflects the common underlying interactions that drive these processes. The operators 
on the far right of the plot, topping off the highest cross-sections, are exactly those 
operators with the largest m^f, now "split" into three different bars. The large bar 
just below 10^ fb is composed of sub-TeV scale operators with tree-level diagram Dg- 
like fermion content. Slightly smaller are the expectations for Oig, again dominated by 
diagram Dg, but characterized by a slightly larger scale (around one TeV). Moving 
down in cross-section, this is followed by the low cutoff scale operators 054^^ and Ojo, 
dominated by a combination of diagrams Dq and D7. On the opposite end of the plot 
I point out the large bar below 10^^^ fb, composed mainly of operators associated to 
high cutoff scales (A^, > 25 TeV). The contributions of these operators are dominated by 
light Majorana neutrino exchange, but their histogram bar contain far fewer models than 
their (3(3Qv counterpart, as many of the latter have been driven up due to new diagram 
^4 and contributions. In general, the large center-of-mass energies tend to magnify 
differences between interaction rates that were not relevant in low-energy observables. 
This naively suggests that high energy probes have a higher potential for distinguishing 
different models. 

There are eleven operators that lead to an observably large (as defined earlier) e~e~ 
qqqq cross-section at the ILC. Note that all of these were already "ruled out" by current 
/?/30z/ searches. As discussed in Sec. 13.2. H however, these bounds only effectively limit 
the couplings of the new physics to first generation of quarks and leptons, and hence, 
if such a scenario is realized in nature, one should still expect large contributions from 
decay modes that lead to second and third generation final state quarks. In fact, even 
one such heavy quark is enough to bypass the constraints from PPOu for several effective 
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operators. Such reasoning implies that constraints on the new physics flavor structure can 
be made quite strong at a linear collider via analyzes of the flavor of the final state quarks. 
By identifying and comparing the outgoing quark flavor one can extract individual limits 
on quark-lepton coupling constants within the operators. Additionally, kinematics can 
be used as a further operator probe. For example, one can potentially determine the 
dominant underlying LNV diagram (say Dg, D-j or Dg) by checking whether the various 
kinematic distributions are characteristic of VT-boson exchange. 

The ILC can cleanly select or discard some LNV scenarios. This characteristic is 
further enhanced by considerations of initial electron polarization. Planned linear colliders 
have the ability to produce partially polarized beams (80% polarization for e~, 40% for 
e"*" |53L 155] ). The power of a high energy polarized e~e~ beam is in model identification 
and rejection. Of all operators that yield observably large cross-sections, the e^e^ mode 
can only probe 053, and therefore any positive LNV signal cleanly identifies this as the 
operator chosen by nature. In a similar way, the ILC running in its e^^e^ mode can easily 
observe LNV from O19, 054^, £^54^, O59 and Oqq] and to a lesser extent, operators 054^ 
and C70, and possibly even 054^. Finally the e^e^ mode can probe operators C36, O37 
and Osg. Within this framework, any LNV detected in one ILC polarization mode will 
generally not be seen in the others. This statement also applies to resonantly enhanced 
low scale operators that lie outside the observability window. 

While e~e~ collisions only probe effective operators that "talk" to first generation 
leptons, there are several lepton collider processes that allow one to explore other members 
of the charged lepton family. Future high energy muon colliders |56] could, in principle, 
also be used to study LNV. In this case, all of the preceding discussions regarding the 
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ILC are applicable. Electron linear collider facilities can also be used to study 76^ and 
77 collisions [57] . 76" collisions can be used to probe 76" ^ £^ + X (and hence the 
"ea" structure of different LNV operators), while 77 ^t^t + ^ probes all the different 
a, (5 charged lepton flavors. For 77 collisions, for example, considering projected ILC-like 
collider parameters, one would expect the same operator distribution as Fig. 13.71 shifted 
down in cross-section by, roughly, a factor of ~ 10^^. Thus, a handful of operators 
should be testable at a future 77 collider assuming 100 fb~^ of integrated luminosity. 

The preceding analyses carry over to the case of hadron colliders, such as the LHC, in 
a relatively straightforward way. The LHC, or Large Hadron Collider, is a proton-proton 
machine that will operate at a center-of-mass energy of 14 TeV and a characteristic 
integrated luminosity around 100 fb~^ |52j (in its high luminosity mode). The relevant 
LNV variants of Eq. (13.11) are dd — > i~i^uu and uu — > i'^i^dd with no missing energy. 
Of course, at center-of-mass energies well above a TeV, the proton-proton collisions are 
dominated by the gluon content of the proton, so most interactions at the LHC will be 
initiated by gluon-gluon and gluon-quark scattering. The dominant LNV subprocesses 
are qg ^t^^^QQfl and gg ^t^^^QQQQ and are illustrated in diagrams (a) and (6) of 
Fig. 13.81 respectively. These are characterized by similar final states as the quark-quark 
scattering reactions but, given that there is no explicit gauge boson field content in the 
LNV operators in question (Table [371!]) . their amplitudes are proportional to unimportant 
order as and coefficients, respectively. The parton level diagram (c) shows the related 
process gg ia^p + 99- The rate for this process can be estimated, relative to its 
four jet cousins, by exchanging a final state phase space suppression for a single loop 
suppression. In all three diagrams depicted in Fig. 13. 8[ the LNV interaction regions 



105 



represented by large grey dots contain all of the diagrams discussed earlier, meaning that 
the operator amplitudes calculated for the ILC can be recycled in this analysis. While all 
three bare diagrams are characterized by rates of the same order of magnitude diagram 
(c) leads to missing transverse energy and potentially undetermined final-state lepton 
number, rendering it a less than optimal experimental search mode. Note that, in all of 
these cases, the external, and internal, fermions outside of the LNV interaction region 
can be of any flavor. Therefore, hadron collider experiments have, in principle, access to 
all LNV operator parameters. Cleanly identifying and constraining all said parameters 
should prove quite difficult for all but the most obvious signatures. The above statements 
regarding signals at the LHC are also applicable at the Tevatron with some minor, but 
important, modifications. The Tevatron's pp collisions are at a much lower center-of-mass 
energy, roughly 2 TeV, while the total expected integrated luminosity, less than 10 fb~^ per 
experiment, is orders of magnitude smaller. These factors lead to much lower amplitudes, 
reduced by approximately a factor of (QTcvatron/QLJ^c)^ ~ 10~^|j The smaller center- 
of-mass energy also limits the Tevatron's ability to directly produce new physics states. 
With this in mind I conclude that the Tevatron has little or no chance of discovering LNV 
(within this minimal framework). 

A detailed set of predictions for the LHC would require a much more refined analysis, 
including the effects of parton structure functions, fiux distributions, and backgrounds, 
and as such is beyond the scope of this general survey. I would, however, like to point out 
that some of the reactions outlined here are subject to large background rates. While SM 



"^Strictly speaking one must also account for the proton's structure functions at the Tevatron's energy 
scale. Unlike the LHC, where collisions are dominated by gluon-gluon interactions, proton collisions 
at the Tevatron are dominated by valence quark interactions. These considerations do not affect our 
conclusions. 
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Figure 3.8. Parton level gluon-gluon and gluon-quark LNV interactions 
relevant at high energy hadron colliders. Each of these yields a same sign 
dilepton signal with jets and no missing energy. Notice that the final state 
fiavor structure is completely arbitrary under the assumption of random 
order one coupling constants. 



processes are lepton number conserving, many can fake the LNV signals in the complicated 
environment of a high energy hadronic interaction. The requirement of no missing final 
state energy is particularly hard to accommodate as some energy is always lost down 
the beampipe. As is typically done, one must rely on the less restrictive conservation of 
transverse momentum in order to constrain invisible states, such as neutrinos. SM same- 
sign dilepton production processes arising from, say, ly-boson pair production, are serious 
potential sources of background. Furthermore, it is impossible to predict correlations 
among final state jets without selecting a particular operator and underlying model of new 
physics, making it difficult to impose general cuts to reduce other hadronic backgrounds. 
Of course, some of the low scale LNV operators yield large enough total cross-sections that 
even crude analyses may suffice to reveal their existence. I conclude by pointing out that 
a large amount of recent work has been dedicated to LNV searches at collider facilities 
[58] . Most of these approach the subject from the perspective of sub-TeV mass, mostly 
sterile Majorana neutrinos that mix with the active neutrinos and are thus related to 
light neutrino masses via the seesaw mechanism [7]. This amounts to one example that 
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leads to the dimension- five operator Oi, but where one assumes that the propagating 
degrees of freedom are twelve or thirteen orders of magnitude lighter than the ultraviolet 
cutoff scale A^^jj In this case, LNV interactions are dominated by diagram Di, of Fig. 13.31 
(where heavy (weak scale) neutrinos are also exchanged), and as such one should make 
use of specific kinematic cuts to reduce background rates. These cuts, however, may also 
remove LNV signals resulting from many of the scenarios explored here, particularly those 
whose rates are dominated by Dg at tree-level. I urge experimentalists to account for this 
possibility while analyzing future data sets. 

3.3. Neutrino mixing 

Table 13.11 contains predictions for all the entries rriais, the Majorana neutrino mass 
matrix. These are computed in the weak basis where the weak interactions and the 
charged-lepton Yukawa couplings are diagonal, so that the eigenvalues of the neutrino 
mass matrix are the neutrino masses (bounded by oscillation experiments and, say, pre- 
cision measurements of tritium beta-decay [23]), while its eigenvectors determine the 
neutrino mixing matrix, constrained mostly by oscillation experiments. Since different 
LNV effective operators predict different flavor-structures for the neutrino mass matrix, 
there is the possibility to constrain the different scenarios with existing oscillation data 
[16] . While I can only predict the values of m^^ within, at best, an order of magnitude, 
it is still possible to extract useful information from the derived large scale structure of 
the expressions. In particular I can test the hypothesis of whether A values associated 
to different lepton flavors are allowed to be of the same order of magnitude. In order to 



^This can be achieved in two different ways. Either the new physics is very weakly coupled, or the new 
physics ~ SM couplings are finely-tuned [29] . In order to observe right-handed neutrinos in colliders, the 
latter must be realized. 
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obtain more accurate predictions and further probe the fine details of lepton mixing one 
must succumb to specific models, beyond the scope and philosophy of this analysis. 

The mass matrix for the three light Majorana neutrinos can be reconstructed from 
nine observables: three masses mi,m2,m3, taken to be real and positive; three (real) 
mixing angles 6^12, 6^23, 6*13; and three CP- violating phases 5,02,03. Here, 5 is a so-called 
Dirac phase that is generally present in the system regardless of the neutrino's nature 
(Majorana or Dirac fermion), while 0i,02 are so-called Majorana phases, only present if 
the neutrinos are Majorana particles (which is the case of all scenarios under consideration 
here). Oscillation data determine with relatively good precision 6*12, ^23, Amfg = m\ — m\ 
and lAm^gl = \m\—mW. I define neutrino masses such that mi < 1712 and Am^2 < l^^^isl; 
so that the sign of Am^g remains as an observables which characterizes the neutrino mass 
hierarchy ("normal" for Am^g > 0, "inverted" for Am^g < 0). See, for example, [14] 
for details. As for the third mixing angle, sin^ 6'i3 is constrained to be less than 0.025 
(0.058) at 2cr (4(t) from a three neutrino global oscillation analysis [16]. A considerable 
amount of uncertainty remains. In particular I have only upper bounds on the absolute 
neutrino mass scale, from kinematical measurements such as tritium beta decay [171 118j . 
plus cosmological observations [19L I20L I21L 122] . Finally, the three CP violating phases 
are completely unconstrained, and I have no information regarding the neutrino mass 
hierarchy. 

The above experimental results allow for several different "textures" for iriap in our 
weak basis of choice (see, for example, [59j). The purpose of this section is to discuss 
whether any of the textures predicted by the different LNV effective operators is "ruled 
out" by current observations. Most of the analyzed operators imply "anarchic" [24] 
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neutrino masses. This simply means that all elements of the neutrino mass matrix are 
uncorrelated and of the same order of magnitude. This hypothesis is known to "fit" the 
current data very well [24j . It will be further challenged by searches for di^ (the anarchic 
hypothesis favors large ^13 values) and probes that may reveal if the neutrino masses are 
hierarchical or whether two or three of the masses are almost degenerate (anarchy naively 
predicts the former). If future data strongly points towards non-anarchic m^,^, I will 
be forced to conclude that there is nontrivial "leptonic" structure in the dimensionless 
coefficients A of most of the LNV operators considered here. 

Many of the operators associated with a low neutrino-mass related cutoff scale (A^, < 
10 TeV), on the other hand, naively predict more structured neutrino mass matrices. 
Operators 

Ot, Os, Ci9, O20, Ou, 054a,6,c,d) ^55aM.cy ^56, C'sT, C^SS, C'sQ, Oqq, O-jq, (^75, (3.1) 

which radiatively generate neutrino mass elements proportional to distinct charged lepton 
Yukawa coupling {ije^y^^yT), yield mass matrices m such that 



Ve y^i Vr 



\ 



m oc 



Vf^ y^i Vr 



(3.2) 



\ Vr Vr Vr J 
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Additionally, models described at low energies by O36, and O^s generate neutrino 
masses proportional to both associated charged Yukawa couplings, such that 



^ VeVe yeVfi VeVr ^ 



m oc 



(3.3) 



VeV^i y^,y^l Vf^Vr 

\ VeVr y^lyT VrVr J 

The strongly hierarchial nature of the charged lepton masses {ye ^ ^ Vt)-, implies 
that the m^^ elements of Eqs. (13. 2p and (13.31) are expected to be hierarchical as well. 
In particular, the ee matrix element, mee, proportional to Ue or y1 is, for all practical 
purposes, negligibly smal|3 in both of these cases. On the other hand, it is well known 
that only a normal neutrino mass hierarchy is consistent with vanishing rriee |60] . so that 
both Eqs. (13. 2p and (13.31) predict the neutrino mass ordering to be normal. In the absence 
of extra structure, scenarios characterized by the LNV operators listed in Eq. (13. ip plus 
C36, C's? and O38 will be ruled out if future data favor an inverted mass hierarchy, or if 
the neutrino masses end up quasi-degenerate (regardless of the hierarchy). As will become 
clear shortly, Eqs. (13.20 and (13.30 predict that the lightest neutrino mass (mi in this case) 
is small (< A/Amf^). 

A more detailed analysis reveals that naive expectations from Eqs. (13.20 are already 
disfavored, while those from Eqs. (13.30 are virtually excluded. Assuming the normal 
hierarchy and very small mee, one can find a relation between the neutrino mass eigenstates 
and the oscillation parameters, thus reducing the number of free parameters in the mass 
matrix by one. Consider the diagonalization of the neutrino mass matrix defined by 



"'^Quantitatively, in the scenarios under investigation, values are, respectively, up to order one cor- 
rections, yelVr ^ 10^^ and yl/yl ^ 10^^ times the characteristic mass scale of the mass matrix. 
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= UM^U^ with = diag(mi, m2e^'^'^^ , m^e^'^'^^) and U the neutrino mixing matrix, 
expressed in the PDG parameterization. In this case, 

ruee = mi cos^ 9i2 cos^ 9i3 + ms sin^ 9i2 cos^ ^ige^^'^^ + mg sin^ ^ige^^^'^^"'^). (3.4) 

Setting rUee = 0, one can solve for mi and one of the Majorana phases. Recalhng that, for 
the normal mass hierarchy, 1712 — ml + Am^2 "^3 — \^^rr^^+~Arn^, and assuming 



small and i] = ^/Arn^^/Am 



2 

13; 



V ,/2nn - ^13 _„o^ COS[2(03 - 5)], 



^/A;^ 'cosV2 2^5 "cos 2^5 



02 ~ - + -arctan -M^^^^^sin[2(03-<^)] . (3.5) 
2 2 \ 77 sm 2^5 / 

One can easily obtain approximate expressions for the other neutrino masses (m2, m^) and 
hence all elements map- Upon substituting the numeric best fit oscillation parameters to 
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avoid introducing a needlessly cumbersome expression, I get 



^ 



Am 



2 

13 



0.5e 



1 1 
1 1 



+ O.710i3e-^('^-2'^^) 



1 1 \ 

1 
1 y 



+ 0.45?7 



-1.3 1 
-1.3 -1 0.61 

1 0.61 -0.36 

/ n 



+ O.9ie^^cos[2(5-03)] 



+ 1.2i^^^sin[2(03-5)] 



v 



1 -0.89 

1 0.12 0.02 
-0.89 0.02 -0.12 ^ 

1 -0.67 

1 1.2 -0.83 
-0.67 -0.83 0.56 



(3.6) 



Eq. (13.71) suggests a clear hierarchy among the mixing matrix elements. The four, lower 
box-diagonal /x — r elements dominate, followed by the off-diagonal e/i and e/i entries, 
and finally the vanishingly small meg. Except for the vanishingly small mee, which was 
required a 'priori^ all of the remaining properties follow directly from the experimentally 
determined mixing parameters. Among the dominant /i — r submatrix, Eq. (13.71) predicts 
that all entries are equal up to small order ri and ^13 corrections. The magnitude, and sign, 
of these "breaking terms" can be tuned with the phases 03 and 5, and to a lesser extent 
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by varying rj and 613 within their allowed ranges. On the other hand, the relative sizes of 
rrie^ and are expected to be similar but not identical, i.e., rrie^ ~ 771^^ ~ (mg^ — mer). 

While some of the gross features of Eq. (13.71) are shared by Eq. ( 13.21) and Eq. (13. 3p . 
a finer analysis reveals several disagreements. The major discrepancy lies in the required 
relations among the matrix elements. Eq. ( 13. 2p predicts that all mar elements are equal, 
while Eq. (13.31) suggests rngT -C -C m^^. Both of these contradict, in different ways, 
the experimental constraint <^ Trifj^r ~ rrirT- Additionally, both Eq. ( 13. 2p and Eq. ( 13.31) 
predict niee <^ -C rnrr, while observations require -C ^ nirT- Similarly, both 
sets of operators suggest nie^ <C while, experimentally, they are constrained to be 
similar. 

In order to quantify how much Eq. (13.21) and Eq. (13. 3p (dis) agree with our current un- 
derstanding of neutrino masses and lepton mixing, I numerically scanned the allowed mass 
matrix parameter space assuming the normal neutrino mass hierarchy and constraining 
I'^eel < He/VT X 1 gV ~ 10~^ cV. It should be noted that, according to this relation, 
rriee is allowed to deviate by nearly a factor of ten above naive expectations from mass 
matrix Eq. (13.21) . thus accounting for the possible order of magnitude uncertainties in 
operator scales and coupling constants. This feature is only included for completeness, 
as one expects that such rriee excursions from zero will generally have negligible effect on 
the mass matrix due to the robust nature of Eq. (13. 7p . Fig. 13. 9[ a scatter plot of mixing 
matrix elements, depicts the result of such a scan. Note that I plot the mass ratios with 
respect to assumed-to-be-dominant m^,- element. The light grey regions of the plot were 
produced allowing all oscillation parameters to vary within their 95% confidance bounds 
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[16] and phases to vary within their entire physical range subject to the constraints dis- 
cussed above. In the purple (dark) region, the phases and reactor mixing angle 6*13 are 
allowed to vary while all other mixing parameters are held fixed at their best fit values. I 
depict the sin^ ^13 variation from zero to 0.06 (4cr upper bound [16]) by varying the purple 
shading from dark to light. It is easy to check that the numeric (Fig. 13. 9p and analytic 
results (Eq. (13.71) ) are consistent both qualitatively and quantitatively. 

Fig. 13.91 also depicts the predictions from Eq. (13.21) and Eq. (13.31) with red (closed) and 
blue (open) dots, respectively. As expected, all the predictions from Eq. (13. 3p fall near 
the origin in each panel and are safely excluded. Because expectations from Eq. (13.31) 
for all map/mrT are much smaller than one, I also include the dot coordinate values for 
both textures within the figure. In order to render the neutrino mass matrix predicted 
from (^36, Oar and O^^ consistent with experimental constraints on neutrino masses and 
lepton mixing, one is required to choose very hierarchical A coefficients. In more detail, 
one needs to choose A values so that all mixing matrix elements are enhanced relative to 
the dominant m^T- oc UtUt by numerical factors that range - for different entries - from 100 
to 10^. A possible mechanism for achieving this is to suppress third generation couplings 
to new physics, thus driving up the ratio rriais/'mrT along with the required cutoff scale 
A^. This procedure would have to be accompanied by a more modest reduction of the 
couplings of second generation fermions. Basically, I need to impose a flavor structure 
that "destroys" the naive flavor structure induced by the charged lepton Yukawa coupling 
hierarchy. I can safely conclude that C36, C37, and O^s, which suggest that the neutrino 
mass matrix has the form Eq. (13. 3p . are strongly disfavored by current neutrino oscillation 
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Figure 3.9. Scatter plots of the symmetric Majorana neutrino mass matrix 
elements normalized to m^,-. Each panel is produced assuming the normal 
mass hierarchy and parameter constraints insuring that rriee ^ 10~^ eV. The 
light grey region is calculated allowing all mixing parameters to vary within 
their respective 95% confidence intervals. In the purple (darker) regions, 
the solar and atmospheric parameters are held constant while all phases are 
scanned within their physical ranges and 613 is varied between zero and its 
4 a upper bound. The sin^ ^13 variation is illustrated by varying the shading 
from dark to light. Also indicated by red (closed) and blue (open) dots are 
the expectations derived from Eqs. (13.21) and (13. 3p . respectively, along with 
a listing of their associated coordinate values. 
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data and, if somehow realized in nature, must be accompanied by a very nontrivial flavor 
structure. 

On the other hand, the operators listed in Eq. (13. ip . which predict Eq. fl3.2p . are not 
quite as disfavored. In this case the hierarchies among different mass matrix elements are 
softer, and one can ask whether the red dots in Fig. 13.91 can move toward the experimen- 
tally allowed regions with order 1-10 relative shifts. Many of the predictions are already 
in agreement with experimental constraints, or at least close enough to be easily "nudged" 
toward acceptable levels with order one coefficients. The figure reveals that only m^^ is 
predicted to be relatively too small. By enhancing it by a factor of order yr/Un ~ 20 one 
obtains moderately good agreement between Eq. ( 13. 2p and experimental requirements. I 
therefore conclude that operators listed in Eq. (13. ip are at least marginally allowed by 
neutrino mixing phenomenology. 

While essential for a complete understanding of neutrino masses and mixing, improved 
measurements of the already determined mixing angles and mass-squared differences will 
not help to further constrain/exclude any of the LNV scenarios in question. Considering 
our parameter flexibility, only future neutrino experiments that provide qualitatively new 
results can aid in this endeavor. In particular, the experimental determination of the 
neutrino mass hierarchy is essential in order to properly test the scenarios highlighted 
in this section, as they all predict, in the absence of very non-trivial flavor structure in 
the LNV sector, the normal hierarchy. Next-generation neutrino oscillation experiments 
are expected to provide non-trivial information regarding the neutrino mass hierarchy. 
Most rely on a neutrino/anti-neutrino oscillation asymmetry via Earth matter effects 
|14|. 161] . and depend heavily on a sufficiently large Oi^ mixing angle. The possibility 
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that 6'i3 is vanishingly small, where the standard approach is ineffective, is addressed 
in |61j considering both oscillation and non-oscillation probes. In that case, one can 
hope to discern the neutrino mass spectrum in future neutrino factory [S6\ / Superbeam 
[35] experiments coupled with improved constraints on the effective masses extracted 
from tritium beta decay [23j and cosmology |19L I20L I21L \22\n Note that these non- 
oscillation probes can be independently used to constrain LNV models, as they provide 
information regarding the the magnitude of the lightest mass eigenstate (mi [m^] in the 
case of normal [inverted] hierarchy). For example, if either cosmological observations or 
tritium beta decay experiments see evidence for non-zero neutrino masses (in more detail, 
they constrain S = Xli"^* ml^ = X^i^^il^eiP respectively) such that S ^ 0.05 eV 
or m^^ ^ 0.01 eV, one would conclude, assuming a normal mass hierarchy, that mi ^ 



-y/Amf^. This would destroy the possibility of negligibly small rriee, and hence disfavor 
the operators that lead to mass matrices of the type Eq. (13. 2p and (13. 3p . Currently, S 
and ruu^ are bounded to be below 0.94 eV and 2.0 eV, respectively, but the sensitivity to 
these observable is expected to significantly improve with next-generation experiments to 
0.1 eV [62J and 0.2 eV [63], respectively. 

3.4. Phenomenologically interesting operators: Sample Renormalizable 

Model 

Having superficially surveyed a large set of LNV operators, I am now in a position to 
identify operators with "interesting" phenomenological features for further detailed study. 



^One traditionally includes the effective P(30iy mass rriee given by Eq. (|3.2p in a neutrino mass hierarchy 
analysis. However, as discussed in Sec. 13.2. II m^f is a potentially convoluted process-dependent quantity 
that generally has little (directly) to do with neutrino masses. For this reason, PPOiy constraints cannot 
be used to determine the neutrino mass spectrum from the point of this analysis. 
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One subset of potentially interesting operators is characterized by those that, when re- 
quired to "explain" the observed neutrino masses, are accompanied by a low cutoff scale of, 
say, less than several TeV. Further requiring a small enough ml^ in order to evade current 
PPOu constraints, this set contains only seven elements: On, Ois, C34, C35, O56, C57, Css- 
Of these, all but operators O^^ and O^s (which lead to the zeroth-order neutrino mass 
matrix Eq. (13.21) and a suppressed mee) should provided a positive LNV signal in the next 
round of double-beta decay experiments, baring specific flavor symmetries or finely-tuned 
couplings. Furthermore, O^q leads to a (3l3Qv rate that is higher than what is naively 
dictated by the values of the neutrino masses. Finally, with the possible exception of 
C56 which may mediate observable LNV processes at high energy colliders, none of the 
seven operators above are expect to mediate LNV violating phenomena (as defined here) 
at accessible rates. 

An "orthogonal" subset consists of the higher dimensional operators already "ex- 
cluded" by (3(3Qv. Not including those operators severely constrained by lepton mixing in 
Sec. 13.31 this hst contains 11 elements: Oie, O19, O53, 054a,fe,c,d, C'sg, C'eo, C^ro, C^75- Most of 
these are associated to cutoff scales of order the weak scale, which are likely to already be 
constrained by different searches for new degrees of freedom with masses around 100 GeV. 
Even if those are considered to be excluded, Cig, C19, C75 are "safely" shielded from direct 
and indirect non-LNV searchesjll while still mediating potentially observable LNV effects 
at colliders as long as the new physics does not couple, to zeroth order, to first generation 
quarks (in order to evade the /?/50z/ constraints). 

Generic new degrees of freedom at the weak scale are constrained by direct and indirect searches at high 
energy colhders {e.g., resonances and effective four-fermion interactions, respectively), flavor-violating 
{e.g., jjL 67), and high precision experiments {e.g., measurements of the anomalous muon magnetic 
moment). 
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Regardless of whether these different options for the LNV sector lead to observable 
LNV phenomena, the low extracted cutoff scale of all the operators highlighted above 
implies that new degrees of freedom should be produced and, with a little luck, observed 
at the LHC or, perhaps, the ILC. Furthermore, the TeV scale has already been identified 
as an interesting candidate scale for new physics for very different reasons, including the 
dark matter puzzle and the gauge hierarchy problem. The fact that, perhaps, the physics 
responsible for neutrino masses also "lives" at the TeV scale is rather appealing. 

In order to study this new physics, as already emphasized earlier, ultraviolet complete 
manifestations of the physics that leads to the effective operators are required. Here 
I discuss one concrete example. Other examples (for different effective operators) were 
discussed in |15] . Given a specific LNV operator, it is a simple matter to write down 
equivalent renormalizable Lagrangians. I briefly illustrate this procedure by constructing a 
renormalizable model that will lead to the dimension-eleven operator O^q. It is among the 
interesting LNV effective operators of the sample highlighted above, since it is currently 
unconstrained by /3/30z/ searches regardless of the quark-fiavor structure of the operator, 
while rrigf ^ niee for [3(3Qv. On the other hand, for O^^ is very low (below 500 GeV), 
so that the new degrees of freedom may already be constrained by, for example, Tevatron 
or LEP data. I will not worry about such constraints henceforth, but will only comment 
on possible phenomenological problems. 

O56 can be accommodated by a wide variety of models, as can be seen from its possible 
Lorentz structures. In terms of scalar/tensor helicity-violating bilinears = l,cr^^, and 
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vector helicity-conserving bilinears = 7^, these are 

C»56 = {{L%Q^){dT,d'){WT,-^), (LT,Q^)(rfT,^)(dT,i^), {LTJ'^){Q^TJ')(d-T,-^), 
{VT^d') {Q'TJ-) (rfTei^) , {UT,d') {Q'T,W) {dTJ-) , {UVj^) {Q^T.d') (dT,i^) , 
(LT,^) (g^ r,i^) {d-VJ') , (LT,i^) iQ'TJ') (dTJ-) , (LT,i^) {Q^TJ-) (rfT^d^) } 
X H^H'eike.i. (3.1) 

It is clear from the chiral field content that these operators depend on combination 
of helicity-conserving and helicity-violating interactions. In particular, it is impossible to 
form any of the operators in this long list with only the addition of vector boson states: 
new heavy scalar and/or tensor particles are probably required if O^^ is the proper tree- 
level manifestation of the LNV physics at low-energieso Furthermore, the couplings of the 
new physics fields with one another must be constrained in order to "block" the presence 
of lower-dimensional tree-level effective operators. This usually implies the existence of 
new exact (broken) symmetries to forbid (suppress) particular interactions. 

Certain Lorentz structures, those containing only bilinears, can be realized assum- 
ing that the LNV ultraviolet sector contains only heavy scalar fields and I concentrate, 
for simplicity, on this possibility [64j . Simple scalar interactions that can lead to O^^ are 
shown in the diagram in Fig. 13.101 Specifically, these yield the effective operator Lorentz 
structure {UQ^){d'^d'^){e'^d'^)H^ eik^ji with the introduction of four charged scalar fields, 
01) 02, 03, 04- The gauge structure is such that, under (5'f/(3)c, SU{2)l, U{1)y]u, 4>i trans- 
forms as a (3, 3, +1/3), 02 as (3, 1, —2/3), 03 as (3, 1, —4/3), and 04 as (3, 1, —2/3). While 



^Other possibilities include heavy vector-like fermions. 

"'in the case of U{1)y, 'transfoms as X' means 'has hypercharge X\ 
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Figure 3.10. Sample scalar interactions that lead to the "interesting" effec- 
tive operator O^q with the Lorentz structure {UQ^){d'^dF){e'^d'^)H^H^eikeji. 



02 and 04 have identical gauge quantum numbers, they have different baryon number (2/3 
versus —1/3). 0i has baryon number —1/3, while 03 has baryon number 1/3. Lepton 
number cannot be consistently assigned as it is explicitly violated by two units. 

04, which does not couple to any of the SM fermions, plays an essential role. It acts 
as a selective "insulator" that connects the various interaction terms in such a way as to 
only alow certain tree-level higher dimensional SM effective operators. All renormalizable 
theories that lead to only very high dimensional effective operators contain one or more 
of these "hidden sector" fields. Note that the new scalar fields should not acquire vacuum 
expectation values in order to avoid the presence of lower dimensional irrelevant operators 
that are likely to dominate low-energy phenomenology and - much more important - to 
prevent the spontaneous breaking of color or electromagnetic charge. 
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Given the scalar field content as well as its transformation properties under SM global 
and local symmetries, it is a simple matter to write down the minimal interaction La- 
grangian density for the system. A candidate renormalizable Lagrangian is 

i 

+ A234M020304 + h.C. . (3.2) 

Each term in Eq. (13.21) . including those involving covariant derivatives Z)^, is implicitly 
assumed to respect the gauge representations of the associated (pi fields, as defined above. 
The Yukawa- type couplings tji, as well as the Aj scalar vertices are dimensionless, and 
assumed to be of order one, while I assume all scalar masses Mj to be of the same order 
of magnitude. In this case, A ~ Mi. In the A234 term, an overall mass scale M has been 
"factored out" and is assumed to be of the same order as the Mj. Note that I neglect 
generation indicies, which are implied. In the case A234 = 0, lepton number is a classical 
global symmetry of Eq. (13.21) . and one can view this three-scalar coupling as the source 
of lepton number violation. One may even envision a scenario where lepton number is 
spontaneously broken by the vacuum expectation value of some SM singlet 05 scalar field, 
(05) = M. 

Provided all Mj are around 0.5 TeV, as required if this Lagrangian is to "explain" 
the observed light neutrino masses, LNV is certainly not the only (or even the main) 
consequence of this model. The yi and terms for example, will mediate /z — e- 
conversion in nuclei at very dangerous levels if their flavor structure is generic. 02 can be 
resonantly produced in (i(i-collisions, while 0i and 03 qualify as scalar lepto-quarks, which 
are constrained by high energy collider experiments, including those at HERA [65] , to 
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weigh more than a few hundred GeV [34]. For more details, I refer readers to, for example, 
the Particle Data Book [34] and references therein. 

I will conclude this discussion by adding that several other effective operator can be 
realized in a very similar way. Oig, for example, if it manifests itself with the Lorentz 
structure {UQ^){d'^d'^){e^u^)eij, can be realized by a Lagrangian very similar to Eq. ( 13. 2p 
where the field in the j/3-coupling interaction is replaced by a u'^ field, and the (pi field 
is replaced by an SU{2)l singlet (it is a triplet in Eq. (13.21) ). Of course, hypercharge 
assignments for the 0j also need to be modified in a straight forward way. The associated 
non LNV phenomenology is similar, except for the fact that for Oig (around 1 TeV) 
is larger than the one for O^^ and hence Oig is less constrained by current experimental 
data. On the other hand, Oig predicts potentially much larger rates for LNV observables 
at coUiders (see Fig. 13. 7p . 

Our definition of "interesting" is arbitrary and motivated only by the fact that the 
physics of the "interesting" operators highlighted earlier in this section will probably be 
explored at next-generation collider and high-precision experiments. One may argue that 
many operators which lead to the observable neutrino masses for high values of are 
interesting on their own right, either due to the theoretically pleasing properties of their 
associated potential ultraviolet completions, or by some observational peculiarity. There 
are many examples of the first type ranging from the different manifestations of the seesaw 
mechanism [7], I25L I26L 127] to the Zee model |66] and the minimal supersymmetric SM 
with R-parity violation [67] . Dedicated analysis of these cases have been widely pursued 
in the literature and will not be discussed here. I would also like to point out that some 
effective operators, like Oj and Og, are, according to our criteria, very "uninteresting." 
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Both O-j and Og predict unobservably suppressed /5/?0i/ rates (both predict small mee) 
and equally hopeless collider prospects given that they are associated to very high cutoff 
scales, Ai. «i 4 X 10^ TeV and A,^ ^ 6 x 10"^ TeV, respectively. If either of these operators 
are responsible for the observed tiny neutrino masses, it is quite possible that I may never 
directly detect LNV. It is curious to consider possible means of indirect detection or other 
observable consequences of the different ultraviolet completions of such scenarios, 3 It 
would also be interesting to ask whether either of these elusive models has any underlying 
theoretical motivation or whether they allow one to solve other outstanding problems in 
particle physics. 

3.5. Conclusion 

If neutrino masses are a consequence of lepton-number violating physics at a very 
high energy scale (higher than the scale of electroweak symmetry breaking), new physics 
effects - including the generation of neutrino Majorana masses - at low enough energies 
can be parameterized in terms of irrelevant operators whose coefficients are suppressed 
by inverse powers of an effective cutoff scale A. As discussed before, A is, roughly, the 
energy scale above which new degrees of freedom must be observed if the new ultraviolet 
physics is perturbative (if the new physics is very weakly coupled, the masses of the new 
degrees of freedom can be much smaller than A). I have explored a very large class of 
such scenarios through 129 irrelevant operators of energy dimension less than or equal to 
eleven that violate lepton number by 2 units. These are tabulated in first two columns of 



Table 13.11 along with a summary of our results. 



^This is very similar to the case of 0\. The main redeeming feature of Oi, other than its simphcity, is 
the fact that many of its ultraviolet completions allow one to explain the matter-antimatter asymmetry 
of the universe 
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Analyzing each effective operator individually, I estimated the predicted general form 
of the Majorana neutrino mass matrix. Our results are listed in the third column of 
Table I3.1[ By comparing each such estimate with our current understanding of neutrino 
masses, I extracted the cutoff scale K^, of each effective operator, assuming that it provides 
the dominant contribution to the observed neutrino masses. These results are listed in 
the fourth column of Table 13.11 assuming light neutrino masses equal to 0.05 eV (the 
square root of the atmospheric mass-squared difference), and are summarized as follows. 
Depending on the field content and dimension of the irrelevant operator, the "lepton 
number breaking scale" Ky is predicted to be anywhere from the weak scale (~ 0.1 TeV) 
all the way up to 10^^ TeV (see Fig. 13. 2p . This means that, depending on how lepton 
number is violated and communicated to the SM, the mass of the associated new degrees 
of freedom is predicted to be anywhere between 100 GeV and 10^^ TeV, even if all new 
physics couplings are order one. I note that in the case of all variations of the seesaw 
mechanism (Oi), neutrino physics constrains = 10^^ TeV such that the new degrees 
of freedom are either unobservably heavy, extremely weakly coupled, or their couplings 
to the SM degrees of freedom are finely-tuned. It is fair to say that this behavior is not 
characteristic of all LNV ultraviolet physics. One sample ultraviolet theory that leads to 
dimension-eleven LNV effective operators was discussed in Sec. 13. 4[ Other examples can 
be found in [15] . and include supersymmetry with trilinear R-parity violation and the 
Zee model. 

Assuming that a particular operator is responsible for nonzero neutrino masses, it is 
straight forward to ask whether it leads to other observable consequences. Here, I con- 
centrated on several LNV observables, and included future LNV searches at the LHC and 
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future lepton machines (like the ILC), along with their ability to directly produce (and 
hopefully observe) new physics states lighter than several hundred GeV. In column five 
of Table 13.11 I list the most favorable modes of experimental observation for each opera- 
tor. The different relevant probes are: neutrinoless double-beta decay (/5/?0z/), neutrino 
oscillation and mixing (mix), direct searches for new particles at the LHC (LHC) and 
ILC (ILC), and virtual LNV effects at collider facilities (HElnv). I find it unlikely that 
other probes of LNV, including rare meson decays, should yield a positive signal in the 
forseeable future. This conclusion is strongly based on the fact that, for all of our analysis, 
I assume that all new physics degrees of freedom are heavier than the weak scale. While 
the vast majority of operators is most sensitive to searches for neutrinoless double-beta 
decay, that is not true of all operators. Some lead to relatively suppressed rates for PPOu 
(mostly because they lead to mass matrices with a very small rriee) even if they are asso- 
ciated to A,^ < 1 TeV, indicating that, for these scenarios, we are more likely to observe 
the physics behind neutrino masses directly at colliders than to see a finite lifetime for 
/?/30z/. Other scenarios naively lead to (3(3Qv rates orders of magnitude higher than what is 
currently allowed by data. If these are responsible for the generation of neutrino masses, 
the new physics is constrained to be somewhat decoupled from first generation quarks 
(for example). In this case, there is hope that LNV phenomena at colliders, which are 
not restricted to first generation quarks, occur with non-negligible rates. 

The sixth column of Table 13.11 lists the current "status" of the operator as either 
experimentally unconstrained (U), constrained (C), or disfavored (D). Such labels are 
assigned based only on the experimental probes reviewed in this work. By arbitrary 
convention, an 'unconstrained' operator can safely accommodate all existing data even 
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if one assumes all its flavor- dependent dimensionless coeflicients to be of order one. A 
'constrained' operator can accommodate all existing data after one allows some of the 
different flavor-dependent dimensionless coefficients to be suppressed with respect to the 
dominant ones by a factor of 100 or so (as described above). 'Disfavored' operators can 
only accommodate all data only if "tuned" much more severely than the 'constrained' 
ones, and are usually in trouble with more than one "type" of constraint. A glance at 
column six reveals that 11 out of the 129 operators are disfavored by current data. The 
most stringent constraints come from /?/50z/, while all 'disfavored' operators are associated 
to cutoffs at or below 1 TeV. Three of the 'disfavored' operators, Osq, O^r, and Oss, are 
also in disagreement with the neutrino oscillation data (see Sec. 13. 3p . 

Our results illustrate that, as far as "explaining" neutrino masses, the model-building 
scene is wide open even if one postulates that neutrino masses arise as a consequence 
of lepton number violating, "heavy" physics. Significant progress will only be achieved 
once more experimental information becomes available. The observation that neutrinoless 
double-beta decay occurs with a nonzero rate will help point us in the right direction, 
but will certainly not reveal much about the mechanism behind neutrino masses. A 
more complete picture can only arise from combined information from several observables, 
including other LNV observables and the search for new physics at the electroweak scale. 
Other important experimental searches, not discussed here, include all lepton-number 
conserving "leptonic" probes, such as precision measurements of the anomalous magnetic 
moment of the muon, searches for leptonic electric dipole moments, searches for charged- 
lepton fiavor violation (see [69] for a model independent discussion of this issue), and 
precision measurements of neutrino-nucleon and neutrino-lepton scattering. 
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CHAPTER 4 

Neutrino Phenomenology of Very Low-Energy Seesaws 

The seesaw mechanism [7] is an appeahng way to generate the observed neutrino 
masses and lepton mixing matrix. The idea is simple. Add an arbitrary number of 
singlet fermion states to the SM matter content. The triviality of their quantum numbers 
allows them to have Majorana masses of magnitude M, as well as couple to the SU(2)l 
lepton and Higgs doublets. The latter vertices become Dirac mass terms of magnitude 
H after electroweak symmetry breaking. The standard theoretical prejudice is that the 
Dirac masses are of order the charged fermion masses, while the Majorana masses are 
at some very high energy scale, M ^ 100 GeV. If this is indeed the case, the resulting 
propagating neutrino degrees of freedom separate into two quasi-decoupled groups: mostly 
active states with very small masses m ~ /i^/M suppressed by the new physics scale, and 
mostly sterile states with very large masses M. In this scenario, the mostly right-handed 
states are not directly observable. Indeed, it is possible that if such a high-energy seesaw 
is realized in Nature, its only observable consequence is that the mostly active neutrinos 
have mass and mix. 

Of course, there is no direct evidence that M - which I refer to as the seesaw scale 
- is large. Large M values are attractive for several reasons, including the fact that one 
may relate M to the grand unified scale. On the other hand, all M values are technically 
natural, given that when M vanishes the global symmetry structure of the Lagrangian 
is enhanced: U{1)b-l is a symmetry of the Lagrangian if M = 0, so that M is often 
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referred to as the lepton number breaking scale. This point was recently emphasized in 
[32] . Recent analyses have also revealed that there are several low-energy choices for the 
seesaw energy scale that allow one to address outstanding problems in particle physics 
and astrophysics. The main reason for this is that, unlike in the high-energy seesaw, 
in a low-energy seesaw the mostly right-handed states do not decouple but remain as 
kinematically accessible sterile neutrinos. 

The data reported by the LSND short baseline neutrino oscillation experiment [70] can 
be explained by postulating the existence of light (m ~ 1 — 10 eV) sterile neutrino states. 
This result is currently being tested by the Fermilab MiniBooNE experiment |71j and, if 
confirmed, will lead the community to seriously contemplate the existence of light, SM 
singlet fermions. It was pointed out in [32] (see also [72j ) that if M ~ 1 — 10 eV, the right- 
handed seesaw neutrinos can easily play the role of the LSND sterile neutrinos. There is 
also evidence for mixed sterile neutrinos at other energy scales: eV sterile neutrinos aid in 
heavy element nucleosynthesis in supernovae, keV sterile neutrinos can help explain the 
peculiar velocity of pulsars, and remain viable warm dark matter candidates. In the past 
several months, it has been shown that the seesaw right-handed neutrinos may play the 
role of all these astrophysically/cosmologically inspired sterile neutrinos |31L 173] . 

In this chapter, I consider the phenomenology of low-energy (M < 1 keV) seesaw 
scenarios, extending the analysis performed in [32] in several ways. In Sec. 14.11 1 review the 
generation of neutrino mass via the seesaw mechanism and apply it to relatively light right- 
handed neutrino states. I pay special attention to the most general active-sterile seesaw 
mixing matrix, whose parameters are the main object of our study. In Sec. 14.21 1 review the 
several different "evidences" for sterile neutrinos, and discuss whether these can be "fit" 
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by the low-energy seesaw. I concentrate on exploring solutions that can accommodate at 
the same time the LSND anomaly and the astrophysical processes outlined above, but 
also discuss different combinations of the seesaw parameters capable of explaining only 
the astrophysics-related observables. In Sec. 14.31 I examine other experimental probes 
that can be used to explore low-energy seesaws - regardless of their relationship to the 
LSND anomaly, pulsar kicks, and warm dark matter. I concentrate on the prospects 
of current /future tritium beta-decay and neutrinoless double-beta decay experiments. I 
conclude in Sec. 14.41 by summarizing our results, commenting on the plausibility of this 
scenario, and offering a general outlook for the future. 

4.1. The Seesaw Mechanism and electron volt neutrino masses: preliminaries 

In order to account for nonzero neutrino masses, I add to the SM particle content 
three SU{3)c x SU{2)l x U{1)y gauge singlet fermion states Ni, conventionally referred 
to as right-handed neutrinos. While sterile under SM gauge interactions, right-handed 
neutrinos may still be charged under new, currently unknown gauge transformations. 
Such interactions, if at all present, are neglected in our analysis. 

The most general renormalizable Lagrangian consistent with SM gauge invariance is 



where £oid is the traditional SM Lagrangian, H is the Higgs weak doublet, L", a = 
e, /X, r, are lepton weak doublets, Xai are neutrino Yukawa couplings, and Mj are Majorana 
masses for the Ni. I choose, without loss of generality, the Majorana mass matrix Mr 
to be diagonal and its eigenvalues Mj to be real and positive. I also choose the charged 
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lepton Yukawa interactions and the charged weak current interactions diagonal so that 
all physical mixing elements are contained in the neutrino sector. 

After electroweak symmetry breaking (when H develops a vacuum expectation value 
f ), Cy will describe, aside from all other SM degrees of freedom, six neutral massive Weyl 
fermions — six neutrinos. The resulting mass terms can be expressed as: 



where /i = Af , and v and are vectors of the three active neutrinos {veiV^iVr) and the 
three right-handed, sterile states, respectively. Each entry in the symmetric mass matrix 
of Eq (14.21) is itself a 3 x 3 matrix of mass parameters. Diagonalization of the mass matrix 
yields eigenstates with Majorana masses that mix the active-active states, related via the 
standard lepton mixing matrix V ^ and the active-sterile states. The physical neutrinos 
are thus linear combinations of all active and sterile states. Throughout I will work in 
the "seesaw limit," defined by <^ Mr. In this case, there are three mostly active light 
neutrinos and three mostly sterile heavy neutrinos where 'mostly' is determined by the 
induced mixing parameters. 
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In the seesaw limit, the diagonahzation is simple. Assuming, for simplicity, that the 
mixing matrices are real: 



fi 



1^^ Mr 
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X 



1 e 

-e^ 1 
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(4.3) 
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+ O(0' 
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where m is the diagonal matrix of light neutrino masses and is a matrix of active-sterile 
mixing angles found from the relations 



QMr = /i. 



QMrQ'^ 



-VmV' 



(4.4) 
(4.5) 



The elements of are small {0{^/ Mr)), and the standard seesaw relation {VmV'^ = 
— jjiM^^ jji^) is easily obtained by combining Eqs. (14.41) and Eqs. (14. 5p . On the other hand, 
using Eq. (14. 5p . I can relate the mixing parameters in to the active-active mixing 
angles contained in V and the neutrino mass eigenvalues. In the case of observably light 
sterile neutrino masses, as considered in our analysis, this equation is very useful, as it 
places testable constraints on observable quantities. The general solution (first discussed 
in detail in [74]) of Eq. ( l45l) is 



e = -Vm}'^OM~^''^, 



(4.6) 
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where O is an orthogonal matrix parameterized by three mixing angles 0i2, 0i3, 023 
Physically, the mixing matrix O is a consequence of our freedom to choose the form of 
Mr. An illustrative example is found by considering the mass ordering of M/j in its 
diagonal form. The reordered matrix Mfi{mi mj) is equivalent to a 7r/2 rotation in 
the i — j plane and therefore represents the same physics as the original matrix, as it 
should. In other words, OMptO^ is the physically relevant object, as opposed to O and 
Mr separately. This object, when constrained to be real, has six free parameters that I 
shall refer to as "heavy parameters:" 0i2, 0i3, 023, Mi, M2, M3. 
Using Eq. (14.61) . the 6x6 unitary neutrino mixing matrix is 



U 



(4.7) 



Note that, up to (9(G^) corrections, V is unitary. U is entirely described in terms of the 
three light mixing angles, six mass eigenvalues and three angles 0jj. Many combinations 
of these have been measured or constrained via oscillation searches, cosmology and as- 
trophysics. In particular, the two active neutrino mass-squared differences and mixing 
angles have been measured j75L I76L 134] , thus leaving free the six heavy parameters and 
the absolute scale of active neutrino masseso With f/, the corresponding neutrino mass 
values and the SM couplings I can calculate all observable quantities and compare them 
with data. 



In general, O is a complex orthogonal matrix. Here, however, I will restrict our analysis to real neutrino 
mass matrices, unless otherwise noted. 

^The active neutrino mass hierarchy, normal vs. inverted, is another (discrete) free parameter. 
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It is natural to wonder how well the seesaw approximation holds once one starts to 
deal with M/j values around 1 eV. From Eq. fl4.4p . it is clear that one can choose for 
the expansion parameter ^JmjM^. In all scenarios considered here, ^Jui/Mr < 0.5 (for 
M ~ 1 eV and m ~ 0.3 eV). In the worst case scenario, therefore, first order corrections are 
55% of the leading order terms, while second order corrections are near 30%. Corrections 
to most observables of interest are much smaller than this because they are suppressed by 
larger right-handed neutrino masses. The first non-trivial correction to Eq. (14.61) occurs 
at second order and I find that, for the ambitions of this chapter, all approximations are 
under control. This simple argument has been verified numerically for the most worrisome 
cases. 

Before concluding this section, I wish to add that operators that lead, after electroweak 
symmetry breaking, to Majorana masses for the left-handed neutrino states (Ml) are also 
allowed if one introduces SU{2)i Higgs boson triplets or nonrenormalizable operators 
to the SM Lagrangian. While I neglect these "active" Majorana masses, I caution the 
reader that the existence of such terms would alter our results significantly. In particular, 
assuming the seesaw approximation holds, Eq. (14. 5 p would read 

VmV^ + QMrQ^ = Ml, (4.8) 

which leads to a relationship between G, m, and different from Eq. (14. 6p . If this were 
the case, for example, it would no longer be true that the largest value (in absolute 
value) is constrained to be smaller than (mmax/A^^R.min)^''^, where mmax is the largest 
element of m, while Mr min is the smallest element of M/j. On the other hand, all objects 
on the left-hand side of Eq. (14. 8 p are observables. Hence, in the case of a low-energy 
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seesaw, one can expect, in principle, to be able to test whether there are contributions to 
the neutrino mass matrix that are unrelated to the presence of right-handed neutrinos. 
By measuring V, m, Mr, and 0, one can establish whether Ml is consistent with zero. 



4.2. Oscillation Phenomenology and Current Evidence for Low-Energy 

Seesaw 

Here I examine a number of experimental and observational anomalies that may be 
explained by light sterile neutrinos. More specifically, I explore what these can teach us 
about the currently unknown parameters of the seesaw Lagrangian, described in detail 
in Sec. 14.11 In all cases I assume 3 mostly active and 3 mostly sterile neutrinos and, 
most of the time, will concentrate on a 3 + 2 + 1 picture of neutrino mass eigenstates, 
that is, three mostly active sub-eV neutrinos, two mostly sterile eV neutrinos and one 
almost completely sterile keV neutrino. The hope is that the heavier state can account for 
warm dark matter (section 14.2.21) or pulsar kicks (Sec. 14.231) . which both require at least 
one keV neutrino, while the other two mostly sterile states help "explain" the existing 
oscillation data where, for all practical purposes, the heaviest neutrino decouples and I 
are left with an effective 3 + 1 or 3 + 2 picture. I remind readers that a third possibility 
(2+2|3) is currently ruled out by solar and atmospheric data [771 I78L 179] and will be 
ignored. 3 + 1 schemes that address the LSND anomaly are also disfavored by global 
analysis of short baseline oscillation experiments \77i I78L I79L 180] and, for this reason, I 
mostly concentrate on 3 + 2 fits to the LSND anomaly [80j . 



^It would have been rather difficult to construct a 2+2 neutrino mass hierarchy using the seesaw La- 
grangian. 
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Our analysis method is as follows: For each experimental probe considered I perform 
a "fit" of the mixing matrix f/, given by Eq. (14.71) . and neutrino masses to the "data", 
and extract the region of parameter space that best explains the data. In most cases I 
allow the light mixing angles and mass squared differences to vary within their la limits 
(according to [76j)|j the angles (pi to vary unconstrained within their physical limits of 
— 27r, and the lightest active neutrino mass eigenvalue mi to vary unconstrained between 
— 0.5 eV. The quotation marks around "fit" and "data" are meant to indicate that our 
methods are crude, in the sense that I am fitting to previously processed experimental 
data, assuming a diagonal correlation matrix with Gaussian uncertainties. In order to 
avoid the subtleties involved in such a "fit to a fit" , I hesitate to mention actual confidence 
intervals, but are compelled to do so for lack of a better measure of an allowed region. I 
sometimes present our best fit parameter points along with confidence intervals, but warn 
the reader to avoid strict interpretations of these numbers. While crude, our methodology 
of error analysis and fitting provides a very useful instrument for identifying whether (and 
how) the low-energy seesaw can accommodate a particular combination of data sets. 

Before proceeding, it is useful to cement our notation. The neutrino masses will be 
ordered in ascending order of magnitude from mi to mg in the case of a normal active 
neutrino mass hierarchy {m\ — mf < m\ — mf), while in the case of an inverted mass 
hierarchy they are ordered < mi < m2 < m^ < < mg (in this case |m| — 
m\\ > m\ — m\). The states with masses mi^2,3 are mostly active, while those with 
masses m4 5 Q are mostly sterile. Elements of the mixing matrix are referred to as Uai-, 
where a = e, /i, r, si, S2, S3 (s's are the right-handed neutrino degrees of freedom) and 



^In the case of 3+1 "fits" (c/. Eqs. (|4.10|4.1ip ). 1 kept the active neutrino parameters fixed at their 
best-fit values. 
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i = 1,2,3,4,5,6. I also define Am^j = — m| and will refer to the lightest active 
neutrino mass as m;. In the case of normal (inverted) active neutrino mass hierarchy 
mi = mi {mi = m^). 

4.2.1. Short baseline oscillation constraints 

Here I analyze the constraints imposed on the unknown mixing parameters by current 
neutrino oscillation data. I will assume that all solar, reactor, long-baseline and atmo- 
spheric data are properly fit with active-active oscillations, and that constraints on the 
other seesaw parameters will be provided mostly by short-baseline accelerator experi- 
ments. It is interesting to note that the inclusion of the angles (pij introduces enormous 
freedom into the system. Any one active-sterile mixing angle contained in can al- 
ways be set to zero by an appropriate choice of O. In fact, all but three elements may 
be set to zero simultaneously, with only a single non-zero element in each row and col- 
umn. In this case, these are constrained to be around ^Jm]JM^^i where i is the column 
of the non-zero element. This is especially true when the mostly active neutrino masses 
are quasi-degenerate. An important "sum rule of thumb" is the following. For a given 
right-handed neutrino mass Mj, the active-sterile mixing angle squared is of order m/Mj, 
where m is a typical active neutrino mass. One can always choose parameters so that, for 
at most two values of a = e, /x, r, Uai are abnormally small. In that case, however, the 
"other" Uoci is constrained to saturate the bound \Ua.i^ < mi /Mi. 

The most compelling evidence for light sterile neutrinos comes from the short baseline 
oscillation experiment by the Liquid Scintillator Neutrino Detector (LSND) collaboration 
at Los Alamos. Using a ~ 30 MeV 17^ beam they observed a better than 3(7 excess 
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of T/g-like events above their expected background at their detector some 30 m away 
from the production point |70] . This evidence of ^ Fg oscillation requires a mass- 
squared difference greater than 1 eV^, clearly incompatible with the small mass-squared 
differences observed between the active neutrinos. Several mechanisms, such as CPT- 
violation [811 177] . Lorentz invariance violation [82], quantum decoherence [83], sterile 
neutrino decay [84] and, of course, oscillation into sterile neutrinos have been proposed 
to explain this result. Here I concentrate on the last possibility. 

In order to take into account all short baseline data I "fit" our mixing parameters and 
masses to the results of the 3 + 2 performed in [80] , which are summarized in Table 14.11 
[85] . Here I assume that the heaviest, mostly sterile state does not participate effectively 
in LSND oscillations. This is guaranteed to happen if mg > 10 eV. On the other hand, 
|f/a6p are partially constrained by our attempts to accommodate LSND data with seesaw 
sterile neutrinos, as will become clear in the next subsections. 

I find that m;, the lightest neutrino mass, is constrained to lie between (0.22 — 0.37) eV, 
with a "best fit" value of 0.29 eV. Thus the active neutrino mass spectrum is predicted 



Table 4.1. Parameter values used in our analysis. These were extracted 
from a fit to all short baseline neutrino oscillation experiments including 
LSND within the 3 + 2 scenario [801 185] . la indicates a rough estimate of 
the 1 sigma allowed range for the different parameters. 





Ue4 




Ue5 




Am|, (eV^) 


Aml^ (eV^) 


Central Value 


0.121 


0.204 


0.036 


0.224 


0.92 


22 


la 


0.015 


0.027 


0.034 


0.018 


0.08 


2.4 
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to be quasi-degenerate. A sample 6x6 neutrino mixing matrix that fits all data is 



3+2 



V 



0.8301 


0.5571 


0.001365 


0.1193 


-0.009399 


-0.006513 


-0.3946 
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0.7072 


0.2016 


0.2262 
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-0.0949 
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0.09514 


-0.4792 


1 








0.1343 


-0.1832 


-0.09284 





1 





0.004963 


0.004295 


-0.001268 








1 



(4.9) 



while the associated masses are mi ^ m2 — = 0.28 eV, 7714 = 1.0 eV, = 4.7 eV, 

and mg = 6.4 keV. Note that the matrix in Eq. (14.91) is only approximately unitary, up to 

corrections of order 25%. This result agrees qualitative with those obtained in |32j . The 

neutrino masses and mixings obtained in this "fit" are depicted in Fig. 14.11 I will use the 

results of "fits" similar to this one throughout the chapter. |-| 

One can also aim at a (currently disfavored) 3 + 1 LSND fito In this case, much 
lower mi values are also allowed, extending well into the hierarchical spectrum range. In 
this case, all mi values above 0.01 eV and 0.03 eV are allowed, assuming an inverted and 
normal mass hierarchy, respectively. This is to be compared with the results found in 
[32] . where only trivial choices for O were considered. Examples that "fit" all oscillation 
data include, for an inverted active mass hierarchy: mi ~ m2 = 0.066 eV, ms = 0.043 eV, 



'This is easily accomplished by requiring 7715 > 10 cV. 
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Figure 4.1. Neutrino mass eigenstate spectrum, along with the flavor com- 
position of each state. This case accommodates all neutrino oscillation data, 
constraints from r-process nucleosynthesis in supernovae, and may help ex- 
plain anomalous pulsar kicks (see text for details). While I choose to depict 
a normal hierarchy for the active neutrino states, an inverted active neu- 
trino mass hierarchy would have yielded exactly the same physics (as far as 
the observables considered are concerned). 
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For a normal mass hierarchy, I find that mi = 0.055 eV, m2 = 0.056 eV , = 0.0744 eV, 
7714 = 0.96 eV, ms = 5 keV, tjiq = 10 GeV, and 
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(4.11) 



"fit" all oscillation data quite well. 

Note that a null result from MiniBooNE is bound to place significant limits on the 
seesaw energy scale. If all right-handed neutrino masses are similar, the effective mixing 
angle that governs — > z/g transitions is sin^ 26'MiniBooNE ^ 4m^/M^. Hence, a null result 
at MiniBooNE would rule out a seesaw energy scale M lighter than 6 eV, assuming all 
active neutrino masses m are around 0.1 eV [71] . This limit is sensitive to the lightest 
neutrino mass mi and can be somewhat relaxed (similar to how I obtain a good 3+2 to 
all neutrino data) by postulating a (mild) hierarchy of right-handed neutrino masses and 
by assuming that sterile-electron and sterile-muon neutrino mixing is suppressed with 
respect to naive expectations for the lightest mostly sterile state(s). For larger values of 
mi, M values around 10 eV are already constrained by — >■ z/g searches at the NuTeV 



and NOMAD |87] experiments, and 



Vr searches at CHORUS 



142 



4.2.2. Cosmological and Astrophysical Constraints, Warm Dark Matter 



Very light sterile neutrinos that mix with the active neutrinos are constrained by several 
cosmological and astrophysical observables. The "seesaw" right-handed neutrinos are no 
exception. Given that active-sterile mixing angles l^/mp < m;/mj (a = e, fi,T, i = 4, 5, 6), 
it turns out that for "standard cosmology," the right-handed neutrinos thermalize with 
the early universe thermal bath of SM particles, as long as the reheat temperature is 
higher than their Majorana masses. For the low seesaw energy scales I am interested in, 
this is a problem. For the values of M/j under consideration here, thermal right-handed 
neutrinos easily overdose the universe. Smaller mi values {mi < 10^^ eV) lead to the 
possibility that right-handed neutrinos are the dark matter, as recently discussed in the 
literature |31l I73j . I comment on this and other possibilities shortly. 

Fig. 14.21 depicts the region of the |?7e6p ^ mg-plane in which the contribution of the 
heaviest neutrino z/g to Q (the normalized energy density of the universe, p/pc) is larger 
than 0.3 (dark region). The same constraint roughly applies for all a = e, /i,r and 
i = 4,5,6. The dashed diagonal lines correspond to If/^eP = ^i/^e, for different values 
of mi. All lines lie deep within the dark Qg > 0.3 region. 

For smaller values of M/j, the situation is also constrained. For values below 
tens of eV, thermal sterile neutrinos contribute to the amount of hot dark matter in 
the universe [901 I21L I22L I20L I91j . Right-handed neutrinos will thermalize as long as 
misin^Oia ^ 5 x 10""^ eV |92L 193] . In low energy seesaws, this roughly translate into 
m > 10^'^ eV, where m is the active neutrino mass scale. For the cases of interest here 
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Figure 4.2. Adapted from [89j . Cosmological and astrophysical constraints 
on the |f/a6p X mg-plane. In the large dark grey region, the density of a 
thermal i^q population is Qg > 0.3, while the light grey 'X-ray' region is 
disfavored by X-ray observations. The regions labeled 1,2,3 are preferred 
if one is to explain anomalous pulsar kicks with active-sterile oscillations 
inside supernovae. Regions 1 and 3 qualitatively extend inside the fi^ > 0.3 
part of the plane as indicated by the horizontal dotted and dash-dotted 
lines, respectively. The regions 'Warm Dark Matter" and "Too Warm Dark 
Matter" are meant to represent the region of parameter space where thermal 
i/g qualifies as a good (or bad) warm dark matter candidate. The region 
above the solid diagonal line is disfavored by the observation of electron 
(anti)neutrinos from SN1987A. The diagonal dashed lines correspond to 
U^Q = mi/rriQ, for different values of m;. Also shown is our "best fit" sterile 
solution for different pulsar kick scenarios, assuming the 3 + 2 LSND fit for 
the lighter states. The regions one and three best fit values are represented 
by circles and a star respectively. See text for details. 
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m > y^Am^ ~ 10 ^ eV, in which case rrii values above somewhere betweero 0.2-2 eV are 
ruled out [901 [941 [211 [221 [201 m]. Note that in our "best fit" 3+2 solution to the LSND 
data, the sum of all active neutrino masses violates slightly some of these constraints. 
This problem can be easily evaded if I choose mi values close to the lower end of the 
"allowed region". Other "mild" non-standard cosmology effects (see, for example, [95] ) 
are also known to alleviate the hot dark matter bound on neutrino masses. Note that 
any sterile neutrino solution to the LSND anomaly faces a similar problem, which must 
be resolved with non-standard cosmology. 

Big-bang nucleosynthesis also proves to be a large obstacle when it comes to the 
existence of light sterile neutrinos in thermal equilibrium in the early Universe (at tem- 
peratures above several MeV). In the absence of "non-standard" assumptions, big-bang 
nucleosynthesis constrains the existence of new thermal relativistic degrees of freedom 
(see, for example, |91L 196]). 

One way to avoid the bounds described above (see also, for example, |97] ) is to con- 
sider that the reheating temperature of the universe is very low. This way, right-handed 
neutrinos, in spite of their "large" mixing angles, never reach thermal equilibrium in the 
early universe and neither overdose the universe nor contribute to the amount of hot dark 
matter. Quantitatively, < 5 MeV is sufficient to avoid eV-mass (or heavier) sterile neu- 
trinos that are allowed to explain the LSND anomaly [98] . Unless otherwise noted, this is 
the assumption I make here. Other possibilities include adding new neutrino interactions 



The upper bound on the sum of neutrino masses from cosmology depends on several assumptions that 
go into analyzing the different cosmological observations. These include the issue of defining the values 
of the concordance cosmological model and deciding which data sets to include in the fit. 
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to lighter scalars, so that neutrinos remain in thermal equilibrium until they are non- 
relativistic [95J. According to [95] . this can even be accomplished for light neutrinos, as 
long as the neutrino-scalar field coupling is finely tuned (see, however, j91j ). Yet another 
possibility is to consider that the lepton asymmetry of the universe is large. The authors 
of |99j have recently studied this issue in detailed, and concluded that a lepton asymmetry 
of order 10~^ is required in order to allow the existence of LSND sterile neutrinos to be 
in agreement with data from large-scale structure and big-bang nucleosynthesis (see also 

Under these circumstances, it is interesting to consider whether light seesaw right- 
handed neutrinos still qualify as good warm dark matter. This could happen if their 
production in the early universe was non-thermal. One concern surrounding warm dark 
matter is whether it is ruled out by large scale structure surveys. Here, I will not add 
to this discussion but refer readers to the recent literature on the subject jlOOj . A brief 
summary of the situation is as follows: constraints on warm dark matter can be translated 
into a lower bound on the mostly sterile neutrino mass. The lower bound has been 
computed by different groups, and lies somewhere between 3 and 14 keV |100j . Different 
lower bounds depend on several issues, including which subset of Lyman alpha-forest data 
is taken into account. 

Another constraint on potential dark matter sterile neutrinos comes from the obser- 
vation of X-rays originating in galactic clusters. Such regions of the universe should be 
overdense with warm dark matter heavy neutrinos, which can be directly observed via 
their radiative decay Ui + j |101] . Bounds from X-ray observations have been sum- 
marized very recently in |102j . Combining the results of [102] and Fig. 14.21 I find that 
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for lightest neutrino masses larger than 10~^ eV, such bounds can only be avoided for 
rrii < 100 eV, where large scale structure constraints on warm/hot dark matter are severe. 
This qualitative analysis indicates that seesaw sterile neutrinos cannot, simultaneously, 
fit the LSND data and serve as cold dark matter. 

On the astrophysics side, the most severe constraint on light, sterile neutrinos is pro- 
vided by the observation of electron (anti) neutrinos coming from SN1987A. The current 
analysis consists of comparing the model-dependent neutrino flux at the surface of the 
neutrinosphere with that detected on Earth. Large sterile neutrino mixing and mass would 
result in modification/depletion of the detected neutrino signal (for a recent detailed dis- 
cussion, see |103] ). Although only twenty neutrinos were observed in this event, one can 
still place bounds on sterile-active neutrino mixing. As far as "LSND" sterile neutrinos 
are concerned, these bounds are still weaker than those obtained by the null short baseline 
oscillation experiments [104] and therefore already accounted for in our analysis. Heavier 
right-handed neutrinos can, however, be excluded by SN1987A neutrino data. Fig. 14.21 
depicts the region of parameter space excluded by SN1987A data (region to the right of 
solid, diagonal line). This bound is defined by mi\/2\f\T^i >0.22 keV |105L 198] . See also 
|93] . According to Fig. 14.21 supernova bounds force the seesaw scale to be below a few 
keV for m; values above 0.01 eV. 

4.2.3. Pulsar Kicks 

Pulsars are born from the gravitational collapse of the iron core of a massive star. These 
core collapse supernova are an excellent source of neutrinos, producing all (active) flavors 
copiously (see |106] for a detailed review). Current observations point to the fact that 
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some pulsars move with peculiar velocities much greater than those expected from an 
asymmetric supernova explosion mechanism. Quantitatively, current three dimensional 
models yield velocities up to 200 km/s |107] while pulsars moving at speeds as high as 
1600 km/s have been observed. I note, however, that some two-dimensional hydrodynamic 
studies |108j indicate that natural anisotropics generated during supernova explosions 
can, in fact, yield large neutron star velocities consistent with observations. However, 
more simulations seem to be required in order to validate this claim. Here, I will operate 
under the hypothesis that new physics, usually in the form of new neutrino physics, is 
responsible for the large pulsar kicks. 

Since roughly 99% of the approximately 10^^ ergs of energy released in a core collapse 
supernova is in the form of neutrinos, it is reasonable that neutrino physics provides a so- 
lution to this anomaly. At these rates a small (1—3)% asymmetry in neutrino emission can 
account for the observed large pulsar velocities. Neutrinos are always produced asymmet- 
rically in the polarized medium of the proto-neutron star, due to the left-handed nature 
of their interactions. Unfortunately, asymmetric production cannot solve this problem be- 
cause the associated medium densities are such that neutrinos undergo multiple scattering 
within the star's interior, eventually diffusing out of an effective surface, called the neu- 
trinosphere, with all initial asymmetries washed away. Several distinct mechanisms have 
been formulated to sidestep this fact. Specifically, the existence of large neutrino mag- 
netic moments has been explored in |109j . and can be tested in next generation neutrino 
scattering experiments [llOL IlllL I112L 1113] . Proposed solutions also exist which utilize 



148 



standard three flavor neutrino oscillations, where and Ur appearing between their neu- 
trinosphere and the larger v^, neutrinosphere can stream unhindered out of the star |114] . 
This solution is, however, currently disfavored by terrestrial oscillation experiments. 

I concentrate on the case of oscillations into sterile neutrinos, which can proceed in 
various ways, depending on the mass and coupling of the relevant neutrinos as well as 
the properties of the collapsing star, including its density and magnetic field. Following 
[89], I separate and analyze these within three distinct categories. Each one requires the 
existence of a keV-scale sterile neutrino with very small couplings to the active flavors, of 
the order 10~^ — 10~^, especially if light sterile neutrinos are thermally produced in the 
early universe. Under these circumstances, if seesaw neutrinos are to play the role of the 
sterile neutrinos responsible for pulsar kicks, \Uai\^ < m/mj {i = 4,5,6) must lie in the 
10~^ range for m, ~ 10^ eV. This implies m ~ 10"^ eV and is only compatible with a 
hierarchical active neutrino mass spectrum and very light mi, as identified in |31L 173] . 

Here, instead, I will concentrate on identifying solutions that will address pulsar kicks 
and the LSND anomaly. According to the discussion in the previous subsection, the 
mostly sterile neutrino masses m4 and are constrained to be less than 10 eV so that 
a 3+2 solution to the LSND anomaly can be obtained from the seesaw Lagrangian. The 
heaviest neutrino mass mg is unconstrained, so I are free to vary it as needed in order to 
attack the pulsar peculiar velocity issueo Naively, the fraction of {a = e, /x, r) in pq is 



expected to be of the order Uas ~ a/O.3 eV/3 x 10'^ eV = 10~^, much too large to satisfy 
the pulsar kick plus cosmology constraints summarized in Fig. 14. 2[ On the other hand, 
once the fls < 0.3 constraint is removed, the 'pulsar kicks' allowed region of the plane 



can neglect the lighter sterile neutrinos (1/4 and 1^5) as they should not alter the kicking mechanism 
significantly due to their small mass and non-resonant production. 
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is significantly enlarged, as qualitatively indicated by the horizontal lines in Fig. 14.21 In 
this case, which I must consider anyway if I am to have agreement between LSND and 
searches for hot dark matter, one can envision explaining pulsar kicks and the LSND data 
simultaneously. Note that once heavy sterile neutrinos are "removed" (so that they do 
not overdose the universe), constraints from X-ray observations (see Fig. 14. 2p are also 
removed. 

In scenario 1, the pulsar kick is produced via an active-sterile MSW resonance in the 
core of the proto-neutron star at large densities, greater than 10^^ g/cm^, and magnetic 
fields, near 10^^ G |115] . The effective neutrino matter potential in material polarized 
by a strong magnetic field contains a term proportional to k ■ B /\k\ [116L I117j . where k 
is the neutrino's three-momentum and B is the local magnetic field vector. Clearly the 
MSW resonance occurs at a radius that depends on k ■ B/\k\, the relative orientation of 
the neutrino momentum and magnetic field. Sterile neutrinos produced at smaller radii 
(higher temperatures) carry greater average momentum than those produced at larger 
radii (lower temperatures), yielding an asymmetric momentum distribution of emitted 
neutrinos. This asymmetry is capable of producing the observed pulsar kicks, in the 
direction of the magnetic field, when the mass and coupling of the sterile state is near 
8 keV and above 1.5 x 10"^, respectively [89]. I found the "best fit" to the LSND data 
(using z/4 and v^) and pulsar kicks (using z/g) and mg > 5 keV. The |?7a6p and mg "best 
fit" values are depicted in Fig. 14. 2[ This solution is strongly disfavored by the observation 
of neutrinos from SN1987A. The fact that If/gel is much larger than the other two active- 
sterile mixing angles is due to the fact that \Uei\ and |[/e5| are constrained by LSND data 
to be much smaller than naive expectations (see Eq. (14.91) ). In order to reduce |t/eg|, one 
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would have to either reduce m; by an order of magnitude - which renders the 3+2 fit to 
oscillation data very poor - or increase mg, which would only push \Uee\ deeper into the 
region of parameter space ruled out by SN1987A. One can however, find 3+1 solutions to 
LSND data where could pose as the sterile neutrino that explains why pulsar peculiar 
velocities are so large (see Eqs. (I4.10|4.1ip ). 

Scenario 2 also relies on a direction-dependent MSW resonance, this time occurring 
outside the core where the matter density and temperature are much lower. Here, both 
the active and sterile neutrinos are free to stream out of the star. The departing active 
flavors still have a small interaction cross- sect ion, cr ~ G'^E^, and can therefore deposit 
energy and momentum into the star's gravitationally bound envelope proportional to the 
matter it transverses. Via the direction-dependent resonance, neutrinos moving in the 
direction of the magnetic field remain active longer, deposit more momentum, and thus 
kick the star forward. The observations can be explained in this case with a smaller sterile 
neutrino mass and larger active-sterile coupling near 4 keV and 4.5 x 10~^ respectively 
|89j . In the case of our LSND "fit" to the data, I can constrain one of |?7^6p ot |f/r6p 
to lie inside region 2. The other |f/a6p (a = e, r or e,fi), however, are constrained to be 
large, thus violating the SN1987A bound in much the same way as the scenario 1 best fit 
results. Another possibility is to choose all \Ua6\^ of the same order of magnitude. I do 
not explicitly consider these points as they reside in the region of parameter space where 
region 2 and 3 overlap, and behave in the same way as the point described below, under 
scenario 3. 

Scenario 3 proceeds through off-resonance production of the sterile neutrino in the 
proto-neutron star core [118] . The amplitude for sterile neutrino production by a weak 



151 



process is proportional to U^q, the effective mixing angle between the heavy mass eigen- 
state and the ffavor eigenstate. Initially, this quantity is very small due to matter effects 
in the dense core. The effective potential in the star's interior is quickly driven to zero 
in the presence of sterile neutrino production by a negative feedback mechanism. If this 
occurs in a time less than the diffusion time-scale for the active neutrinos, approximately 
(3 — 10) s, the mixing angle will reach its vacuum value [119] . The sterile neutrinos will 
then be produced and emitted asymmetrically and thus kick the pulsar to large velocities. 
Lower limits on the vacuum mass and mixing values are derived by requiring that the 
off-resonance time scale (inversely proportional to mg sin^ 29^^) for the evolution of the 
matter potential to zero be less than about ten seconds. This places the sterile mass 
and mixing at approximately 1 keV and above 5 x 10~^, respectively. Since all three 
active flavors are present in equal abundances, and all contribute to the effective matter 
potential, the mixing angle in question is not any particular [/^e- Rather it is the angle, 
6q associated with the projection of z/g onto the space spanned by i^e, I'n and i^r, that is 
6l = f/gQ + f/^g + U^Q. From Eq. ( 14. 6p I see that 6e = \/rail mg up to corrections due to the 
non-unitarity of Y and active neutrino mass differences. This is independent of mixing 
angles, and therefore cannot be tuned to be small. Our "best flt" region-3 solution is de- 
picted in Fig. 14.21 by a star. It turns out that f/ag have very similar values for a = e, /z, r. 
In order to evade the SN1987A constraint, I were forced to pick mi values close to lower 
bound of our 3+2 LSND "flt" (m/ = 0.22 eV), so that \UeA\-, and |?7^5| are close to 

the low-end of the allowed range in Table 14.11 
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4.2.4. Supernova Nucleosynthesis 

Core collapse supernova are believed to produce the observed heavy element {A > 100) 
abundance through the r-process, or rapid neutron capture process. Here I briefly re- 
view this mechanism (see |120] for a comprehensive review), as well as its facilitation 
by the addition of active-sterile neutrino oscillations |121j . This scenario begins in the 
neutrino driven wind; that is, the wind of ejected nucleons driven by neutrinos radiated 
from the cooling proto-neutron star. The maintenance of equilibrium among neutrons, 
protons, and electron (anti) neutrinos in neutrino capture processes leads to a neutrino- 
rich environment. As the wind propagates, it cools enough for all free protons to bind 
into a particles. In the ideal r-process picture, as the wind cools further these a particles 
bind into intermediate size seed nuclei which later undergo neutron capture to form the 
observed heavy r-process elements. 

This ideal scenario is dampened by the large number of electron neutrinos present at 
the stage of a particle formation. These will capture on the free neutrons, converting 
them to protons, which in turn will fuse to make more a particles. The end result is a 
very small free neutron to a particle ratio, conditions unfavorable for r-process element 
formation. This is known as the a effect and must be circumvented to produce the correct 
distribution of heavy elements. A clear solution to this problem is to reduce the number of 
electron neutrinos present at this stage, which can be accomplished by resonant — > z/^ 
conversioio [121] . 



''in this mechanism the effective matter potential, which depends on the number of electrons, positrons 
and neutrinos, varies wildly as a function of distance from the core. Along this radial direction there are 
three relevant MSW resonant conversions that must be tracked and understood: ~^ '^s, I'e — * i^s and 
i^s i^e- See [121) for more information. 
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The sterile neutrino solution to the r-process mechanism is modeled and fit to the 
data in reference |122j including the effects of relevant nuclear physics and additional 
neutrino oscillations in the star's envelope. The analysis is expanded in [123j with the 
inclusion of fission cycling of the produced heavy elements. The analysis indicates the need 
for an eV-scale sterile neutrino with an allowed parameter space much larger than that 
constrained by LSND. By itself, the requirement of successful r-process in supernovae only 
weakly constrains the light neutrino mass scale to be greater than 10~^ eV and 10"^ eV 
for a 1 eV and 10 eV sterile neutrino, respectively. With regard to the LSND results, 
it has been demonstrated that the 3 + 1 oscillation scenario fits within this parameter 
space |123] . Considering that the best fit mass-squared difference and mixing angles for 
the fourth mass eigenstate, which makes up most of the lightest sterile neutrino, is very 
similar between the 3 + 1 and 3 + 2 case |85] , it is reasonable to conclude that z/g ^ Usa 
resonant conversion will also fit within this scenario. Even oscillations into the heavier 
z/s5 state can potentially solve this anomaly if the neutrino driven wind expansion time- 
scale is sufficiently small, < 0.1 sec. To conclude this section I note that, although the 
sterile neutrino solution to the supernova nucleosynthesis problem fits well within our 
seesaw framework, it adds no additional constraints, and therefore does not increase the 
predictability of our scenario. 

4.3. Other Probes of the Seesaw Energy Scale 

Here I survey other existing and future probes of light sterile neutrinos. As opposed 
to the previous cases, these probes are perfectly consistent by themselves. That is, extra 
heavy neutrinos are not required to solve problems within the system. However, their 
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addition can lead to large modifications to the outcome of such experiments, thus ren- 
dering the eV-scale seesaw scenario testable. Specifically, I consider bounds from tritium 
beta-decay and neutrinoless double-beta decay. Observations in all of these areas have 
already yielded useful constraints on sterile neutrinos, and the situation is expected to 
improve in the next few years. 

4.3.1. Tritium Beta-Decay 

The endpoint of the electron energy spectrum in the beta-decay of tritium is a powerful 
probe of nonzero neutrino masses. This results from the decay kinematics of the system 
which is necessarily modified by the presence of a massive neutrino. The nonzero neutrino 
mass effect can be understood almost entirely from the analysis of the phase space distri- 
bution of the emitted electrons, and is therefore quite model independent. Existing beta 
decay experiments extract limits on an effective electron neutrino mass m^^ = \ Uei\'^m'f 
[23], provided that the neutrino masses are smaller than the detector energy resolution. 
Currently the most stringent bounds on ml^ are (2.3)^ eV^ at 95% confidence from the 
Mainz experiment |17] and (2.5)^ eV^ at 95% confidence from the Troitsk experiment 
|18] . In the next few years the Katrin experiment should exceed these limits by nearly 
two orders of magnitude, probing down to (0.2)^ eV^ at the 90% confidence level |63j . One 
might naively compute this effective mass for the "best fit" mixing parameters obtained 
in the previous section. In this case, I expect a keV seesaw neutrino to contribute to m^^ 
by a huge amount Sml^ = Ul^ml ~ ^''^1 = "^/"^6 [32j so that it would be excluded by 
the current precision measures of tritium beta-decay for m/ > 10^'^ eV. This is clearly an 
incorrect treatment of the physics. As pointed out in, for example, [124] . the existence of 
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a heavy neutrino state would produce a kink in the electron energy spectrum of size |?7eiP 
at an energy Eq — rrii as well as a suppression of events at the endpoint of order 1 — |t/ejp. 
Here Uei (i = 4, 5, 6) is the mixing between the electron neutrino flavor eigenstate and the 
heavy mass eigenstate, while Eq = 18.6 keV is the endpoint energy of tritium beta-decay. 

Fig. 14.31 depicts 1 — S/Sq, where S is the /3-ray energy spectrum obtained assuming 
three mostly active, degenerate neutrinos with mass m = 0.1 eV and one mostly sterile 
neutrino Ui with various masses rrii and mixing angle U^^ = m/rrii, while 5*0 is the spectrum 
associated with massless neutrinos. One can readily observe "kinks" in the spectrum above 
rrii. For /3-energies above Eq — rrii, the impact of the sterile state is to "remove" around 
1 — |t/ejp of the /3-rays from spectrum. This is most significant between Eq — rrii and Eq 
minus the mass of the active neutrinos. For energies below Eq — rrii, the spectrum agrees 
with that obtained from the emission of one effective neutrino with mass-squared m^^. 

I estimate the sensitivity of future tritium beta-decay experiments to the emission of 
one heavy state by considering the ratio between the number of electrons with energies 
above Eq — AE in the case of one heavy massive neutrino Ui and in the case of massless 
neutrinos 



where dN/dE is the energy distribution of /3-rays, which depends on the neutrino mass rrij. 
This expression can be easily generalized for more than one heavy neutrino. The advantage 
of using the ratio above is that potential systematic uncertainties and normalization effects 
can be safely ignored. An experiment is sensitive to a massive neutrino state if it can 
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Figure 4.3. 1 — S/Sq as a function of the /?-ray energy, where S is the 
/?-ray energy spectrum obtained assuming three mostly active, degenerate 
neutrinos with mass m = 0.1 eV and one mostly sterile neutrino z/j with 
rrii = 0.1, 0.5, 1, 5, and 10 eV. The mixing angle is given by = m/rrii. Sq 
is the spectrum associated with massless neutrinos. See text for details. 

distinguish R from unity, a determination that should be limited by statistics due to the 
very low /3-ray flux in the high-energy tail of the electron spectrum. 

In order to compute R, I use an analytic expression for Eq. fl4.12p . which exists pro- 
vided that one neglects nucleon recoil in the decay. Fig. 14.41 depicts constant i?-contours 
in the \Uei\ x rrii plane. Contours were computed for AE = 25 eV, in order to allow one to 
easily compare our results with the sensitivity estimates of the Katrin experiment. After 
data-taking, Katrin is expected to measure R at the 0.1%-1% level (lightest grey region). 



Its sensitivity is expected to be \/fn^^ > 0.2 eV at the 90% confidence level. This can be 
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extracted from the plot by concentrating on the Uei = 1 hne. Note that while the expected 
energy resolution for Katrin is of order 1 eV, the expected number of signal events above 
Eq — 1 eV is minuscule (both in absolute terms and compared with expected number of 
background events), so that most of the sensitivity to nonzero neutrino masses comes from 
analyzing the shape of the electron spectrum in the last tens of electron-volts. A larger 
"window" would suffer from increased systematic uncertainties, so that AE ~ 25 eV is 
representative of Katrin's optimal reach [63J. 

The shape of the constant R contours is easy to understand. As already discussed, 
for rrii > AE, the effect of the right-handed neutrinos is to reduce the spectrum in an 
energy independent way by 1 — |t/e.jp, while for rrii < AE, states with the same effective 
mass-squared mf\Uei\'^ produce the same effect in tritium beta-decay so that the diagonal 
lines coincide with lines of constant mf\Uei\'^- 

A more sensitive approach would be to "bin" the last tens of eV of the "data" into 
1 eV bins, and fit the distribution to a massless neutrino hypothesis. For the values of 
the parameters in which I am interested, I find that one 25 eV bin yields roughly the 
same sensitivity to nonzero neutrino masses as twenty five 1 eV bins for large masses and 
small mixing angles. For smaller masses and larger mixing angles, a "binned" analysis 
should be sensitive to effects which are localized in individual bins (such as "kinks"). 
Another recent estimate of the sensitivity of tritium beta-decay experiments to heavy, 
sterile neutrinos can be found in [93j. Our results agree qualitatively. 

Fig. 14.41 also depicts the loose upper bound for Uei = ^f^ilf^i as a function rrii, 
for vTii = 0.32 eV and mi = 0.01 eV. For mi > 0.1 eV, Katrin should be sensitive to 
Mji < 1 keV while for mi > 0.01 eV (the solar mass scale) Katrin should be sensitive to 
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Figure 4.4. Contour plot of constant i?, as defined by Eq. ( 14.120 . assuming 
an energy window AE = 25 eV. The solid (dashed) line corresponds to 
^^/mi/rrii, a naive upper bound for \Uei\, for m/ = 0.3 eV (0.01 eV). The 
circles correspond to Uei for the three mostly sterile states obtained by our 
"fit" to other neutrino data, Eq. (14. 9p . See text for details. 



Mpi > 10 eV and Mn < 100 eV, where here I assume that all right-handed neutrino masses 
are of order Mn. In the case of seesaw parameters that fit the LSND data with a 3+2 
neutrino spectrum (see Eq. (14. 9p ). expectations are high as far as observing a kinematical 
neutrino mass effect at Katrin, in spite of the fact that the fit to LSND data requires 
UeA and Ue5 to be "abnormally" low. Fig. 14.41 depicts Uei and m, values for the heavy 
neutrinos (open circles). The contribution of the heaviest of the two LSND- related sterile 
neutrinos is of order the Katrin sensitivity, while the active contribution itself, which leads 
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to ml = 



X]i=l,2,3 l^ei^i 



mf is already within the Katrin sensitivity range, given that 



large rrii > 0.22 eV values are required by our 3+2 LSND "fit". The effect of z/g is small 
if tuq is larger than 1 keV (required if one takes the "pulsar kicks" hint into account), but 
would be very significant if mg were less than 1 keV. 

4.3.2. Neutrinoless Double-Beta Decay 

If the neutrinos are Majorana fermions - as predicted in the case of interest here - lepton 
number is no longer a conserved quantity. The best experimental probe of lepton number 
violation is the rate for neutrinoless double-beta decay. This process, which violates lepton 
number by two units, is currently the subject of intense search |125[ , 140] . If neutrino 
masses are the only source of lepton number violation, the decay width for neutrinoless 
double-beta decay is 



where Ai is the relevant nuclear matrix element and ~ 50^ MeV^ is the relevant 
momentum transfer. In the limit of very small neutrino masses {mf <^ Q^), ^oujsp is 
proportional to an effective neutrino mass \mee\, 



The sum is over all light neutrino mass eigenstates. In the case of a low-energy seesaw, 
when all mj , i = 1 , . . . , 6 are much smaller than , it is easy to see that rriee vanishes |32] . 
The reason for this is that, in the weak basis I are working on (diagonal charged- lepton 
and charged weak- current), rriee is the ee-element of the neutrino mass matrix, as defined 




2 



(4.13) 



n 




(4.14) 
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in Eq. (14.21) . One can trivially check that, by assumption, not only does rriee vanish, but 
so do all other m^/j, a, P = e, /z, r. Note that this result does not depend on the fact that 
I have been assuming all elements of the neutrino mass matrix to be real jl26j . 

For heavy z/j neutrinos, U^iirii no longer captures the dependency of Foy/?/? on the ex- 
change of Ui. For mf ^ Q^, instead, the dependency on neutrino exchange is proportional 
to U^Jnii. If this is the case, the overall contribution (including all heavy and light states) 
is no longer proportional to rriee but, instead, can be qualitatively expressed as a function 
of an effective m't^, 



The approximation Fojy/3/3 oc \ml^\ is good as long as one can neglect the dependency of 
Ai on rrii and is not expected to be a great approximation when mj ~ Q^. Nonetheless, 
rrigf still qualitatively captures the behavior of Fojy/3/3 as a function of the sterile neutrino 
masses and studying its behavior is sufficient for our ambitions in this discussion. 

Fig. 14. 51 depicts m^f for our "best fit" 3 + 2 LSND solution (see Sec J4.2.1[) . as a function 
of the unconstrained mg. As advertised, rrigf vanishes for mg <^ Q^. The figure also 
depicts the "active only" value of m^f"^ = ^,=1 2 3 Even in the limit ml ^ Q^, 

there is still partial cancellation between the mostly active and mostly sterile "LSND" 
states. This is a feature of the Lagrangian I am exploring here, and is not in general 
observed in other scenarios with light sterile neutrinos tailor-made to solve the LSND 
anomaly. 

Currently, the most stringent limits on this effective mass comes from the Heidelberg- 
Moscow experiment |127j where they find rriee < 0.91 eV at 99% confidence. In the near 
future, experiments aim to reach down to ruee values close to 10~^ eV |125[ , 140] . A signal 
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Figure 4.5. Effective rriee for neutrinoless double-beta decay as a function of 
me, tlie lieaviest riglit-lianded neutrino mass, assuming tlie existence of only 
light, active, neutrinos (magenta curve), with a degenerate mass spectra, 
and for our "best fit" 3 + 2 LSND sterile neutrino solution (blue curve). 
See text for details. Also indicated is the parameter region preferred by 
astrophysical hints of sterile neutrinos. I assume Q = 50 MeV. In the case 
of a low-energy seesaw, uiee vanishes as long as mg <^ Q. 

would rule out a seesaw scale below tens of MeV. On the other hand, if I were to conclude 
that the neutrino masses are quasi-degenerate (through, say, a signal in tritium beta- 
decay) and if the LSND 3+2 solution were experimentally confirmed, a vanishing result 
for could be considered strong evidence for a very small seesaw scale. On the other 
hand, if this were the case {mee zero, large active neutrino masses), it would also be very 
reasonable to conclude that neutrinos are Dirac fermions. Distinguishing between the two 
possibihties would prove very challenging indeed. It is curious (but unfortunate) that in a 
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low energy seesaw model the neutrinos are Majorana fermions, but all "standard" lepton- 
number violating observables vanish, given that their rates are all effectively proportional 

to TTlap^. 

4.4. Conclusions 

The "New Standard Model" (equal to the "old" Standard Model plus the addition of 
three gauge singlet Weyl fermions) is, arguably, the simplest extension of the SM capa- 
ble of accommodating neutrino masses. This Lagrangian contains a new dimension-full 
parameter: the right-handed neutrino mass scale M, which must be determined experi- 
mentally. Unlike the Higgs mass-squared parameter, all M values are technically natural 
given that the global symmetry of the Lagrangian is enhanced in the limit M —>■ 0. 

Very large M values are theoretically very intriguing, and have received most of the 
attention of the particle physics community. There are several strong hints that new 
phenomena are expected at the electroweak breaking scale ~ 10^ GeV, the grand unified 
scale ~ 10^^^^^ GeV, and the Planck scale ~ ig^*^"^^ GeV, and it is tempting to associate 
M to one of these energy scales. Furthermore, large M values provide an elegant mech- 
anism for generating the matter antimatter asymmetry of the Universe [68]. Of course, 
large M values are experimentally very frustrating. It may ultimately prove impossible 
to experimentally verify whether the New Standard Model is really the correct way to 
describe Nature. 

Here, I explore the opposite end of the M spectrum, M < 1 keV. Such values are phe- 
nomenologically very intriguing, given that small M values imply the existence of light 
sterile neutrinos that mix significantly with the active neutrinos and can potentially be 
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directly observed. Furthermore, there are several experimental and astrophysical phenom- 
ena that are best understood if one postulates the existence of light, moderately mixed 
sterile neutrinos. I find that by requiring all three right-handed neutrino masses to be less 
than a few keV I can simultaneously explain all neutrino oscillation data, including those 
from LSND, explain the large peculiar velocities of pulsars, and accommodate the produc- 
tion of heavy elements in supernova environments. Our fit also provides constraints for 
the "active" neutrino oscillation parameters, most strongly to the lightest active neutrino 
mass. All successful parameter choices that accommodate the LSND data require mi to 
be "large" {mi > 0.1 eV), and the "best fit" requires all active neutrino masses to be 
quasi-degenerate. It is important to emphasize that the presence of light sterile neutrinos 
that mix relatively strongly with the active neutrinos is only in agreement with cosmolog- 
ical data (especially large scale structure and big-bang nucleosynthesis) if non-standard 
cosmological ingredients are present. 

Fig. 14.11 depicts such a scenario. This six neutrino mass spectrum (including mixing 
angles) fits all neutrino oscillation data (including those from LSND), provides a sterile 
neutrino solution to the pulsar kick puzzle, and contains all the necessary ingredients 
for heavy element nucleosynthesis in supernovae. The heaviest of the neutrinos does not 
qualify as thermal warm dark matter (in the absence of new cosmological ingredients, its 
presence would overdose the Universe). Note that even if it were non-thermally produced, 
constraints from the observation of X-rays from the center of the galaxy would rule out 
as a good dark matter candidate. Lighter z/g masses will evade X-ray constraints, but 
would render i/g too "hot," and hence not a good dark matter candidate. 
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On the negative side, low M values are, theoretically, rather puzzling. In order to 
obtain the observed light neutrino masses, neutrino Yukawa couplings are required to be 
much smaller than the electron Yukawa coupling, and it is tempting to believe that such 
small numbers are proof that a more satisfying understanding of fermion masses must 
exist. Furhermore, thermal leptogenesis is no longer an option (see, however, |128j ). 
Finally, the fact that M is naively unrelated to other mass scales can also be perceived as 
disheartening, but, in our opinion, should be interpreted as evidence that there is more 
to the lepton sector than meets the eye. 

Regardless of one's preference for a high or low seesaw energy scale, and independent 
of whether the data from LSND and the astronomical observables discussed above have 
anything to do with sterile neutrinos, our main point is that the determination of M 
is an experimental issue. In the near/intermediate future, low energy seesaw scales will 
be probed by several experiments, most importantly measurements of the end-point of 
tritium beta-decay, the MiniBooNE experiment, searches for neutrinoless double-beta 
decat and, if we get lucky, the detection of neutrinos from a nearby supernova explosion. I 
find, for example, that Katrin should be sensitive to seesaw energy scales below tens of keV 
if all right-handed neutrino masses are similar, while null results from MiniBooNE would 
severely constrain right-handed neutrino masses below several eV. I conclude by pointing 
out that larger (but still "small") values of M are much harder to constrain. For GeV 
sterile neutrinos, typical active-sterile mixing angles are f/^j < 10^^°, probably too small 
to observe in particle physics processes. It is frustrating (and, I hope, ultimately false) that 
we seem to be unable to experimentally distinguish M ~ 1 GeV from M ~ 10^^ GeV. . . 
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APPENDIX A 



Standard Model of Particle Physics 



The Standard Model (SM) is a local Lorentz invariant quantum field theory defined 
by its gauge symmetry, field content and symmetry breaking scale. The gauge symmetry 
is 



Each group has its own distinct gauge coupling and is necessarily associated with spin one 
gauge boson fields transforming under the adjoint group representations. These are the 
eight gluons (G), three weak isobosons (W), and single hypercharge boson (B), respec- 
tively. The fermion field content, along with its associated quantum numbers, is listed in 
Table \M Here, i is a fiavor label that runs over the three known generations. In other 
words, there are at least three copies of each fermion generation identical in interactions, 
which are only mass scale. Additionally, there is a single complex Lorentz scalar, the 
Higgs boson H, with quantum numbers (1, 2, —1/2). 

From here it is a simple matter to write down the most general Lagrangian incorpo- 
rating these principles. This consists of three parts; the gauge sector 



Ccauge = L^l ^.D^' U + Qil ^.D''Qi + el-i ^D^' ei + da ^D^' di + ul^i ^D^'ui + \D>' H\\ (A.2) 



SU{?>)c X SU{2)l X f/(l) 
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Field 


SU{3), SU{2)l U{l)y 


el 

Q^ 

Wi 
di 


1 2 -i 

1 1 1 

3 2 i 
3 1 -| 
3 1 i 



Table A.l. Summary table of the Standard Model fermion field content 
along with their associated gauge quantum numbers. The fields are written 
in terms of their chiral left handed projections which may be conjugated 
to yield the right handed projections. The boldface entries describe the 
representaion of the field under the nonabelian symmetry groups. 



the Yukawa sector 

^Yukawa = v'/ LiHcj + yjQiHdj + y'^QiH^Uj + h.cU, (A. 3) 

and the Higgs sector 

jCmggs = f^'H^H + r]{H^Hf. (A.4) 

Here, = d^ — igsG^\ — igWfj,T — ig'BfjY is the covariant derivative containing the gauge 
fields with the matricies A and r, the generators of SU{3)c and SU{2)i in the appropriate 
representation. Notice that all fermions and gauge bosons are strictly massless in this 
theory. This follows from the chiral nature of the fermions where the left and right 
handed fields transform differently under SU{2)l x f/(l)y thus forbidding the required 
ipL'^R mass term. 



The simple multiplicative notation employed here is schematic. Since all terms are weak isosinglets it is 
clear that some spinor products between doublets are nontrivial involving the antisymmetric tensor eaf3- 
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The standard treatment now postulates that the /i^ < in Eq. (]A.4p so that H acquires 
a vacuum expectation value (vev). This breaks electroweak symmetry 

SU{2)l X U{1)y ^ U{1)q (A.5) 

at the energy scale characterized by v. Once done, only one abelian gauge symmetry 
remains to mediate electromagnetic interactions. Upon expanding the SM Lagrangian 
about V, one finds that the charged fermions (now under Q) have acquired masses as well 
as three of the electroweak gauge bosons. These extra longitudinal degrees of freedom 
arise from the Goldstone Higgs bosons "eaten" by the gauge fields. One can easily rewrite 
the Lagrangian terms expanded about the new vacuum and in the basis of definite gauge 
boson mass, but this is not needed for what follows. For completeness, I do point out that 
orthogonal linear combinations of Wi and W2 combine to form the massive bosons, 
and similarly, combinations of W3 and B combine to form the massive Z and (photon) 
A bosons. The 5*0(2) transformation that facilitates the rotation to this mass basis is 
defined by 6^^, the weak mixing angle. One should also notice that the massless photon 
is now the field associated with the remaining unbroken symmetry. 

Prom the point of view of this thesis, the most important feature of the SM Lagrangian 
is its accidental symmetries. The gauge and Higgs sectors contain a huge global symmetry 
S'[/(3)j corresponding to the freedom to rotate each of the five field families in generation 
space. The Yukawa terms substantially break this but still leave some residual symmetries. 
In particular, I find that if I assign identical U(l) charges, called baryon number (B), to 
Qi, Ui and di while holding the rest of the field content neutral, B is conserved. Similarly, 
I see that Lepton number (L) is conserved provided that I assign identical charges to 
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both Li and Cj with the remaining content neutrao This is not the whole story, as both 
B and L are broken by nonperturbative instaton effects. At the end of the day, only 
one non-anomalous symmetry remains, conserving the quantity B — L, U{1)b-l- It is 
important to understand that this symmetry is completely accidental within the SM and 
that any additional new physics is likely to break it. This implies that searches for B 
and L violation are ideal laboratories for studying physics BSM. The search for Lepton 
Number Violation (LNV) is the major unifying theme of this work, as it is intimately 
related to the Majorana nature of neutrinos. 



^Of course, in the SM with no means of generating neutrino mass, a separately conserved lepton number 
may be assigned to each lepton family. This is broken by the observed neutrino mixing phenomena 
induced by BSM physics 
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APPENDIX B 



Extended SU{2) neutrino mixing matrix notation 

Neutrino mixing and all of its invariance may be expressed as products of three distinct 
matrix classes in addition to the identity I □• Two are defined in terms of continuous 
variables. These are real orthogonal rotations in the a — b plane R'^''(^^) and single phase 
rotations P'*(0) defined by 



and 



ab 



ab 



= [R'^''(0)]fefe = cos^ 
'[R"^^)]ba = -sin0 



S,j (ij ^ a,b) 



(B.i: 



(B.2) 



respectively. The discrete transformation 



^ab _ ^ba _ l^ab^^) {P^(vr) - P"(vr)} 



(B.3) 



completes the set. It is a simple matter to write an arbitrary mixing matrix with such com- 
ponents. For example, the complex rotation matrix typically assigned to the 1 — 3 plane in 



"'^It is also true that one may use the full SU{n) algebra to parameterize n neutrino mixing in a more 
elegant way than in my extended SU{2) parameterization. These algebras are well known, but their 
use would require new commutator relations for each separate case, which is not conducive to my global 
approach. 
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standard three neutrino mixing analysis is R^''(^, S) = P^{-S)P^{S)R''^{e)P''{S)P^{-S). 
Additionally, discrete permutations of ab vector elements are accomplished with S^^ = 
R^''(7r/2)P''(7r). Manipulations of these matrices utilize their commutation relations de- 
fined by 



[R^''(^),R^'=(^')] = sin^(l-cos^')A^'' + sm^'(l-cos^)A^'= (B.4) 

+isin^sin^'{R-(-|)-R«'=(^)} 

[R^^{9),P^{(f))] = sin^(e^'^- 1) A^^ (B.5) 

[R^*'(^),A'^*'] = -sin^{P'^(7r) -P^(7r)} (B.6) 

[R"''(^),A"'^] = sin^A'''= + i(cos^- 1) |r"'=(-|) -R"'=(|)} (B.7) 

[A^P-(0)] = l(e^<^-l)|R^Vf)-R^'(^)} (B.8) 

[Aab^A-] = i{R^'^(-|)-R^'=(|)}. (B.9) 

All other commutators vanish. This matrix set closes upon itself facilitating a simple 
platform to perform manipulations. In that spirit, the following identities are easy to 



prove and useful in calculations: 



p-(0)p-(-0) 


= I 


(B.IO) 




= R'^'^(^)R^^(^) = / 


(B.ll) 


R"''(^ + 9') 


= R^^(^)R^''(^') 


(B.12) 


P"(0 + 0') 


= P^(0)P^(0') 


(B.13) 


P^(7r)P^7r) 


= R^'^(7r) 


(B.14) 


P^(7r)R"''(^) 


= R"''(-^)P''(7r) 


(B.15) 


R^^(^)P"(0) - P^(0)R"^(^) 


= cos^(e^'^- 1) A^^ 


(B.16) 


A^^P^(0) 




(B.17) 


R^''(^)R'=^(-) 
2 


= R'=^(-)R'^*'(^). 
2 


(B.18) 



